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FOREWORD 

A  short  history  of  the  development  of  the  prediction  methods  in  this  Technical 

Note  will  permit  the  reader  to  compare  them  with  earlier  procedures.     Some  of  these 

miethods  were  first  reported  by  Norton,    Rice  and  Vogler  [  1955],     Further  development  of 

forward  scatter  predictions  and  a  better  understanding  of  the  refractive  index  structure  of 

the  atmosphere  led  to  changes  reported  in  an  early  unpublished  NBS  report  and  in  NBS 

Technical  Note  15  [Rice,    Longley  and  Norton,    1959].     The  methods  of  Technical  Note  15 

served  as  a  basis  for  part  of  another  unpublished  NBS  report  which  was  incorporated  in 

Air  Force  Technical  Order  T.  O.    31Z-10-1     in  1961.     A  preliminary  draft  of  the  current 

technical  note  was  submitted  as  a  U.S.   Studv  Group  V  contribution  to  the  CGIR  in  1962, 

Technical  Note  101  uses  the  metric  system  throughout.     For  most  computations 

both  a  graphical  method  and  formulas  suitable  for  a  digital  computer  are  presented.     These 
include  simple  and  comprehensive  formulas  for  computing  diffraction  over  smooth  earth  and 
over  irregular  terrain,    as  well  as  methods  for  estimating  diffraction  over  an  isolated  rounded 
obstacle.     New  ennpirical  graphs  are  included  for  estimating  long-term  variability  for  sev- 
eral climatic  regions,    based  on  data  that  have  been  made  available  to  NBS. 

For  paths    in  a  continental  temperate  climate,    these  predictions  are  practically 
the  same  as  those  pu"blished  in  1961.     The  reader  will  find  that  a  number  of  graphs  have  been 
simplified  and  that  many  of  the  calculations  are  more  readily  adaptable  to  computer  program- 
ming.    The  new  material  on  time  availability  and  service  probability  in  several  climatic  re- 
gions    should  prove  valuable  for  areas  other  than  the  U.S.A. 


Note;     This  Technical  Note  consists  of  two  voluines  as  indicated  in  the  Table  of  Contents. 
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Annex    I 

AVAILABLE  DATA  AND  STANDARD  CURVES 
The  simplest  way  to  predict  long-term  median  transmission  loss  values  would  be  to 
use  a  best-fit  curve  drawn  through  measured  data  (represented  by  their  overall  median  values), 
plotted  as  a  function  of  path  length.     Such  a  method  ignor';s  essentially  all  of  our  understanding 
of  the  physics  of  tropospheric  propagation,    is   subject  to  especially  large  errors  over  rough 
terrain,    and  such  empirical  curves  represent  only  the  conditions  for  which  data  are  available. 

Curves  that  may  be  useful  for  establishing  preliminary  allocation  plans  are  presented 
in  section  I.  2  of  this  annex.      These  "standard"  curves  were  prepared  for  a  fixed  combination 
of  antenna  heights  and  assume  propagation  over  a  smooth  earth.      The  curves  are  not  suitable 
for  use  on  particular  point-to-point  paths,    since  they  make  no  allowance  for  the  wide  range 
of  propagation  path  profiles  or  atmospheric  conditions  that  may  be  encountered  over  particular 
paths. 

I.  1    Available  Data  as  a  Function  of  Path  Length 

Period-of-record  median  values  of  attenuation  relative  to  free  space  are  plotted  vs. 
distance  in  figures  I.  1  to  I.  4  for  a  total  of  750  radio  paths,    separating  the  frequency  ranges 
40-150  MHz,     150-600    MHz,    600-1000    MHz,    and  1-10    GHz.     Major  sources  of  data  other 
than  those  referenced  by  Herbstreit  and  Rice   [1959]   are  either  unpublished  or  are  given  by 
Bray,    Hopkins,    Kitchen,    and  Saxton  [1955] ,    Bullington  [19  55],    du  Castel  [19  57b] ,    Crysdale 
[1958],    Crysdale,    Day,    Cook,    Psutka,    and  Robillard  [19  57] ,    Dolukhanov  [1957],    Grosskopf 
[1956],    Hirai  [1961a,  b],    Josephson  and  Carlson  [19  58] ,    Jowett  [19 58  ] ,    Joy  [19 58a,  b] ,    Kitchen 
and  Richmond  [19  57],    Kitchen,    Richards,    and  Richmond  [19  58  ] ,    Millington  and  Isted  [19  50  ] , 
New^ton  and  Roger s   [1953],    Onoe,    Hirai,    and  Ni'wa  [1958],    Rowden,    Tagholm,    and  Stark 
[1958],    Saxton  [1951],    Ugai[l96l],    and  Vvedenskii  and  Sokolov  [19  57] . 

Three  straight  lines  were  determined  for  each  of  the  data  plots  shown  in  figures  I.  1 
to  1.4.     Near  the  transmitting  antenna,     A  =  0     on  the  average.      Data  for  intermediate  dis- 
tances,    where  the  average  rate  of  diffraction  attenuation  is  approximately    0.09  f^  db    per 
kilometer,    determine  a  second  straight  line.      Data  for  the  greater  distances,    where  the  level 
of  for-ward  scatter  fields  is  reached,    determine  the  level  of  a  straight  line  with  a  slope  varying 
from  1/18  to  1/14  db  per  kilometer,    depending  on  the  frequency. 

The  dashed  curves  of  figures  I.  l-I.  3   show  averages  of  broadcast  signals  recorded  at 
2500  random  locations  in  six  different  areas  of  the  United  States.      The  data  were  normalized 
to  10-meter  and  300-meter  antenna  heights,    and  to  frequencies  of  90,    230,    and  750   MHz. 
For  this  data  sample  [TASO  19  59],    average  fields  are  low  mainly  because  the  receiver  loca- 
tions were  not  carefully  selected,    as  they  were  for  most  other  paths  for  which  data  are  shown. 

The  extremely  large  variance  of  long-term  median  transmission  loss  values  recorded 
over  irregular  terrain  is  due  mainly  to  differences  in  terrain  profiles  and  effective  antenna 
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heights.      For  a  given  distance  and  given  antenna  heights  a  v/ide  range  of  angular  distances  is 
possible,    particularly  over  short  diffraction  and  extra-diffraction  paths.     Angular  distance, 
the  angle  between  radio  horizon  rays  from  each  antenna  in  the  great  circle  plane  containing 
the  antennas,    is  a  very  important  parameter  for  transmission  loss  calculations,    (see  section 
6).      Figure  I.  5  shows  for  a  number  of  paths  the  variability  of  angular  distance  relative  to 
its  value  over  a  smooth  spherical  earth  as  a  function  of  path  distance  and  antenna  heights. 

Most  of  the  "scatter"  of  the  experimental  long-term  medians  shown  in  figures  I.  1   - 
1.4  is  due  to  path-to-path  differences.     A  small  part  of  this  variation  is  due  to  the  lengths  of 
the  recording  periods.      For  all  data  plotted  in  the  figures  the  recording  period  exceeded  two 
weeks,    for   630  paths  it  exceeded  one  month,    and  for  90  paths  recordings  were  made  for  more 
than  a  year. 

An  evaluation  of  the  differences  between  predicted  and  measured  transnnission  loss 
values  is  discussed  briefly  in  annex  V.     In  evaluating  a  prediction  method  by  its  variance 
from  observed  data,    it  is  important  to  remember  that  this  variance  is  strongly  influenced  by 
the  particular  data  sample  available  for  comparison.      Thus  it  is  most  important  that  these 
data  samples  be  as  representative  as  possible  of  the  wide  range  of  propagation  path  conditions 
likely  to  be  encountered  in  the  various  types  of  service  and  in  various  parts  of  the  world. 

To  aid  in  deciding  whether  it  is  worthwhile  to  use  the  point-to-point  prediction  method 
outlined  in  sections  4  -   10,    instead  of  simpler  methods,    figure  I.  6  shows  the  cumulative  dis- 
tribution of  deviations  of  predicted  from  observed  long-term  median  values.      The  dash-dotted 
curve  shows  the  cumulative  distribution  of  deviations  from  the  lines  drawn  in  figures  I.  1   - 
1.4  for  all  available  data.      The  solid  and  dashed  curves  compare  predictions  based  on  these 
figures  with  ones  using  the  point-to-point  method  for  the  same  paths.      Note  that  the  detailed 
point-to-point  method  could  not  be  used  in  many  cases  because  of  the  lack  of  terrain  profiles. 

Figure  I.  6  shows  a  much  greater  variance  of  data  from  the  "empirical"  curves  of  fig- 
ures I.  1  -I.  4  for  the  sample  of  750  paths  than  for  the  smaller  sample  of  217  paths  for  which 
terrain  profiles  are  available.      The  -wide  scatter  of  data  illsutrated  in  figure  1.4  for  the  fre- 
quency range   1-10  GHz    appears  to  be  mainly  responsible  for  this.      Figure  1.4  appears  to 
show  that  propagation  is  much  more  sensitive  to  differences  in  terrain  profiles  at  these  higher 
frequencies,    as  might  be  expected.      The  point-to-point  prediction  methods,    depending  on  a 
number  of  parameters  besides  distance  and  frequency,    are  also  empirical,    since  they  are 
made  to  agree  with  available  data,    but  estimates  of  their  reliability  over  a  period  of  years 
have  not  varied  a  great  deal  with  the  size  of  the  sample  of  data  made  available  for  compari- 
son with  them. 
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I.  Z    standard  Point-to-point  Transmission  Loss  Curves 

A  set  of  standard  curves  of  basic  transmission  loss  versus  distance  is  presented  here. 
Such  curves  may  be  useful  for  establishing  preliminary  allocation  plans  but  they  are  clearly 
not  suitable  for  use  on  particular  point-to-point  paths,    since  they  make  no  allo-wance  for  the 
■wide  range  of  propagation  path  terrain  profiles  or  atmospheric  conditions  which  may  be  en- 
countered.    Similar  curves  developed  by  the  CCIR  [l963g;  1963h]  are  subject  to  the  same 
limitations. 

The  standard  curves  show  basic  transmission  loss  versus  path  distance  for  p  =  0.01 
to  99.99  percent  of  the  time.  These  curves  vv^ere  obtained  using  the  point-to-point  predic- 
tions for  a  smooth  earth,  N  =301,  antenna  heights  of  30  meters,  and  estimates  of  oxygen, 
water  vapor,  and  rain  absorption  explained  in  section  3.  Cumulative  distributions  of  hourly 
median  transmission  loss  for  terrestrial  links  may  be  read  from  figures  I.  7  to  I.  17  for  dis- 
tances from  0  to  1000  kilometers  and  for  0.  1,  0.2,  0.5,  1,  2,  5,  10,  22,  32.5,  60  and 
100    GHz.      The  same  information  may  be  obtained  from  figures  I.  18  to  I.  20. 

For  earth-space  links,    it  is  important  to  know  the  attenuation  relative  to  free  space, 

A,    between  the  earth  station  and  space  station  as  a  function  of  distance,    frequency,    and  the 

angle  of  elevation,      9,   ,    of  the  space  station  relative  to  the  horizontal  at  the  earth  station 

[CCIR  19631;  1963j].     Using  the  CCIR  basic  reference  atmosphere",    [CCIR  Report  231,     1963e] 

standard  propagation  curves  providing  this  information  for  2,    5,    10,    22,    32.  5,    60  and  100 

GHz,    for  values  of    p    ranging  from    0 .  0 1  to  99 .  99  per  cent,    and  for     6=0,    0.03,    0.1,    0.3, 

1.0,    and  TT  /  2    radians  are  shown  in  figures  I.  21  -  I.  26,    where    A     (p)     is  plotted  against  the 

straight-line  distance,      r,    between  antennas.      The  relationship  between    A         and    L       is 

m  b 

given  by 

L        (p)  +  A     (p)  +  L,  ^  =  A     (p)  +  32.  446  +  20  log  f  +  20  log  r      db  (I.  1) 

bm  m  bf  m 

■where    f    is  the  radio  frequency  in  megahertz  and     r     is  the  straight-line  distance 
between  antennas,    expressed  in  kilometers. 

The  curves  in  figures  I.  7  -  I.  26  provide  long-term  cumulative  distributions  of  hourly 
median  values.     Such  standard  propagation  curves  are  primarily  useful  only  for  general 
qualitative  analyses  and  clearly  do  not  take  account  of  particular  terrain  profiles  or  particular 
climatic  effects.      For  example,    the  transmission  loss  at  the  0.1%  and  0.01%  levels  will  be 
substantially  smaller  in  maritime  climates  ■where  ducting  conditions  are  more  common. 


The  transmission  loss  predictions  for  this  atnnosphere  are  essentially  the  same  as 


predictions  for     N     =301, 
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ANGULAR    DISTANCE  VERSUS  DISTANCE  FOR  THE  290  PATHS   FOR  WHICH 
TERRAIN  PROFILES  ARE  AVAILABLE 

THE  CURVES  SHOW  ANGULAR   DISTANCE,  9,  AS  A  FUNCTION  OF  DISTANCE 
OVER  A  SMOOTH  EARTH  OF  EFFECTIVE  RADIUS  =  9000  KILOMETERS 
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THE  WIDE  SCATTER  OF  THE  DATA  ON  THIS  FIGURE  ARISES  ALMOST  ENTIRELY 

FROM  DIFFERENCES  IN  TERRAIN  PROFILES,  AND  ILLUSTRATES  THE 

IMPORTANCE  OF  ANGULAR  DISTANCE  AS  A  PREDICTION  PARAMETER 


f 

/> 

V 

K- 

lOm,  hr  = 

IOm-~^^> 

// 

„ 

' 

// 

. 

' 

/^ 

A^h, 

=  l50m,  h 

r  =  IOm 

y 

■; 

a. 

. 

/ 

/  ' 

' 

- 

' 

• 

°-§    Y 

y-- 

% . 

% 

■0 

rg  9                     ° 
\ 

, 

/^ 

* 

200  300  400  500  600  700  800  900  1000 

PATH    DISTANCE,  d,  IN  KILOMETERS 


1200 


1300 


Figure  1.5 


1-8 


CUMULATIVE    DISTRIBUTION    OF   DEVIATIONS    OF   OBSERVED 
FROM    PREDICTED    VALUES   OF    TRANSMISSION    LOSS 
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0> 


mr 

::  £:: 

XI 

Tm 

rTTTTTTm 

rTTVTTn 

■■■ 

... 

..!.. 

-j- 

• 

... 

■■ 

.,... 

-- 

-■ 

1 

1 

1 

_L 

> 

ID 

i 

' 

ID 

Ul 

- 

3) 
o 

t 

""  ID    ■ 

■"     O 

1 

' 

H 

O 

1 

Tl 

/ 

T 

■  ■  --JI 

O 

71 

■Til"' 
3)    -- 

m  -■ 
O  .. 
c 

m  "" 
z    - 
o  -- 
-<  .. 

O)  "" 

o  -■ 

CD  "" 

I ;; 

3j 

""  "" 

O) 

TT 

J 

iLL 

m 

m 

:g 

m 

— 

o  CO 


c 


tyj^mj^l  Ml 

ilfilli 

:^:=-'y|H''::::::::::::::: 

1 

T 

T 

■ 

■  -  K  ^Jjm 

1               11 '" 

.■■!!■''■ 

; :' '  '         s  " !  •  •'' 

T 

,. -J. -"T. ,.-■:,  !- 

'•■h'-   -  ■  •  7 

""                          I 

i 

^0 
(X> 

s 

HMJy 

■'i'''  t''  ■■' ,' '.''  " 

/ 

itlit- 

:..::;::,i=::::^:;:;^:::;?  /■    ": 

^iS[_ 

""TTT  [rriTlTTT 

::::::;?:::;:::■;?!::•!■;'';::::  :::" 

■    ,'    .(    ■'  ,'  -' 

.     . 

.,:..,/.,(..-,.. ::::: 

/...>....,.'...,'...■..  ,•.:.. 



I 

>'    .     .'I'  t 

1 

...,..,..,..,.-, 

, 

^::;;.:::^;lv:::::::::::::::: 

/  ,!,/,, 

, 

" ;■"  /' 

1    /  '  1 

;  ., .. 

^:::;!.;i;^:::::::::::::::::::::::::::": 

,'  .    , 

/  /  ,' 

,..:,: 

;,.^^  :::::::: :::::::::::::::::::::::: 

1 

;.,.,.. 

., ., 

JJJ, 
-n 

3)  - 
m  - 

O    - 

c 
m 

2  - 
o  - 

K    . 

O   - 
o  . 

0  - 

1  - 

.     ,'  (    1    I 

,'  /  i  •-' 

, .,  .,.,.,. 

.!.... 

-D 

.,'  ;^  y  ' 

....    /  (   ,  '   ;  ( 

T 

n 

,.'...  ,l-'-i 

1 

2 

H 
O 

I 
n 

...,.,'./..   '.t.,! 

..._.  o_ 

..,(..<.-.,,,.,.. 

^.^.. 

,,/..(..,.., 

'    <   ,     ' 

T    L 

J" 

/             '   *    i'        — 1  TT 

.^..f-.,!.'  ^TtVi-S  .. 

n 

.<.,(,:,  ..S^.II.. 

1   ,'  '  '     ~ 

" 

'   II. 

!,!!::::::::::: 

t! J 

. 

- 

n 

D 

n 
n 

0 

□ 
n 

'"       -: 

c 

I 

Ll 

> 

o 

"D 

o 


^  ^  O 

>  en 

p  CO  O 

]>  o  c: 

C  r  < 

H  o  m 

^  CO  CO 

CO 


r-19 


BASIC  TRANSMISSION  LOSS  IN  DECIBELS 


31 


*-•        c 


CO 


... 

.. 



■- 

<''"J'i:'':'';''|:::^-'::'''" 

09 

,( t.i  ..1  ../..i/.'.'/.i'.'. 

a'  /■■''' 'I'r/^' 

'''Z----mm ■::::::: 

..■■•■■.,,::: 

'■'■■'''''■'■''''■'''' ■'''X'Y 

'                                     ** 

yf/jrv ;;;;;;;;;;;;;;;; 

• /' : 

: , 'L;  ■  i // 

t 

*    1              / 

/   1!;, 

t     1 
1 

/ 

..>.. 

1    III 

...  .,...t..A...i.^ 

'f              '                    

/ 

1 
i 

i\ 

1 
1 

1 

1  "1           '    , 
.  ./ ...J.. .[ . 

■1  1 

/ 

* 

'      /  i 

./    .1' 

1 

1 
1 

I 

1 

1 

\            1     1 

1 

[ ...    '  Jj.. 

1 

0 

/ 

V 

1 

' ..  ..  ...j..., .... 

1 
1 

Si 

^ 

•     ' 

1 
1 

I                |l     /        1 

1  ! 

/ 

frt 

f    1 

1      1 

/ 

ill  ii-  g  J 

rv,  .0     X 

-n 

1 



..    ..' 11 

ze ' 

m 

2 

1 

1  . 

1 

II 

■< 

1   ..    .  L 

/ 

1 

' 

1 
1 

i 

1     1    1 

I 

/     / 

1     ...  .. 

/ 

1/ 

i 

(                1 

/ 

1 

,                  ,     ..      1...      . 

1 

/ 

I 

[ ...j 1| 

■ /■ 

/ 

/ 

j 

1.1 ,.. 

1 

'               f      ^ 

: 

...1.  ..  ,,j 

..<..L. . 

.....  .1 

// 

J 

1. , 

'/ 

1 

/ 

,  1 

1 ,  r 



,...' [, 

g 
o 

H 

m 
Sen  3> 

-^p 

°^^ 
§lm 
5c«  w 


1-20 


LONG  TERM  VARIABILITY  ABOUT  THE   MEDIAN.  IN  DECIBELS 
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STANDARD  PROPAGATION  CURVES  FOR  EARTH-SPACE  LINKS 


NO  ALLOWANCE  HAS  BEEN  MADE  FOR  GROUND   REFLECTION 
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STANDARD   PROPAGATION   CURVES     FOR   EARTH    SPACE   LINKS 

00  =   0.03  RADIANS 

NO  ALLOWANCE  HAS  BEEN  MADE  FOR  GROUND  REFLECTION 
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STANDARD  PROPAGATION  CURVES  FOR  EARTH-SPACE  LINKS 

00  =  0,3  RADIANS 

NO  ALLOWANCE  HAS  BEEN  MADE  FOR  GROUND  REFLECTION 
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STANDARD     PROPAGATION      CURVES    FOR     EARTH -SPACE    LINKS 

6^0=    1.0  RADIAN 
NO    ALLOWANCE     HAS     BEEN     MADE    FOR    GROUND    REFLECTION 
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STANDARD   PROPAGATION    CURVES    FOR  EARTH-SPACE   LINKS 

NO  ALLOWANCE  HAS  BEEN   MADE   FOR   GROUND  REFLECTION 
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Annex    II 

BEAM  ORIENTATION,    POLARIZATION,   AND  MULTIPATH  COUPLING  LOSS 
Ordinarily,    to  minimize  the  transmission  loss  between  two  antennas,    they  are  oriented 
to  take  advantage  of  maximum  directive  gains  (directivity)  and  the  polarizations  are  matched. 
This  maximizes  the  path  antenna  gain.     With  a  single  uniform  plane  wave  incident  upon  a  re- 
ceiving antenna,    there  will  be  a  reduction  in  the  power  transferred  if  the  antenna  beam  is  not 
oriented  for  maximum  free  space  gain.     If  the  polarization  of  the  receiving  antenna  is  matched 
to  that  of  the  incident  wave,    this  loss  in  path  antenna  gain  is  due  to  "orientation  coupling  loss", 
and  if  there  is  a  polarization  mismatch,    there  will  be  an  additional  "polarization  coupling  loss". 
In  general,    more  than  one  plane  vi/ave  ■will  be  incident  upon  a  receiving  antenna  from  a  single 
source  because  of  reflection,    diffraction,    or  scattering  by  terrain  or  atmospheric  inhomoge- 
neities.     Mismatch  between  the  relative  phases  of  these  waves  and  the  relative  phases  of  the 
receiving  antenna  response  in  different  directions  will  contribute  to  a  "nnultipath  coupling  loss" 
•which  will  include  orientation,    polarization,    and  phase  mismatch  effects.     If  multipath  propa- 
gation involves  non-uniform  waves  whose  amplitudes,    polarizations,    and  phases  can  only  be 
described  statistically,    the  corresponding  loss  in  path  antenna  gain  will  include  "antenna-to- 
medium  coupling  loss",    a  statistical  average  of  phase    incoherence  effects. 

This  annex  indicates  how  multipath  coupling  loss  may  be  calculated  when  incident  waves 
are  plane  and  uniform  with  known  phases,    and  when  the  directivity,    polarization,    and  phase 
response  of  the  receiving  antenna  are  known  for  every  direction.     It  is  assumed  that  the  radia- 
tion resistance  of  the  receiving  antenna  is  unaffected  by  its  environment,   and  that  the  electric 
and  magnetic  field  vectors  of  every  incident  w^ave  are  perpendicular  to  each  other  and  per- 
pendicular to  the  direction  of  propagation. 

II.  1     Representation  of  Complex  Vector  Fields 

Studying  the  response  of  a  receiving  antenna  to  coherently  phased  plane  waves  with 
several  different  directions  of  arrival,    it  is  convenient  to  locate  the  receiving  antenna  at  the 
center  of  a  coordinate  system.     A  radio  ray  traveling  a  distance    r   from  a  transmitter  to  the 
receiver  may  be  refracted  or  reflected  so  that  its  initial  and  final  directions  are  different. 
If   -f   is  the  direction  of  propagation  at  the  receiver,      r  =  r  r     is  the  vector  distance  from  the 
receiver  to  the  transmitter  if  the  ray  path  is  a  straight  line,    but  not  otherwise. 

A  paper  by  Kales  [  1951]    shows  how  the  amplitude,    phase,    and  polarization  of  a  uniform, 

monochromatic,    elliptically  polarized  and  locally  plane  wave  may  be  expressed  with  the  aid 

of  complex  vectors.     For  instance,    such  a  wave  may  be  expressed  as  the  real  part  of  the  sum 

of  two  linearly  polarized  complex  plane  waves    \J~Z  e     exp(iT)     and    i\r2"  e .  exp  (i  t).      These 

r  1 

components  are  in  time  phase  quadrature  and  travel  in  the  same  direction   -f,    where    i  =  \l  -  1 

and   e      and   e.    are  real  vectors  perpendicular  to   r.      The  vector    e     +  ie.    is  then  a  complex 
ri  "-       '■  ,  ri 

vector.     Field  strengths  are  denoted  in  volts /km    (10     microvolts  per  meter)    and  field  inten- 

2  -3 

sities  in  watts /km       (10        milliw^atts  per  square  meter),     since  all  lengths  are  in  kilometers. 
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The  time-varying  phase 

T  =  k(ct  -  r)  (n.  1) 

is  a  function  of  the  free- space  wavelength  \,  the  propagation  constant  k  =  Ztt/X,  the  free- 
space  velocity  of  radio  waves  c  -  299792.5  ±0.3  km/sec,  the  time  t  at  the  radio  source, 
and  the  length  of  a  radio  ray  between  the  receiver  and  the  source. 

Figure  II- 1  illustrates  three  sets  of  coordinates  which  are  useful  in  studying  the  phase 
and  polarization  characteristics  associated  with  the  radiation  pattern  or  response  pattern  of 
an  antenna.     Let    r  =  f  r     represent  the  vector  distance  between  the  antenna  and  a  distant  point, 
specified  either  in  terms  of  a  right-handed  cartesian  unit  vector  coordinate  system    x   ,    x   , 
X      or  in  terms  of  polar  coordinates    r,    9,    <\>: 

—-  2  2  2  2 

r  =  f  r  =  X    X     +  X    X     +  X    x   ,  r  ■    =  x      +  x      +  x  (II.  2a) 

0OII22  0x2 

X     =  r  cos  9,  X     =  r  sin  9  cos  c|),  x     =  r  sin  9  sin(j)  {II.  2b) 

r=r(9,  <t))=x    cos  9  +  (x    coscb  +  x^  sincb)  sin9  (n.  2c) 

o  1  2 

Asa  general  rule,    either  of  two  antennas  separated  by  a  distance    r    is  in  the  far  field  or 
radiation  field  of  the  other  antenna  if    r  >  2D    /\,     where    D   is  the  largest  linear  dimension 
of  either  antenna. 

The  amplitude  and  polarization  of  electric  field  vectors    e      and   e  , .    perpendicular 
to  each  other  and  to    r,     is  often  calculated  or  measured  to  correspond  to  the  right-handed 
cartesian  unit  vector  coordinate  system    r,    e    ,    e      illustrated  in  figure  II- 1.     The  unit  vector 
e,    is  perpendicular  to    r    and   x   ,    and   e„   is  perpendicular  to    e,    and   r.     In  terms  of  vector 
cross-products : 

e  ,    =  (r  X  X    )/sin  9  =  x,  sin  6  -  x.  cos  6  (II.  3a) 

90  12 

e.   =  e     X  f  =  (x     -   f  cos  9)/sin9  (II.  3b) 

9         cj)  o 

The  directive       gain     g,     a  scalar,     may  be  expressed  as  the  sum  of    directive    gains      g        and 

0 

g        associated  with  polarization  components      e     see        and     e     =  e    e    ,     where  the  coeffi- 
(p  9  6     9  cp         (j)    9 

cients      e        and     e        are  expressed  in  volts /km : 
9  c{) 

g^gg    +     g^  (II.   4) 

Subscripts    t   and    r    are  used  to  refer  to  the  gains    g     and    g      of  transmitting  and  receiving 

^  ^  2 

antennas,   while    g    is  the  ratio  of  the  available  mean  power  flux  density  and    e     /r|    ,    where 


n-2 


e       as  defined  by    (II.5)    is  the  free  space  field  strength  at  a  distance     r     in  kilometers  from 

o 
an  isotropic  antenna  radiating     p       watts  : 

2     1/ 
e      =  [ti    p    /(4uT    )]'^     volts /km  (H.  5) 

-7 
Here,     ri     =    4itc  .  10        =  376.7304  ±  0.0004  ohms  is  the  characteristic  impedance  of  free 
o 

space.     The  maximum  amplitudes  of  the  9  and   c(>   components  of  a  radiated  or  incident  field 

are        e       ^/X     and        e       \f^ ,     where 
9  <p 

\^a\-^a-^    Sq      volts /km,  |T  |    =  e      =  e    g'f^     volts /km  (11.6) 

If  phases    t      and    t      are  associated  with  the  electric  field  components    e      and   e    ,    which 
9  (p  9  cp 

are  in  phase  quadrature  in  space  but  not  necessarily  in  time,    the  total  complex  wave  at  any 
point   r    is 

sfT  (  e     +  ie.)  exp(i  t)  =  \/T  [  e„  exp(i  t    )  +  e  ,  exp(i  t   )  ]   exp(i  t)  (IL  7) 

r  1  '■     9  9  (p  <p 

From  this  expression  and  a  knowledge  of    e_  t        ,    we  may  determine  the  real  and  imaginary 

_^  _^  9,  c|),      9,  4) 

components    e      and   e.,    which  are  in  phase  quadrature  in  time  but  not  necessarily  in  space: 

e     see     =  e„  cos  t.  +  e  ,  cos  t  ,  (11.  8a) 

r  r     r  9  9         (()  (j) 

e.  =  e.  e.  =  e„  sin  t„  +  e  ,  sin  t  (II,  8b) 

1110         9       4>         4> 

The  next  section  of  this  annex  introduces  components  of  this' wave  Which  are  in  phase  quadrature 
in  both  time  and  space. 
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n.  2    Principal  and  Cross-Polarization  Components 

Principal  and  cross-polarization  components  of  an  incident  complex  wave 

\fZ' (  e     +ie.)exp(iT)    may  be  defined  in  terms  of  a  time-independent  phase    t.    which  is  a 
r  i  1 

function  of   r    [Kales,    1951].     If  we  write 

~e     +  i~e.  =  (  e     +  ie    )  exp(i  t.)  (II.  9) 

r  1  1  2  1 

and  solve  for  the  real  and  imaginary  components  of  the  complex  vector    e     +  ie    ,     we  find  that 

e     see     =  e    cos  t.  +  e.  sin  t.  (II.  10a) 

1  1     1  r  11  1 

e      =  e,  e,  =  e.  cos  t.  -  e    sin  t.  (II«  10b) 

2  2     2         1  1         r  1 


Whichever  of  these  vectors  has  the  greater  magnitude  is  the  principal  polarization  component 

e    ,    and  the  other  is  the  orthogonal  cross-polarization  component   e    : 
p  "  ^  ^  c 

2  2  2  2         2  — ►      — ► 

e,    =  e     cos     T.  +  e.    sin    t.  +  e     •  e.  sin(2T.)  (II.  11a) 

Ir  11  iri  1 

2  2         2  2  2  — ►      — ► 

e^    =  e      sin    t.  +  e.    cos     t.  -  e     •  e.  sin(2T,)  (II.  lib) 

2r  11  iri  1 


The  phase  angle    t.    is  determined  from  the  condition  that    e     '6     =0: 

tan(2T.)  =  2^    .7./(e^-e.^)  (11.12) 

1  r        1        r  1 

Any  incident  plane  wave,    traveling  in  a  direction   -f   is  then  represented  as  the  real 
part  of  the  complex  wave  given  by 

\f2"Texpri(T  +  T.)l    =  \rT  (T    +  iT  )expri(T  +  T.)l  (n.  13) 

—  "•  i-"  p  c  1 

\ 
The  principal  and  cross-polarization  directions    e      and   e      are  chosen  so  that  their  vector 


product  is  a  unit  vector  in  the  direction  of  propagation: 


e     xe     =  -r  (H.  14) 

P         c 
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A  bar  is  used  under  the  symbol  for  the  complex  vector    e  s  e     +  ie     in  (11.13)    to  distinguish 

it  from  real  vectors  such  as    e,.,    e,,    e    ,    e.,    e    ,    and    e    .     The  absolute  values  of  the  vector 

9        <p        r        1        p  c 

coefficients    e      and    e      may  be  found  using  (11.11). 
p  c 

As  the  time   t  at  the  transmitter  or  the  time    t  at  the  receiver  increases,   the  real 
vector  component  of   (11.13),    or    "polarization  vector"  , 


\fZ'  fe     cos(T+T.)-e     sin(T+T.)l 
■■    p  1  c  1  -■ 

describes  an  ellipse  in  the  plane  of  the  orthogonal  unit  vectors    e     =  e    /e       and    e     =  e    /e    , 
^  ^  ^  pp'pccc 

Looking  in  the  direction  of  propagation    -r  (9,<j))    with  e      and   e      both  positive  or  both  negative, 
we  see  a  clockwise  rotation  of  the  polarization  vector  as    t   increases. 

Right-handed  polarization  is  defined  by  the  IRE  or  IEEE  and  in  CCIR  Report  321  [  1963m]    to 
correspond  to  a  clockw^ise  rotation  of  a  polarization  ellipse,   looking  in  the  direction  of  prop- 
agation with  r   fixed  and   t   or    t  increasing.     This  is  opposite  to  the  definition  used  in  classical 
physics. 

The    "axial  ratio"    e^/e       of  the  polarization  ellipse  of  an  incident  plane  wave 
\r2~_e  exp  [i(T  +  T   )]    is  denoted  here  as 

\^\^%  (n.  15) 

and  may  be  either  positive  or  negative  depending  on  whether  the  polarization  of  the  incident 

wave  is  right-handed  or  left-handed.     The  range  of  possible  values  for   a      is    -1  to  +1. 

x 
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a       = 

-e       i 

'e      , 

e        = 

-e       a 

xr 

cr 

pr 

cr 

pr    xr 

II.  3     Unit  Complex  Polarization  Vectors 

If  the  receiving  antenna  were  a  point  source  of  radio  waves,    it  would  produce  a  plane 

wave    \rz'  e    exp[i(T  +  t    )]     at  a  point    r    in  free  space.     The  receiving  pattern  of  such  an 

antenna  as  it  responds    to  an  incident  plane  wave    \fZ'  eexp[i(T  +  t    )]     traveling  in  the  opposite 

direction   -f   is  proportional  to  the  complex  conjugate  of   e    exp  (i  T   )     [S,   A.   Schelkunoff  and 

— r  r 

H,    T.   Friis.    1952]: 

[T   exp(iT   )]*  =  (T       -  iT     )exp(-iT   )  (n.  16) 

— r  r  pr  cr  r 

The  axial  ratio    e       /e  of  the  type  of  wave  that  would  be  radiated  by  a  receiving  antenna 

cr      pr 

is  defined  for  propagation  in  the  direction   r.     An  incident  plane  wave,    however,    is  propagating 
in  the  direction   -r,    and  by  definition  the  sense  of  polarization  of  an  antenna  used  for  reception 
is  opposite  to  the   sense  of  polarization  when  the  antenna  is  used  as  a  radiator.     The  polarization 
associated  with  a  receiving  pattern  is  right-handed  or  left-handed  depending  on  whether   a 
is  positive  or  negative,   -where 

{II.  17) 

The  amplitudes    |e       |    and    |e       |     of  the  principal  and  cross-polarization  field  components 

e         and    e         are  proportional  to  the  square  roots  of  principal  and  cross-polarization  directive 
pr  cr 

gains    g         and   g      ,     respectively.     It  is  convenient  to  define  a  unit  complex  polarization 
vector    p     which  contains  all  the  information  about  the  polarization  response  associated  with 
a  receiving  pattern: 

p     =  (e       +ie       a      )(l+a^    )"'^  (H,  18) 

— r  pr  cr    xr  xr 

a        =  g       /g  (n.  19) 

xr       "cr      pr  ^  ' 

The  directions    e         and   e         are   chosen  so  that 
pr  cr 

e       xe        =r  (H.  20) 

pr  cr 

In  a  similar  fashion,    the  axial  ratio   a      defined  by  (TT    1  S^  and  the  orientations     e     and  e 

X  '   ^     •       '  p  c 

of  the  principal  and  cross- polarization  axes  of  the  polarization  ellipse  completely  describe 
the  state  of  polarization  of  an  incident  wave    \r2~£exp[i(T  +  t.)],     and  its  direction  of  propa- 
gation   -f  =  e     X  e    .     The  unit  complex  polarization  vector  for  the  incident  wave  is 
p         c 
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^=T/|T|=(e     +  i?    a    )(1  +  a^)"'/^  (n.21) 

—      —    '  —  '  p  ex  X 

The  magnitude  of  a  complex  vector    e  =  e     +  ie       is  the  square  root  of  the  product  of   e    and 
its  complex   conjugate     e     -  ie    : 

ITI    =  (T.  T*/'   =  (e  ^  +  e  ^'/^        volts /km  (n.  22) 

—  —    —  p  c 
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n.4    Power  Flux  Densities 

The  coefficients    e      and   e      of  the  unit  vectors    e      and   e      are  chosen  to  be  r.m,  s. 
p  c  p  c 

values  of  field  strength,    expressed  in  volts/km,    and  the  mean  power  flux  densities    s      and 

s      associated  with  these  components  are 
c 

?  7  2  Z 

s      =  e     /ri  watts/km    ,  s      =  e     /ri  watts/km  (11.23) 

p         p      'o  ceo 

The  corresponding  principal  and  cross-polarization  directive  gains     g^    and    g^     are 

g     =  4iTr    Sp/p^,  g^  =  4iTr    sjp^  (U.  24) 

where     p       is  the  total  power  radiated  from  the  transmitting  antenna.     This  is  the  same  rela- 
tion as  that  expressed  by    (11.6)    between  the  gains     g   ,     g    ,     and  the  orthogonal  polarization 

9         cp 

components     e       and     e    . 

9  9 

The  total  mean  power  flux  density    s    at  any  point  where    e    is  known  to  be  in  the  radi- 
ation field  of  the  transmitting  antenna  and  any  reradiating  sources   is 

I  — ►.  2  2  2  2 

s=e       /ri     =ge     /ri     =s     +s      =(e      +e     )/r| 
'  —  '         'o  o      'o         p         c  p  c        'o 

=  (e  ^  +  e,^)/r|     =  (e^^  +  e^/r]  watts/km^  (H.  25a) 

r  1  o  0  (|)         o 

a  -  s     +g     =g„  +  g,=  4T:r    s  /p^    =  s  t]    /e  (!!•  ^^b) 

^       ^p      ^c       ^9       ^(f)  t  'o'    o 

where     e       is  given  by    (II.5).     The  power  flux  density     s     is  proportional  to  the  transmitting 
o 

antenna  gain     g  ,     but  in  general     g     is  not  equal  to     g      as  there  may  be  a  fraction     a       of 
t  t  p 

energy  absorbed  along  a  ray  path  or  scattered  out  of  the  path.     We  therefore  write 

2.  ,,    .       2,  (n.  26) 


,<^  +  ^x^=%St^'^"x)  =  ^pgt 


The  path  absorption  factor    a       can  also  be  useful  in  approximating  propagation  mechanisms 

ir 

which  are  more  readily  described  as  a  sum  of  modes  than  by  using  geometric  optics.     For 
instance,   in  the  case  of  tropospheric  ducting  a  single  dominant  TEM  mode  may  correspond 
theoretically  to  an  infinite  number  of  ray  paths,   and  yet  be  satisfactorily  approximated  by 
a  single  great-circle  ray  path  if     a       is  appropriately  defined.     In  such  a  case,      a       will 
occasionally  be  greater  than  unity  rather  than  less. 
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Orienting  a  receiving  dipole  for  maximum  reception  to  determine    s      and  for  minimum 

reception  to  determine    s      will  also  determine    e      and    e    ,     except  in  the  case  of  circular  polar- 
^  c  p  c 

ization,    where  the  direction  of   e      in  the  plane  normal  to    r    is  arbitrary.     In  the  general 

case  where     |a    ]    <  1,     either  of  two  opposite  directions  along  the  line  of  principal  polarization 

is  equally  suitable  for    e    . 

p 

Reception  with  a  dipole  will  not  show  the  sense  of  polarization.     Right-handed  and  left- 
handed  circularly  polarized  receiving  antennas  will  in  theory  furnish  this  information,    since 
e    may  also  be  written  to  correspond  to  the  difference  of  right-handed  and  left-handed  circu- 
larly polarized  waves  v/hich  are  in  phase  quadrature  in  time  and  space: 

I-  (% '  ^%)  (^4^  -  ^(\  -  ^  V  fv^j  ^"-''^ 

The  mean  power  flux  densities    s      and    s      associated  with  right-handed  and  left-handed  polar- 
izations are 

s     =  (e     +e    )^/(2ti   )    watts/km^  (n.  28a) 

r  p         c  o 

s.   =  (e     -  e    )^/(2ti    )    watts/km^  (H.  28b) 

i  p         c  o 

so  the  sense  of  polarization  may  be  determined  by  whether    s    /s      is  greater  than  or  less  than 
unity.     The  flux  densities    s      and    s     are  equal  only  for  linear  polarization,    where     e      =0, 
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II.  5    Polarization  Efficiency 
The  polarization  efficiency  for  a  transfer  of  energy  from  a  single  plane  wave  to  the 
terminals  of  a  receiving  antenna  at  a  given  radio  frequency  may  be  expressed  as  a  function  of 
the  unit  complex  polarization  vectors  defined  by    (11.18)    and    (11.21)    and  the  angle    ^i       between 

ir 

principal  polarization  directions  associated  with     e_    and     e_  .     This  polarization  efficiency  is 


A     A 


2                                 2              2  2 

cos    di    {a    a       +  1)     +  sin    ^Jj    (a     +  a      ) 
2  _     p      X    xr p      X         xr 


£-£r'      -  ,     2   .    ,.,     2 


(n.29) 


(a      +  l)(a       +  1) 
X  xr 


where 


e     •  e        =  -e     •  e        =  cos  Jj   ,  e     •  e        =  e       •  e      =  sindi  (11.30) 

p       pr  c        cr  p  p        cr  pr        c  p 

As  noted  in  section  2  following  (2.27),    any  receiving  antenna  is  completely    "blind" 
to  an  incoming  plane  wave    n/^  e  exp[i(T  +  t.)]     which  has  a  sense  of  polarization  opposite 
to  that  of  the  receiving  antenna  if  the  eccentricities  of  the  polarization  ellipses  are  the  same 
Ma    I    =  ja      1 1      and  if  the  principal  polarization  direction   e      of  the  incident  wave  is  perpen- 
dicular to   e      .In  such  a  case,     cos  Uj     =  0,     a     =  -a      ,     and  (11.29)   shows  that  the  polarization 
pr  p  X  xr 

efficiency     |p'  P    |       is  zero.     As  an  interesting  special  case,    reflection  of  a  circularly 

polarized  wave  incident  normally  on  a  perfectly  conducting  sheet  will  change  the  sense  of 

polarization  so  that  the  antenna  which  radiates  such  a  wave  cannot  receive  the  reflected  wave. 

In  sue  ha  case    a     =-a        =±1,    so  that     ^  •  p  =0    for  any  value  of   iJj    . 

x  xr  '—    — r'  p 

On  the  other  hand,    the  polarization  efficiency  given  by  (11.29)  is  unity  and  a  maximum 

transfer  of  power  will  occur  if    a      -a         and    Jj     =0,     that  is,    if  the  sense,    eccentricity, 

X         xr  p 

and  principal  polarization  direction  of  the  receiving  antenna  match  the  sense,    eccentricity, 
and  principal  polarization  direction  of  the  incident  wave. 

For  transmission  in  free  space,    antenna  radiation  efficiencies,    their  directive  gains, 
and  the  polarization  coupling  efficiency  are  independent  quantities,    and  all  five  must  be 
maximized  for  a  maximum  transfer  of  power  between  the  antennas.     A  reduction  in  either  one 

_,     ir)    or  a  reduction  in  the  polarization  etticiencv     IP*  D    I 
will  reduce  the  transfer  of  power  between  two  antennas. 

With  each  plane  wave  incident  on  the  receiving  antenna  there  is  associated  a  ray  of 
length   r    from  the  transmitter,    an  initial  direction  of  radiation,    and  the  radiated  wave 
e    exp[i(T  +  T  )]     which  would  be  found  in  free  space  at  this  distance  and  in  this  direction. 
When  it  is  practical  to  separate  antenna  characteristics  from  environmental  and  path  charac- 
teristics,   it  is  assumed  that  the  antenna  phase  response    t  ,    like    T   ,    is  a  characteristic 

t  r 

of  the  antenna  and  its  environment  and  that 


of  the  directive  gains     g(-r)    and    g    (r)    or  a  reduction  in  the  polarization  efficiency     jp-  p    | 
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T.  =  T  +  T  {n.3i) 

1        t        p 

where    t      is  a  function  of  the  ray  path  and  includes  allowances  for  path  length  differences  and 
diffraction  or  reflection  phase  shifts. 

Random  phase  changes  in  either  antenna,    absorption  and  reradiation  by  the  environ- 
ment,   or  random  fluctuations  of  refractive  index  in  the  atmosphere  will  all  tend  to  fill  in  any 
sharp  nulls  in  a  theoretical  free- space  radiation  pattern   e    or    e    .     Also,    it  is  not  possible 
to  have  a  complex  vector  pattern   e/r   which  is  independent  of   r    in  the  vicinity  of  antenna  nulls 

unless  the  radiation  field,    proportional  to    1/r,    dominates  over  the  induction  field,    which 

2 
is  approximately  proportional  to    1/r    . 
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n.  6    Multipath  Coupling  Loss 

Coherently  phased  multipath  components  from  a  single  source  may  arrive  at  a  receiving 

antenna  from  directions  sufficiently  different  so  that    t.    and    t      vary  significantly.     It  is 

then  important  to  be  able  to  add  complex  signal  voltages  at  the  antenna  terminals.     Let    n  =  1, 

2, ,    N  and  assume   N   discrete  plane  waves  incident  on  an  antenna  from  a  single  source. 

The  following  expressions  represent  the  complex  open-circuit  r.m.s.    signal  voltage    v 

n 

corresponding  to  a  radio  frequency    v   cycles  per  second,    a  single  incident  plane  wave 

n/"?  e    exp[i(T  +  t.    )] ,     a  loss-free  receiving  antenna  with  a  directivity  gain    g 

and  an  effective  absorbing  area  a      ,   matched  antenna  and  load  impedances,   and  an  input 

en 

resistance     r      which  is  the  same  for  the  antenna  and  its  load: 


V     =  (4r     s    a      )'^  (p    •  p      )  exp  ri(T+T       +t      -t      )1     volts  (11.32) 

n  V    n    en         — n    — rn  *•  pn         tn         rn  ■■ 

s     =   I  e    I     /ri     =  p    a       g^     /(4Trr     )    watts/km  (11.33) 

n       '  — n'         'o         r    pn    tn  n 

a        =  g       k^/{4n)    km^  (n.34) 

en      °rn  *  ' 

p    .p       =  [(1  +  a^    )(1  +  a^      )l''^[(l+a       a        )  cos  +       +i(a       +a        )sin4j      1     (n.35) 
— n    — rn       "•  xn  xrn  ■■  ■•  xn    xrn  pn  xn         xrn  P^i 


If  the  polarization  of  the  receiving  antenna  is  matched  to  that  of  the  incident  plane  wave,    then 

a       =a        ,+       =0.     P'P       =1>     3-nd 
xn         xrn         pn  — n        rn 

v=r4rpa       g^g       \    /(4TTr    )    V^expri(T+T       +t       -t      )]     volts  (n.  3^) 

n       ■•        V    r     pn    tnrn  n     ^  *■  pn         tn        rn  -i 

If  the  coefficient  of  the  phasor  in    (11.36)    has  the  same  value  for  two  incident  plane  waves,   but 
the  values  of     t.     -  r         differ  by     tt     radians,   the  sum  of  the  corresponding  complex  voltages 
is  zero.     This  shows  that  the  multipath  coupling  efficiency  can  theoretically  be  zero  even 
when  the  beam  orientation  and  polarization  coupling  are  maximized.     Adjacent  lobes  in  a 
receiving  antenna  directivity  pattern,   for  instance,   may  be      180      out  of  phase  and  thus  cancel 
two  discrete  in-phase  plane -wave  components. 

Equation  (2.6)  in  section  2  shows  the  relation  between  the  total  open-circuit   r.m.s. 
voltage 


volts  (n.  37) 


N 

% 

V        = 

V 

'    Y    V    V*  " 

^       n    m 

n=  1 

m  =  1 
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and  the  power     p       available  at  the  terminals  of  a  loss -free  receiving  antenna: 

3, 

p     =  V    /(4r     )    watts  (JI.  38) 

a         V  V 

In  writing     p       for     p  in    (11.37),   the  subscript     v     has  been  suppressed,   as  with  almost 

a  av 

all  of  the  symbols  in  this  annex.     Studying    (11.32)  -  (11.35),     (11.37),     and    (11.38),     it  is  seen 

Z 
that  the  expression  for     p       is  symmetrical  in  the  antenna  gains     g    ,     g      ,     and    g     =  a     g    , 

g       =  a       g      ,      and  that     p       is  a  linear  function  of  these  parameters,   though     v       is  not. 
°cr         xr  °pr  ^a  '^  '  "  v 

From  this  follows  a  theorem  of  reciprocity,   that  the  transmission  loss     L  =  -  10  log  (p    /p  ) 

a      L 

is  the  same  if  the  roles  of  the  transmitting  and  receiving  antennas  are  reversed. 

The  basic  transmission  loss     K      is  the  system  loss  that  would  be  expected  if  the 
actual  antennas  were  replaced  at  the  same  locations  by  hypothetical  antennas  which  are: 

(a)  loss-free,    so  that     L  ^  =  L       =0  db.       See    (2.3), 
^    '  et  er  *        ' 

(b)  isotropic,    so  that     g.  =  g     =1     in  every  direction  important  to  propagation  between 
the  actual  antennas. 

(c)  free  of  polarization  coupling  loss,    so  that      |p  •   P    |     =1     for  every  locally  plane 
wave  incident  at  the  receiving  antenna, 

(d)  isotropic  in  their  phase  response,    so  that     t    =  t     =0     in  every  direction. 

The  available  power     p  ,      corresponding  to  propagation  between  hypothetical  isotropic 
antennas  is  then 


x2     N      ,  ,%  , 

p    X       ^^      (a       a        )'     cos(t       -  t 


Pab 


t Y  pn    pm pn         pm  ^^^  ^ 


(4t7)       ^-',  n    m 

n  =  1 

m  =  1 


The  basic  transmission  loss    L      corresponding  to  these  assumptions  is 

b 

L^=-10  1og(p^^/p^)  =  P^-P^^      db  (n.40) 

The  basic  transmission  loss  in  free  space,    L,  .,    corresponds  to   N  =  1,     a        =  1,      t        =0, 

^  bf  ^  pi  pi 

"l   ="  = 

L       =   -  10  log      [\/(4TTr)]^  =  32.45  +  20  log  f  +  20  log    r      db  (n.  41) 

where   f   is  in  megacycles  per  second  and   r    is  in  kilometers.      Compare  with  (2.31). 

As  may  be  seen  from  the  above  relations,    only  a  fraction    s      of  the  total  flux  density 

s      per  unit  radiated  power    p      contributes  to  the  available  received  power    p     from   N  plane 

2  2 

waves.     While    s      is  expressed  in  watts/km    ,     s      is  expressed  in  watts/km     for  each  watt 
n  e 
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of  the  power    p      radiated  by  a  single  source: 

For  each  plane  wave  from  a  given  source,     e     expfir.    )    or    p       exp(-iT      )    may  sometimes 
^  ^  — n       ^         m'  -rn       ^  rn'  ^ 

be  regarded  as  a  statistical  variable  chosen  at  random  from  a  uniform  distribution,    with  all 

phases  from    -it   to   it    equally  likely.      Then  real  power  proportional  to     I  e     •  e  *   |        may  be 

— n     — rn 

added  at  the  antenna  terminals,    rather  than  the  complex  voltages  defined  by  (11.32)  -(11.35). 

For  this  case,    the  statistical    "expected  value"     <s    >     of     s        is 

e  e 

N 

<s>=Va       gg        |p-p       I     /(4TTr     )  (n.43) 

e  ^_^        pn     tn     rn    — n     — rn'  n 

n=  1 

In  terms  of    s    ,    the  transmission  loss    L    is 
e 

L  =  21.46  +  20  log  f  -   10  log     s-        db  (H.  44) 

Substituting     <s    >     for      s        in    (11.44),    we  would  not  in  general  obtain  the  statistical  expected 

e  e 

value      <L>     of     L,    since      <L>     is  an  ensemble  average  of  logarithms ,    which  may  be  quite 

different  from  the  logarithm  of  the  corresponding  ensemble  average      <Sg>.      For  this  reason, 

median  values  are  often  a  more  practical  measure  of  central  tendency  than     "expected" 

values.      With     p        and    \     fixed,    median  values  of     s        and     L,     always  obey  the  relation 

t  e 

(11.44),    while  average  values  of     s        and     L     often  do  not. 

The  remainder  of  this  appendix  is  concerned  with  a  few  artificial  problems  designed 
to  shov/  how  these  formulas  are  used  and  to  demonstrate  some  of  the  properties  of  radiation 
and  response  patterns.     In  general,    information  is  needed  about  antenna  patterns  only  in  the 
few  directions  ■which  are  important  in  determining  the  amplitude  and  fading  of  a  tropospheric 
signal.     Although  section  II.  7  shows  how  a  complex  vector  radiation  or  reception  pattern  may 
be  derived  from  an  integral  over  all  directions,    it  is  proposed  that     the  power  radiation 
efficiencies  and  the  gains    g    (f)    or    g    (-r)   for  actual  antennas  should  be  determined  by  meas- 
urements in  a  few  critical  directions  using  standard  methods  and  a  minimum  of  calculations. 
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II.  7    Idealized  Theoretical  Antenna  Patterns 
Consider  a  point  source  of  plane  waves,    represented  by  complex  dipole  monaents 
in  three  mutually  perpendicular  directions,     x   ,    x   ,    and  x    .      These  three  unit  vectors, 
illustrated  in  figure  II.  1,    define  a  right-handed  system,    and  it  is  assumed  that  the  corre- 
sponding elementary  dipoles  support   r.m.s.     currents  of   I   ,    I    ,    and  I^    amperes,    respectively. 
The  corresponding  peak  scalar  current  dipole  moments  are    \r2  I     i    ampere-kilometers, 

where    m  =  0,    1,    2,    and  the  sum  of  the  complex  vector  dipole  moments    x      \rz  1     Jlexpiir     ) 

m  m         m 

may  be  expressed  as  follows: 

a"  =  a"    +  i  a"  _  (H.  45a,) 

r  =  »/tm(x  c^  +  X  c  +  5  cj,     r  =  Nr2n(x  s^  +  X  s  +  5c  s.)  cl 45b) 


I^  =  I„^  +  I,^  +  I,  ,  c       =  {I     /I)  COST     ,  s       =  (I     /I)sinT     ,  m  r:  0,   1,  2.       (11.46) 

012  mm  m  mm  m 

Here,     t_^,    t   ,     and    t      represent  initial  phases  of  the  currents  supported  by  the  elementary 
dipoles.     The  time  phase  factor  is  assumed  to  be     exp(ikct). 

Using  the  same  unit  vector  coordinate  system  to  represent  the  vector  distance    r    from 
this  idealized  point  source  to  a  distant  point: 

r=xx^  +  xx     +x^x     =rr  (11. 47) 

o     (J  11  2     2 

where     x„,   x   ,     and     x       are  given  by    (n.2b)    as  functions  of     r,    9,    cj).     The  complex  wave  at 
T     due  to  any  one  of  the  elementary  dipoles  is  polarized  in  a  direction 

r  X  (x       X  r)   =  X       -   r  X     /r  (11.  48) 

m  mm 

which  is  perpendicular  to  the  propagation  direction   f   and  in  the  plane  of   x       and   r.     The 
total  complex  wave  at    r   may  be  represented  in  the  form  given  by  (11.8)  : 

\rz    e  (r)  exp(i  t)  :=  n/^  (e     +  i  e.)  exp  (i  t)  =  n/"!"  (e     +  i  e    )  exp  [  i  (t  +  t  )] 

=  [r  X  (Tx  r)]  [Ti^/(2>.r)]  exp(iT)  (n.  49) 

T  =  k(ct  -  r)  +  tt/4  (n.  50) 

\fTT    =  [a",    -  r("r    •  r)l  n    /(2Xr)    volts/km  (H.  51a) 

r  1  i  o 
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^fZ    e.   =  [a     -   f  (a     .  r)]  r]    /(Z\r)    volts/km 


(n. 51b) 


The  total  mean  power  flux  density    s(r)    at    r    is  given  by  (11.13) 


s  (T)  =  (e^^  +  e.^/Ti^  =  [a^^  -  U^  o  f)^  +  a.^  -  (T^ .   f)^]  /^^ 


'±K    r 


{n.52) 


c         =  (I     I   /I   )  cos{t      -  t   ) 

mn       m  n  m        n 


(II.  53) 


The  total  radiated  power     p       is  obtained  by  integrating     s(r)     over  the  surface  of  a  sphere 
of  radius     r,    using  the  spherical  coordinates     r,    0,    cj)     illustrated  in  figure  II,  1  : 


iTT  „  TT 


f  "  2      — 

\        de  r     s(r)  sine 


2-Tr  Ti  (li) 


3\ 


watts 


{U.  54) 


From    (11.54)    it  is  seen  that  the  peak  scalar  dipole  moment     nTSIj?     used  to  define      a^^     and 
"a        in    (11.45)    may  be  expressed  in  terms  of  the  total  radiated  power  : 


^r^li  -     \n/  3p  TIttti   )     ampere -kilometers 


(n.  55) 


The  directive  gain     g(r)     is 


2      —  ,  3 

g(r)  =  4-n-  r     s(r)/p^   =    j 


1  -     (  —r-   I      COS 


"M  2  „  2 

sin     9  COS    i 


Z\  2  2         f  \  2 

—  )     sin    0  sin   tj)-(ccoS(j)  +  c       sincj)  )  sin(2  0)  -  c       sin    0    sin(2(})) 


(IL  56) 


This  is  the  most  general  expression  possible  for  the  directive  gain  of  any  combination  of 
elementary  electric  dipoles  centered  at  a  point.  Studying  (II. 56),  it  may  be  shown  that  no 
combination  of  values  for     I   ,   I   ,    I   ,    t   ,     t   ,     t-,     will  provide  anisotropic  radiator.     As 

(J  1  ^  U  i.  ^ 

defined  in  this  appendix,   an  isotropic  antenna  radiates  or  receives  waves  of  any  phase  and 
polarization  equally  in  every  direction. 
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For  the  special  case  where     I^  =  I,   =  I,  =  I/nTT,      t     =  tt/Z,      t,   -  0,      and     t,  -  ir, 
(11.57)    shows  that 

g(r)  =   1  +  sin    9    sin  ej)    cos  ^  (IX,  57) 

With  these  specifications,  (H..46)  shows  that     c     =  0,      c     =  -  c     =  l/sTT,      s     =  l/sTJ, 
s^  =  s^  =  0,      and  (11.45)    with    (11.55)    shows  that 

a     =   (x    -  X  )b  ,        a-,  =   X   b  (11.  58a) 


b  =  \ 


P^A-^^o)   I  (IL58b) 


Substituting  next  in    (11.51)    with  the  aid  of    (II.2): 


\r2~e     =  e    (x,  -  X.  -  r  b,)   ,         n/T  e.  =  e    (x    -  r  cos  6)  (IL  5i9) 

rolzZ  100  *' 


e     = 
o 


2    ./^ 


2    1 ' 
r|     p    /  (4'ir  r    )  ,         ^?  ~  ^"^  ®  (cos  cj)  -  sin  <j) ) 


(11.60) 


The  principal  and  cross -polarization  gains  deternained  using    (11.24)    and    (11.2  5)    are 

2  2 

g   (r)  =    1  +  sin    e(sincj)  cos  i^  -  Vz)  ,         g   (r)  =  ^/^  sin     9  (11.  61) 

The  subscripts     p     and     c     in    (IL61)    should  be  reversed  whenever     g(9,  <j))     is  less  than 

2 

sin   9.     Minimum  and  maximum  values  of     g     are      1/2     and     3/2      while     g       ranges  from 

ir 

1/3     to     1     and     g       from     0     to     1/2. 

The  importance  of  phases  to  multipath  coupling  is  more  readily  demonstrated  using 
a  somewhat  more  complicated  antenna.     The  following  paragraphs  derive  an  expression  for 
a  wave  which  is  approximately  plane  at  a  distance     r     exceeding     200     wavelengths,    radiated 
by  am.  antenna  composed  of  two  three-dimensional  complex  dipoles  located  at     -5    X.  x       and 
+5  X  X       and  thus  spaced     10     wavelengths  apart.     When  the  radiation  pattern  has  been 
determined,   it  will  be  assumed  that  this  is  the  receiving  antenna.     Its  response  to  known 
plane  waves  from  two  given  directions  will  then  be  calcixlated. 

With  the  radiated  power     p       divided  equally  between  two  three-dimensional  complex 
dipoles,     a    is  redefined  as 

T=(b/\rT)a'     ,        "a     =x,  -x_+ix  (H.  62) 

—  — O  — ff.  ICO 
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Since     5\    is  negligible  compared  to     r     except  in  phase  factors  critically  depending  on 
r     -  ^-,>     the  exact  expressions 


r,=r-5Xx      ,         r-=r+5Xx 
1  o  2  o 


(n,  63) 


lead  to  the  following  approximations  and  definitions: 


r^  =  r(l-0.      r2  =  r(l+6),      e   =  5{\/r)  cos  9 


^  ^  . 


r,   =r(l+e)  -x  e   sec  9,       r=r(l-6)+xe   sec  9 
1  o  ^  o 


^       ^ 


r  =  x^  cos  9  +  (x     cos  c|)  +  x     sine)))  sin  9 


(n.  64) 
(11.65) 
(n.  66) 


For  distances     r     exceeding     200  wavelengths,       [e  |    <     0.025     and     e       is  neglected  entirely, 
so  that 


r,r,=r-5Xx      ,         r^r,  =  r  +  5X.x 
11  o  2     2  o 


(n.67) 


At  a  point     r,   the  complex  wave  radiated  by  this  antenna  is  approximately  plane  and 
may  be  represented  as 


sir  f  e     +  i  e.  j    exp{i  t)  =  \[T       e    exp(i  t  )  +  £2  ®^P(i  "^2^ 


(n.  68) 


where 


T  =  k(ct-  r)  +  Tr/4 


(n.  69) 


T,     =T+T       ,  T^    =    T-T       ,  T       =10  IT    COS   9 

1  a  2  a  a 


(n.  70) 


As  in    (11.49),   the  waves  radiated  by  the  two  main  elements  of  this  antenna  are  represented  in 
(11.68)    as  the  product  of  phasors     exp(iT.)     and     exp(iT-)     multiplied  by  the  complex  vectors 
sTZ  e        and     n/~2  e     ,      respectively: 


nTS"  e^  =      ?^X  (  a  X  r^ 


T,   /(2Xr)  =  (e    /2) 
o  o 


a     -  r     I  a     •    r 
—  o         1  V—  o        1 


(n. 71a) 


^±2 


r^X 


ax?^)!  Ti^/(2Xr)  =[ejZ) 


a     -  r     1  a     •    r 
—  o  2  V-o  2 


(n.  71b) 
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Evaluating    T   •?,,     7   .  r„     7^,      and    7^     with  the  aid  of    (n.62),   (n.65),    (n.66)    and    (11.71): 


a"    •    r,  =  b-,(l+€)  +  i    cos  9  -  €  (sec  9  -  cos  9) 
—  o  1         ^ 


a     .    r^  =b   (l-€)  +i 
—  o         2         ^ 


cos  9  +  £(sec  9  -  cos  9) 


(n.72a) 
(n. 72b) 


7^  =  (e^/2)  /Tx^  -S^  -?b2(l  +  26)  +^^^2  ^^^^^ 


+  i 


X  ( 1  +€  )  -  r  cos  9  +  r  e  (sec  9-2  cos  9)       I 


(Jl.  73a) 


e     =  (e    /2) 

—  ii  o 


X,   -  x_  -  rb-  (1  -  2t)  -  X    b-  €  sec  9 
1         ^  Z  o    2 


+  i 


X  ( 1  -  €  )  -  r  cos  9  -  r  €  (sec  9-2  cos  9) 
o 


(n.  73b) 


Since  the  sum  and  difference  of     exp(iT,)     and     exp(iT-)     are       2  cos  t       exp(iT)     and 

12a 

Zisinr     exp(iT),     respectively,      e       and     e.     as  defined  by    (11.68)     are 


'o{ 


A  ^  ^  _ 

^1-^2  "^^2 


cos    T      -  € 

a 


x     +  r (sec 
o 


2  cos  9)       sin  t 


(n. 74a) 


e.    =    e 
1  o 


X     -  r  cos  9      cos  t     -     b- 

o  1  a  2 


r  -  X     sec  9 
o 


(U.  74b) 


The  complex  wave      '/2(  e     +  ie.   1    is  a  plane  wave  only  when     e       and     e.     are  both  per- 
pendicular to  the  direction  of  propagation,      r,      or  when 


r  •   I  e     +  le.    I    =  6   sin  t 
r  1/  a 


(cos  9  -  sec  9)  +  i  sin  9  (cos  cj)  -  sin  (j>) 


=  0 


(n.  75) 


which  requires  that     e  =  0,      sin  t     =  0,      or      9=0.      If    e      is  negligible,   the  total  mean 
power  fliix  density  in  terms  of  the  directive  gain     g(r)     is  given  by 


-*■        /     2  z\  A        2 

5(  r  )  =  (^e^    +  e.    J/r^^  =  g(r)  e^  / x]^ 


(11.76) 
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2  2 

g(r)  =2(1+  sin     9  sin  (j>  cos  cj))  cos     t  (n.  77) 


That     p       is  the  corresponding  total  radiated  power  may  be  verified  by  substituting    (11.76) 
and    (no77)    in    (n.54),   with  the  aid  of    (11.60)    and    (n.70). 

Now  let  this  antenna  be  a  receiving  antenna,   and  suppose  that  direct  and  ground- 
reflected  waves  arrive  from  directions     r(9,  cj))     equal  to 


rf\         ^ 


r^(0.32,  Tr/4)  =  0.9492  x^  +  0.2224(x^  +  x^)  (n.78a) 

r   (0.28,    0.75)  =  0.9611  X    +  0.1430  x    +  0.1332  x  (n.78b) 

"  O  i.  A 


Note  that     r       and     r       in    (11.78)    are  not  related  to     r       and     r       in    (11.65),     but  are  two 

particular  values  of     r.     Corresponding  values  of     t   ,      cos  t   ,      and     sin  t       are 

a  a  a 

T       =  29.82111   ,         cos  T       =  -  0.0240   ,         sin  t       =  -  0.9997  (H.  79a) 

ai  ai  ai 

T       =  30.19245   ,         cos  T       =   0.3404      ,         sin  t       =  -  0.9403  (H.  79b) 

az  a.2  az 

The  incoming  waves  in  the  two  directions     r       and     r       are  assiuned  to  be  plane,   and  the 

distances     r       and     r^     to  their  source  are  assumed  large  enough  so  that       e  ,  sin  t 

12  i  ai 

and     €      sin  t         are  negligible  compared  to     cos  t         and     cos  t      ,   respectively.     The  plane 
2  az  ai  az 

wave  response  of  the  receiving  antenna  in  these  directions  may  be  expressed  in  terms  of  the 
complex  vectors  associated  with     r       and     r    : 


e      +ie.     =  -  0.024  e 
rl  XI  o 


(x,   -  X,)  +  i(0.099  X     -  0.211  X,   -  0.211  x   ) 
^12  o  1  i- 


(II.  80a) 


e"     +iT      =  0,340  e    T  (-0.013  S     +  0.998  x,   -  1.002  x.)  +  i(0.076  x     -  0.137  x^   -  0. 128  x) 
rziz  o  o  i  c.  o  i  ^ 


(n. 80b) 


Since    T      •  T      =  0  ,  (11. 11)  with  (II.9)  and  (11. 10)  shows  that     t      =0,     so  that     e  and     e. 

ri  rz  ^1  ^1  ^1 

are  principal  and  cross -polarization  components  of  the  complex  vector  receiving  pattern: 


T         +  i T         =T       +  iT.  (II.  81a) 

pri  cri  ri  il 
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These  same  equations  show  that     t      =  -  0.005     and  that 


rz 


e        +  i  e         =  0.340  e 
pr2  crz  o 


(-0,013  X     +  0.999  S^  -  1.001  x^)  +  1(0.076  x^  -  0.132  x^  -  0.113  x^)! 


(n. 81b) 


differing  only  slightly  from  (Il.SOb),  since  t  ^      is  almost  zero. 

The  axial  ratios  of  the  two  polarization  ellipses,   defined  by    (n.l7),    are 


a  =-0.222,  a  =-0.143  (U.  82) 

xri  xrz 


and  the  unit  complex  polarization  vectors    J         and    J  defined  by    (IL18.)    are  therefore 

p       =  0.674(x,   -  xj  +i(0.067x     -  0.142  x     -  0.142  x)  (H.  83a) 

p      =  (-0.009  Xq  +  0.692  x^  -  0.694  x^)  +  i(0.053  x^  -  0.095  x^  -  0.089  x^)  (n.83b) 


The  antenna  gains     g    (r    )     and    g    (r    )     are  given  by  (11.77): 
r      1  r    2 

g(r)=  0,0021,       g(r)=  0.241  (11.84) 

r     1  r     J 

which  shows  that  the  gain    G    (r   )  =   10  log  g    (r   )  associated  with  the  direct  ray  is     29.2  db 

r      1  r     1 

below  that  of  an  isotropic  antenna,    while  the  gain    G   (r   )     associated  with  the  ground-reflected 

r     2 
ray  is     -6.  2  db  .     It  might  be  expected  that  only  the  incident  wave  propagating  in  the  direction 

-r       would  need  to  be  considered  in  determining  the  complex  voltage 'at  the  receiving  antenna 

2 

terminals.     Suppose,    however,    that  the    ground-reflected   ray  has  been  attenuated  considerably 

more  than  the  direct  ray,    so  that  the  path  attenuation  factor     a  is    0.01,    while    a        =1. 

P2  Pi 

Suppose  further  that  the  transmitting  antenna  gain  associated  with  the  ground-reflected  ray  is 

6  db  less  than  that  associated  with  the  direct  ray.      Then  the  mean  incident  flux  density    s^ 

associated  with  the  ground-reflected  ray  will  be  26  db  less  than  the  fLux  density    s       associated 

with  the  direct  ray. 

In  order  to  calculate  the  complex  received  voltage    v      given  by  (11,37),    then,      the 

following  is  assumed: 


2  2 

r^  =  52  ohms  ,  ^\~^  watt/km    (  =  -  30  dbm/m   ) 
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2 
s     =  0.0025  watts /km     ,         \  =  0.0003  km(  f  =  1000  MHz  ) 

a       =  0.2  ,         a      =  0.4  ,         4^      =    7r/2  ,         \\i      =1.5 
XI  X2  pi  pz 

T.      =    T  +    T         =    0    ,  T.       =    T  +    T         =     TT  (II.  85) 

11  Pl  tl  12  pZ  t2  *  ' 

It  will  be  seen  that  these  assumptions  imply  a  more  nearly  complete  polarization  coupling  loss 
between  the  direct  ^vave  and  the  receiving  antenna  than  betw^een  the  ground-reflected  w^ave  and 
the  receiving  antenna.      The  effective  absorbing  area  of  the  receiving  antenna  for  each  wave, 
as  given  by  (11.34)  is 

a^^   =  1.504  X  10"       km     ,         a       =  1.726  X  lO"^  km^  (11.86) 


ei 


The  polarization  factors  are 

(Pl  '  iri  ^  "  ■  °-°^^  ^  •         (Pz  °  irz  ^  "  °-°^^  +  °-^^^  ^  <^-  ^"^^ 

and  the  phase  factors  are     exp  (i  t)     and     exp  [i(T+  3.137)  ],     respectively.     Substituting  these 
values  in    (n.32),   the  complex  voltages  are 

Vj^  =  -  1.175(10"    )i  exp(iT)  ,         v^  =  -(1.887  +  7.071i)(10"^)  exp(iT)  (n.  88) 

The  real  voltage  at  the  antenna  terminals,  as  given  lay  (11.37)     is 


\  =  (V]^  +  v^)  (v^  +  v^)*  =  8.33  X  10"^  volts  =  8.  33  microvolts  (n.  89) 


and  the  corresponding  power     p^     available  at  the  terminals  of  the  loss -free  receiving 
antenna  is 

'  -12 

^a  "  °*^^^  ^  ^^         ^^**^  '         P     =  -  125  dbw  =  -  95  dbm  (n.  90) 

as  given  by    (11.38). 
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n.  8     Conclusions 
The  foregoing  exercise  demonstrates  that: 

(1)  Even  small  changes  in  antenna  beam  orientation,   transmission  loss,   polarization 
coupling,  and  multipath  phasing  may  have  a  visible  effect  on  the  available  power  at  the 
terminals  of  a  receiving  antenna. 

(2)  If  the  formulation  of  the  general  relationships  for  a  completely  polarized  wave 
is      programmed  for  a  digital  computer,    it  may  be  feasible  to  estimate  the  complete  stat- 
istics of  a  received  signal  whenever  reasonable  assumptions  can  be  made  about  the  stat- 
istics of  the  parameters  described  in  this  annex. 

(3)  The  measurement  of  antenna  characteristics  in  a  few  critical  directions  will  often 
be  sufficient  to  provide  valuable  information  to  be  used  with  the  relationships  given  here. 
The  measurement  of  Stokes'  parameters,   for  instaxLce,   will  provide  information  about     a      , 

g   ,    ip   ,     and  both  the  polarized  field  intensity     s       and  the  unpolarized  field  intensity     s    . 
These  parameters    [Stokes,  1922   ]     are 

s     +s     =  total  mean  field  intensity  (11.91) 

r         o 

Q  =  s    cos(2B)  cos(2  i)j   >  (n.  92) 

r  p 

U=s    cos(2(3)  sin(2i|j  )  (U.  93) 

V  =  s    sin(2S)  (n.  94) 

r 

where 

P  =  tan"^a  (11,95) 

xr 

The  unpolarized  or  randomly  polarized  field  intensity     s       is  determined  from  (11.91)    and  the 


identity 


/a 


2  2  2/2 

s     =  (Q     +  U     +  V    )  (11.96) 


r 

Using  standard  sources  and  antenna  model  ranges,   the  gain     g       may  be  determined  from 

r 


't  =  ^J^^o^\^  '         ^o    "  \^J^'^''''^^  <J"-9^) 


assuming,    if     e       is  measured,   that  any  power  reception  efficiency     1/i         less  than  unity 

o  2  er 

will  affect     s       and     e         alike, 
r  o 
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Finally,    a  method  for  measuring  relative  phase  responses     t       is  also  needed.     In 

r 

individual  cases,    multipath  coupling  loss  may  be  instifficient  to  provide  adequate  \mwanted 

signal  rejection.     Variations  of     t       may  lead  to  phase  interference  fading  of  wanted  signals, 

r 

just  as  variations  of    a  are  associated  with  long-term  power  fading.     Because  of  the 

ir 

complexity  of  these  phenomena,   they  are  usually  described  in  terms  of  cximulative  distri- 
butions of  signal  amplitudes  or  fade  durations.     Fortionately,   even  crude  measurements  or 
simple  theories  may  then  suffice  to  provide  statistical  infornciation  about     t  . 


11-24 


COORDINATE     SYSTEMS     FOR 
STUDYING     ANTENNA     PATTERNS 


A 


Figure  II.  I 
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Annex  III 
FORMULAS,    COMPUTER  METHODS,    AND  SAMPLE  CALCULATIONS 
The  material  of  this  annex  is  organized  into  the  following  sections: 
1.     Line -of- sight 
Z.     Diffraction  over  a  single  isolated  obstacle 

3.  Diffraction  over  a  single  isolated  obstacle  with  ground  reflections 

4.  Diffraction  over  irregular  terrain 

5.  Forward  scatter 

6.  Forward  scatter  with  antennas  elevated 

7.  Long-term  variability 

8.  Sample  calculations 

Section  1  lists  geometric  optics  formulas  for  computing  transmission  loss  over  a 
smooth  earth,    for  determining  the  magnitude  and  phase  of  the  reflection  coefficient,    and  for 
computing  a  first  Fresnel  zone  along  a  great  circle  path.     Graphs  of  the  magnitude    R    and 
phase    c    of  the  reflection  coefficient  are  included.     Section  2  gives  mathematical  expres- 
sions that  approximate  the  curves    A(v,  0),     A(0,  p)     and    U(vp)     for  convenience  in  using 
a  digital  computer.     Section  3  lists  geometric  optics  formulas  used  to  compute  diffraction 
attenuation  when  several  components  of  the  received  field  are  affected  by  reflection  from  the 
earth's  surface.     Section  4  defines  the  parameters    K    and    b    for  both  horizontally  and 
vertically  polarized  radio  waves.     Section  5  shows  the  function    F(0d)     for    N     =  250,    301, 
3  50,    and  400,    and  for  values  of    s    from  0.01  to  1.     Curve  fits  to  the  function    F(0d)  and 
equations  for  computing    H   (r|     =0)     are  included.     Section  6  suggests  modifications  of  the 
prediction  methods  for  use  when  antenna  beams  are  elevated  or  directed  out  of  the  great 
circle  plane.     Section  7  shows  diurnal  and  seasonal  changes  in  long-term  variability. 
Mathematical  expressions  used  to  compute  predicted  distributions  are  shown  and  a  method 
of  mixing  distributions  is  described.     Section  8  illustrates  the  application  of  prediction 
methods  to  specific  problems: 

Example   1.     Line-of-sight  communication. 

Example  2.     Diffraction  over  an  isolated  rounded  obstacle. 

Example  3,     A  transhorizon  path:     diffraction  and  scatter  losses  are  computed  to 
obtain  a  reference  median  transmission  loss.     The  predicted  distribution  of  hourly  median 
values  of  transmission  loss  is  compared  with  observed  values. 
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m.  1     Line -of- Sight 

Simple  formulas  for  line-of- sight  propagation  which  suffice  for  most  applications,    are 
given  in  section  5  of  the  report.     Formulas  for  geometry  over  a  smooth  earth  and  for  deter- 
mining the  magnitude  and  phase  of  the  reflection  coefficient  are  given  here.      These  formulas 
may  be  used  when  the  great  circle  path  terrain  visible  to  both  antennas   will  support  a  substan- 
tial amount  of  reflection,   and  it  is  reasonable  to  fit  a  smooth  convex  curve  of  radius   a  to  this 
portion  of  the  terrain. 

Figure  5.  lb  illustrates  the  geometry  appropriate  for  reflection  of  a  single  ray  by  a 
smooth  earth  of  effective  radius    a.     In  the  figure,    t|j   is  the  grazing  angle  at  the  geometrical 
reflection  point  located  at  a  distance    d      from  an  antenna  of  height   h     and  at  a  distance   d 
from  an  antenna  of  height   h   .     The  total  path  distance    d  =  d     +  d      is  measured  along  an  arc 
of  radius   a.     The  difference,  Ar,  between  the  reflected  ray  path  length    r     +  r       and  the  length 
of  the  direct  ray,    r    ,    is  calculated  to  find  the  phase  of  a  radio  field  which  is  the  sum  of  ground- 
reflected  and  free  space  fields.     If  Ar    is  less  than  0,06k,    these  ray  optics  formulas  are  not 
applicable.     For  almost  all  cases  of  interest  the  angle    ^jj   is  small  and  the  straight  line  distances 
r    ,    r     and  r      are  very  nearly  equal  to  the  mean  sea  level  arc  distances    d    ,    d     and  d.     The 
geometric  optics  formulas  given  below  usually  require  double-precision  arithmetic. 

.1  \  ^1 

tan  4i  =  cot(d,  /a)  -  (1  +  h,  /a)       esc  (d,  /a)   ^  ^—    -    .r— 

1  1  Id  2a 

h^  d 

-1  2  2 

tan  4j  =  cot(d    /a)  -  (1  +  h   /a)       csc(d    /a)   ^  -y-    -    -r— 

2 
r^=  a|(h^/a)^+(h2/a)^-  2(h^/a)(h2/a)  +  2[1  +  h^/a  +  h^/a  +  (h^/aXh^/a)]  [  1  -  cos{d/a)]  1 


^1 


"2 


2 

(a  sin  qj)     +  h   (2a  +  h   ) 


-  a  sin  ij" 


(III. 

1) 

(m. 

2) 

)" 

(III. 

3) 

(ni.4) 

2  1  y^ 

(a  sin  4j)     +  h   (2a  +  h    )    \        -  a  sin  i|j  (m.  5) 


2 

Ar  =  r+r-r     =4rrsin     4'/(r+r+r)  (III.  6) 

J.^O  i.        ^  xuO 

Equating  (III.  1)  and  (III.  2)   and  substituting    d-d       for     d       in  (III.  2),    the  distance    d       may 
be  determined  graphically  or  by  trial  and  error,   and    tan  i|j    is  then  calculated  using  (HI.  1). 
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Using  double  precision  arithmetic,      (III.l)     through    (IU.6)    give  an  accurate  estimate  of  the 
path  difference     Ar     for  reflection  of  a  single  ray  from  a  smooth  earth.     This  value  is  then 
used  in    (5.1)    of  section    5    to  obtain  the  long-term  median  reference  value  of  transmission 

loss     L 

cr 

If  either     h^     or     h^     greatly  exceeds  one  kilometer,   and  if  it  is  considered  worth- 
while to  trace  rays  through  the  atmosphere  in  order  to  determine      i|j     more  accurately, 
values  of     d^      or     d^,     tabulated  by  Bean  and  Thayer    [  1959  ],    may  be  used.     Given     h   ,     h   , 

and  the  surface  refr activity,      N   ,     select  trial  values  for    tjj,     calculate     d       and     d         and 

s  12 

continue  until     d^  +  d^  =  d.     Then    (in.l)    and    (JILZ)     must  be  solved  for  new  values  of     h 

and     h^     if     (IU.3),     (in.4),     and    {ni.5)     are  used  to  obtain  the  path  difference, 

Ar  =  r,    +r^-r     . 
1  Z  o 

The  symbols     R     in    (5.1)    and     c     in    (5.4)    represent  the  magnitude  and  the  phase 
angle  relative  to     ir,     respectively,    of  the  theoretical  coefficient     R  exp[-i  (ir-c)  ]     for  reflec- 
tion of  a  plane  wave  from  a  smooth  plane  surface  of  a  given  conductivity      o-     and  relative  di- 
electric constant    e  .     Values  of    R     and    c     as  a  function  of  the  grazing  angle    ijj    are  shown  in 
figures    III.l    to    III. 8    for  vertical  and  horizontal  polarization  over  good,   average,   and  poor 
ground,   and  over  sea  water.     The  magnitude     R     of  the  smooth  plane  earth  reflection  co- 
efficient is  designated     R        or       R       for  vertical  or    horizontal  polarization  respectively, 
and  is  read  on  the  left-hand  ordinate  scale  using  the  solid  curves.     The  phase  angle  relative 
to     TT,     is  designated     c        or       c,      for  vertical    or    horizontal  polarization  respectively,   and 
is  read  in  radians  on  the  right-hand  scale  using  the  dashed  curves.     As  seen  from  these 
figures  in  most  cases  when  the  angle     \\i     is  small,      R     is  very  nearly  unity  and     c     may  be 
set  equal  to  zero.     A  notable  exception  occurs  in  the  case  of  propagation  over  sea  water 
using  vertical  polarization. 

In  preparing  figures    III.l    to    III. 8,     the  following  general  expressions  for  the  mag- 
nitudes    R       and     R,      and  lags  (ir  -  c    )     and     (it  -  c,  )     were  used.     In  these  equations, 
V                   h                  °  V  h 

€     is  the  ratio  of  the  surface  dielectric  constant  to  that  of  air,     or     is  the  surface  conductivity 
in  mhos  per  meter,      f     is  the  radio  frequency  in  megacycles  per  second,   ajid     4'     is  the 
grazing  angle  in  radians, 

X  =  1.  80  X   10^  0- /f,     q  =  x/(2p)  (IIL  7) 


z,'- 


% 


2       2         21  2 

(e    -  cos     4j)     +  X  +  (e    -  cos     ^)  (UI.  8) 


2         2. 

€      +  X 


V 

"       2         2    • 

p    +  q 

2(p€   +  qx) 

V 

2           2 
p      +  q 

\--T^  (in.  9) 

p      +  q 


m^  =       /P    ,  {in.  10) 

h  2.2 


p    +  q 


m-3 


R 


1+b     sin     41-m     sin 


,] 


2  --1 

1+b     sin     ijj  +  m     sini|i 


J  L      ^ 


(in.  11) 


R, 


1+b     sin     ijj  -  ni,    sin  4^       ] 
h  h 


1+b,    sin     di  +  m,    sin  d' 
h  h 


(in.  12) 


tan 


-If     X  sin  t];  -  1 
e   sin  ijj  -  P 


tan 


-1 


X  sin  4'  +  q 
€   sin  tjj  +  p 


(in.  13) 


tan 


sin  iji 


tan 


^ 


siniiJ  +  P    y 


(HL  14) 


The  angle    c^    is  always  positive  and  less  than  it,    and    c       is  always  negative  with  an  ab- 
solute magnitude  less  than    it.      The  pseudo  Brewster  angle,    where    c       suddenly  changes  from 


near  zero  to    tt  /  2,    and  where    R       is  a  minimum,    is     sin       \fTJh~ . 

V  V 

For  grazing  angles  less  than    0.1    radian,   for  overland  propagation,  and  for  fre- 
quencies above     30  Mc/s,    excellent  approximations  to    (III.ll)    and    (III.12)  are  provided  by 
the  f ollov/ing  formulas : 

R     =  exp(-m    x\i)  (in.  15) 


R^  =  exp  ( -  mj_^  i(;) 


(in.  16) 


The  assumption  of  a  discrete  reflection  point  with  equal  angles  of  incidence  and 
reflection  as  shown  in  figure    5,1    is  an  oversimplification.     Actiially,    reflection  occurs  from 
all  points  of  the  surface.     For  irregiolar  terrain,   this  is  taken  into  account  by  a  terrain 
roughness  factor      o"   ,     (subsection  5.  1 ),    which  is  the  r.  m.  s.    deviation  of  terrain  relative  to 
a  smooth  curve  computed  within  the  limits  of  a  first  Fresnel  zone  in  a  horizontal  plane.     The 
outline  of  such  a  Fresnel   ellipse  is  determined  by  the  condition  that  the  length  of  a  ray  path, 
r       +  r      ,     corresponding  to  scattering  from  a  point  on  the  edge  of  the  ellipse  is  half  a  wave 
length  longer  than  the  geometrical  ray  path,      r    +  r    ,     where  the  angles  of  incidence  and 
reflection  are  equal. 

The  first  Fresnel  ellipse  cuts  the  great  circle  plane  at  two  points,      x       and     x 

3.  D 

kilometers  from  the  transmitter.     Th«  distemces     x       and     x^     are  defined  by  the  relation    ' 


2         2  2  2         2  2 

r      sin    v|j+x    +or      sin    <\i  +  [  (r    +  r    )  cos  ^  -  x]      =r+r+   \/2  (HL  17) 


in-4 


The  exact  solution  for     x     is 


2x^^^(1  +  5)=  [(r^  +  r^Kl  +  S)  +  (r^-r2)]   cos  ^  ±  (t^  +  t^+ \/Z)    5\ll+4:T_T^/[{r^  +  r  ^)     5] 


where 


2         ,2  2  2         ,2  2  ,  ,2  ,,22 

r^    =h'      +d^     .         r^    =h'      +d2     ,  (r^+r^)     -    (h^ +h^)     +d 


cos  ^  =  d^/r^  =  d^/r^    ,         sin    ^Jj  =  h'^  /r^     =  h^  /i 


(ni.  18) 


41^+^2 


^1+^2 


/sin    ij; 


'^l'     *^2     3-^®  "defined  by    (5,7),     and     \     is  the  radio  wavelength  in  kilometers. 

As  an  alternative  method,   the  points     x       and     x^     may  be  computed  in  terms  of  path 
distance,   the  heights     h       and     h' ,      and  the  radio  frequency.     In  this  method,   the  distance 
X       to  the  center  of  the  first  Fresnel  zone  is  first  computed,    then  the  distance     x      from  the 
center  to  the  margin  of  the  zone  is  subtracted  from    x       to  give     x   ,     and  added  to  give     x_  . 

(in.  19) 


where 


Xq  =  d/2     1  +  B  f{h'^  -  h^  ) 


km 


0.3d(l+2h'^h^/d^)  +f(h|^+h^)^ 


-1 


(in.  20) 


X     =  0.548  B  d^l  Tf  h'h^/d  +  0.075  (1 +  2  hV  /d^) 


l+(h|+h^)^/d^ 


% 


-l  +  2h|h2/d' 


(ni.  21) 


X     =  x„  -  X,       km   ,         X,    =  x„  +  X,       km 
a         0         1  b         0         1 


III- 5 


The  method  given  in    (III.19)    to    (III. 21)    is  applicable  whenever     d   >>  \.     If  in  addi- 

I      1  2 

tion,      h    h     <  <    d   ,     the  computation  of     B,      and     x       may  be  simplified  as  follows: 

L       O  1 


B  = 


0.3d  +f(h'^  +  h^)^'' 


-1 


(in.  22) 


X     =  0.548  B  d 


f h'h^/d  +  0.075 


1  +(h'^+h2)^/d^ 


(m.  23) 
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in.  2    Diffraction  over  a  Single  Isolated  Obstacle. 

The  theoretical  diffraction  loss  curves  on  figures    7.1    to    7.4    have  been  fitted  by- 
arbitrary  mathematical  expressions  for  convenience  in  using  a  digital  computer. 

The  diffraction  loss  for  an  isolated  rounded  obstacle  and  irregular  terrain  is  given  in 
section  7  as: 

A(v,  p)  =  A(v,  0)  +  A(0,  p)  +  U(vp)      db  (7.7) 

v/here  the  parameter     v     is  defined  as 


V  =  ±  2V  Ar/X     S   ±  V  2d  a   6    /X  (7.1a) 

o   o 


V  =  ±  2.583  e«/Td^d^7d"  (7.1b) 


and     p      an  index  of  curvature  of  the  rounded  obstacles  is  defined  as: 

p    =  0.676  r-^^^  f"^^^  [d/r^r^)]'"^  (7.8) 


For  an  ideal  knife  edge,    (p  =  0),    the  diffraction  loss  is     A(vj  0)     and  is  shown  on 
figure  7.1.     For  values  of    v     from  minus  one  to  large  positive  values,   this  curve  may  be 
approximated  using  the  following  mathematical  expressions: 

For      - 1   <   V    <    0  , 
A(v,  0)  =  2,78v     +  9.95  v+ 6.17        db  .  (in.  24a) 


For        0    <    V    <    3  , 
A(v,  0)  =  6.08  -  0.005  v"*  +  0.159  v^  -  1.7  v^  +  9.3  V        db.  (in.  24b) 


For        V   >   3> 
A(v,  0)  =  12.953  +  20  log  v      db   .  (in.  24c) 

The  theoretical  curve  for     A(0,p).   is  approximated  by  : 

A(0,  p)  =  5.86  +  1.49  p^  +  1.46  p^  +  6.49  p         db   ,  (UI.  25) 

and  the  curve     U(vp)     is  approximated  as  follows: 


ni-15 


For    vp     <   2,  U(vp)  =  (43.6  +  23.5  vp)  log  (1  +  vp)  -  6  -  6,7  vp       db  .  (lUo  26a) 


For    vp    >  2,  U(vp)  =  22  vp  -  20  log  (vp)  -  14.13        db.  (IE.  26b) 

An  average  allowance  for  terrain  foregroiind  effects  may  be  made  by  adding  a  term 

10  exp(-2.3p)     to     A(0,  p).     This  term  gives  a  correction  which  ranges  from    10  db    for 

p    =  0     to      1  db    for     p    =   1. 

When  reflections  from  terrain  on  either  or  both  sides  of  the  obstacle  should  be  con- 
sidered,  the  method  given  in  the  following  section  may  be  used.     This  method  considers  the 
diffraction  loss  and  phase  lag  over  the  diffracting  obstacle,   and  the  path  length  differences 
and  reflection  coefficients  of  the  reflected  waves. 
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III.  3     Diffraction  over  a  Single  Isolated  Obstacle  with  Ground  Reflections 
Diffraction  over  an  isolated  obstacle  is  discussed  in  section  7,   where  ways  of  approx- 
imating the  effects  of  reflection  and  diffraction  from  foregroiind  terrain  are  indicated.      Where 
the  effects  of  reflection  are  expected  to  be  of  great  importance,    such  as  in  the  case  of  prop- 
agation over  a  large  body  of  water,    the  following  geometric  optics  method  miay  be  used. 

Figure    111.9    illustrates  four  distinct  ray  paths  over  a  knife  edge;  the  first  ray  is  not 
reflected  from  the  ground,    the  second  and  third  are  each  reflected  once,    and  the  fourth  ray 
is  reflected  once  on  each  side  of  the  knife  edge.     Each  ray  is  subject  to  a  diffraction  loss     f . 
and  a  phase  lag        ^.     at  the  knife  edge,   where     j  =   1,    Z,    3,    4.       Both     f.     and     $.     depend 
on  the  knife-edge  parameter     v     given  in  section    111.2.     When  the  isolated  obstacle  is  roxxnded, 
rather  than  an  ideal  knife  edge,    the  diffraction  loss  depends  on     v     and     p  ,   where     p      is  the 
index  of  curvature  of  the  crest  radius,    defined  in  section    111.2.      The  parameter     v     may  be 
written: 

V.  =  ±  2n/  A./\    S     ±  '/Td"^7F~.TX~  (in.  27) 

J  J  oj  ^oj  ^  ' 

where     A.     is 
J 

^1  =  ^1+^2-^0'     ^2  =  ^11+^12+^2-^02 
^3  =  ^1+  ^21  +  "22  -  ^03'     ^4=^1+  "12  +  ^21  +  ^22  "  "o4  ^^"-  '^^ 

Path  differences      A.     used  to  calculate     v.     in    (111.27)    are  closely  approximated  by 


the  following  formulas  : 


A.  =  d     0.^,     d     =dd^/(2d),      G.  =  6  +  e. 
J  r     J  r  1     2  J  jr 


^r  =  °'     ^2r  =  ^^l^l/'^l'      ^3r=  ^^22^2/^2'      ^r  =   ^2r  +  ^r  ^"^'^'^ 

The  total  phase  change     $  (v,  p)     at  an  isolated  rounded  obstacle  is 

2 
$.  =  ®.(v,  p)  =  90  V     +   c()(v,  0)  +  (j>(0,  p)  +  cj)(vp)  (in.  30a) 

where  the  fxinctions      cj)(v,  0),     cj)(0,  p),      and     cJ3(vp)     are  plotted  as  dashed  curves  on  figures 

7.1,     7.3,     and    7.4.     For  an  ideal  knife  edge,   where  the  radius  of  curvature  of  the  crest  is 

zero,      p    =  0,      and    (III.  30a)     reduces  to 

2 
for    v>0  <&.(v,  0)  =  90v     +  c|5(v,  0)  (in.  30b) 


for    v<0  $.(v,  0)   =  4)(v,  0)  (in.  30c) 
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The  three  components  of  the  received  field  which  are  affected  by  reflection  from  the 

earth's  surface  depend  also  upon  effective  ground  reflection  coefficients    R       exp  [-i(T7  -  c    )] 

and    R       exp  [-i(ir  -  c    )  ] ,   defined  in  section  5,   and  upon  ray  path  differences     A         and    A       : 
eo  J  2r  3r 

^2r  =  ^11  -^^12-^1^    ^^I'^l^z/^l 


^3r  =  "21  +  "22  -  "2  =    '  +2   ^21  ^^22/^2 


(in.  31) 


Usually,   it  may  be  assumed  that    c-  =  c,  =  0    so  that  the  reflection  coefficients  are    -R    _ 

2  3  e2 

and    -R    ,  . 
e3 

Introducing  the  propagation  constant    k  =  2  ir/X,   the  attenuation  relative  to  free  space 


is  then 


A  = 


-20  log  I    f^exp(-i$^)  -  R^2  f^  exp  ^i($^  +'^^2^1 


\3^3^"P 


-i($3+kA3)|    +R^2R^3f4expp($^  +  kA2+kA3) 


db      (in.  32) 


where 


C    -  S 

f.  =  n^l(l-C.-S.)^+(C.-S.)^     tan..=^-X-L. 


C.  -     C         cosT  ^  j    dt.     S.  -     r         sin(^-J    dt  (in.33) 

Pearcey  [  1956  ],   and  the  NBS  AMS  55  Handbook  of  Mathematical  Furtctions  [  1964]  give 

complete  tables,    series  expansions,   and  asymptotic  expressions  for  the  Fresnel  integrals 

C       and     S  .     The  magnitude     f.sf(v.)   for   v  =  v.   may  also  be  determined  from  figure  7.1 

j  j  J  J  3 

and  the  expression 

(m.  34) 


log  f(v.)  =  -  A(v.)/20 


and  where    v    is  larger  than  3; 


f.  =  0.22508/v.  ,     $.  =    —     1  +  2  v.         radians 
J  J  J        4  V  J     ' 


(ILL.  35) 


Figure    III.  10    is  a  nomogram  which  may  be  used  in  the  determination  of     f(v.)     and     $(v.) 
for  both  positive  and  negative  values  of     v.     This  nomogram  is  based  on  the  representation 
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of  Fresnel  integrals  by  the  Cornu  spiral. 

The  general  problem  requires  calculating      9,     d,     d   ,     d    ,     d,,>     '^■i->>     '^yi'     ^77' 

\\i   ,      and     iji   ,      as  shown  in  figure    111,9. 
1  ^ 

1.  Calculate     9.     and     A.     for     j  =  1,     2,     3,     4,      using    (111.29). 

J  ] 

2.  Calculate     v.,     C,     S.,     f.,      and     $.,     using    (IILZI),     (111.33),     and  figure    7.1. 

J  J  J  J  J 

3.  Calculate     A         and     A  from    (m.31), 

zr  3r 

4.  Calculate     R  and     R  from    (5.1),     or  assume  that     R       =  R       =1. 

ez  es  ez  es 

5.  Substitute  these  values  in    (111.32). 

To  check  the  calculation  of  each    v.,   the  approximation  given  in  (111.27)    may  be  used,   with 

the  following  formulas  for     a    .  =  d^  9./d    and    S    .  =  d     9  /d: 

oj  2     J  oj  1     j 

«o^  =  dz  e/d  Poi  =  '^i  Q/^ 

"02  =  -01  +  ^^1^1^2/<'^1^)  P02  =  Pqi  +  ^^1  ^/^ 

^^03  =  '^Ol  +  ^^22  ^2/^  P03  =  PqI  +  ^'^22^2^/<'^2^) 

^^04  =  °'02  +  ^^03  -  '^Ol  P04  =  P02  +  ^03  "  ^01  ^^-  3^) 

Two  special  cases  will  be  described  for  which    (III.29)    and    (111,31)    may  be  simplified. 

First,   assume  that  each  reflecting  surface  may  be  considered  a  plane.     Let     h      and     h 

t  tm 

be  the  heights  of  the  transmitting  antenna  and  the  knife  edge  above  the  first  plane,    and  let 

h  and     h       be  the  heights  of  the  knife  edge  and  the  receiving  antenna  above  the  second 

rm  r 

reflecting  plane.     Assume  that     Ar     is  very  small  for  every     A,     In  terms  of  the  heights     h  , 

h       ,     h       ,     h   ,      the  parameters     9,      d,,      and     d.     and  the  parameter     d     =  d,  d^/(2d): 
tm         rm         r  12  r  1     2 


A,     =  2h  h       /d,  ,     A-     =  2h    h        /d- 
2r  t    tm      1  3r  r    rm      2 

A,  -  d     9^,     A,  =  d    (9  +h        A,    )^ 
1         r  2         r  tm      2r 

A_  =  d    (9  +  h         A_    )^,     A^  3  d    (9  +  h^      A,     +  h         A_    )^  (in.  37) 

3         r  rm      3r  4         r  tm      2r  rm      3r  *  ' 

The  second  special  case  assumes  a  knife  edge  over  an  otherwise  smooth  earth  of 

effective  radius     a,   with  antenna  heights     h      and    h       small  compared  to  the  height  of  the  knife 

edge.     In  this  case,     h      and    h       are  heights  above  the  smooth  curved  earth.      The  angle  of 

elevation  of  the  knife  edge  relative  to  the  horizontal  at  one  antenna  is     9        and  relative  to  the 

ht 

hr 
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^zr-\\y'L  +^V(^^)  +  \t] 


3r 


e^^+4h^/(3a)  +  e^^ 


(in.  38) 


For  this  special  case,   the  forrmilas  (III. 29)    for     A.     may  be  simplified  by  writing 

e,     =  aJ     d/(2hj,     e.     =  A.^    d    /(2h    ) 
2r  2r      1  t  3r  3r     2  r 


(in.  39) 
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BEYOND-HORIZON    KNIFE-EDGE 

DIFFRACTION     WITH      GROUND 

REFLECTIONS 

(KNIFE-EDGE   NORMAL  TO   RAY  PATH) 

d    


—  d,    — H 
* — d,o — A 


Figure    111.9 
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III.  4    Parameters     K     and     b        for  Smooth  Earth  Diffraction 

o 

In  section  8,    the  parameters     K     and     b        are  shown  on  figures     8.1    and    8.2    for 
horizontally  and  vertically  polarized  waves  for  poor,    average,    and  good  groixnd,    and  for  sea 
water. 

Assiaming  a  homogeneous  ground  in  which  the  relative  dielectric  constant    e      and 
conductivity     cr      of  the  ground  are  everywhere  constant,      K     and     b     are  defined  as  follows: 


For  horizontal  polarization. 


-2  „      .-4     r,       ,.2  2      "^^ 


K     =   1.7778  X    10"     C     f"^     [(e-1)     +x    ]  (in.  40a) 

h  o  ■- '         '  •■ 


b     -  180°    -  tan"       (^^-^ j         degrees  .  (111.400) 


For  vertical  polarization, 


2         Z    Vz 
K     =  (e     +x   )  ^    K^  (in.  41a) 

v  h 


b     =  2  tan"     (e/x)  -  tan"      (^~      )         degrees.  (LEI.  41b) 


where     x     depends  on  the  gro\and  conductivity     cr      in  mhos  per  meter  and  the  radio  frequency 
f     in  megacycles  per  second,   and  has  been  defined  by  (III.  7)  as 

X  =  1.8  X  10^  (r/f 


C       is  defined  in  section  8  as 
o 

C     =  (8497/a)^ 
o 

where     a     is  the  effective  earth's  radius  in  kilometers. 

_4 
When     (r/f    >>  (e/2)  X  10        ,   the  parameters     K     and    b"  may  be  written  as 


.    1.325  X  10"^  C     f^/^-^/^      , 


K     ~    1.325X10       C     r       0-  ,  b,    ^    180  (in.  42) 

h  o  h  ' 

1/2     -5/6  o 

K    £=i    2.385  C     (T         f      '       ,         b     =^  0  (in.43) 

V  O  V 

and  when      (r/f  <   <  (e  /2)  X  10       ,     the  parameters     K     and    b     may  be  written  as 

K     ~    1.7778  X  10"^  C     f"^  (e  -  1)"^^    ,         b     --90°  (IH.  44) 

h  oh 

K^=:i€K^,         b^=^90°  (in.  45) 
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III.  5    Forward  Scatter 

The  attenuation  fiinction     F(ed)     for     N     =250,     301,     350,     and    400,     shown  in 
figure    9.1,     may  be  used  for  most  land-based  scatter  links.     When  a  path  is  highly  asym- 
metrical,  the  attenuation  for  a  given  value  of     0d     is  less  than  it  would  be  for  a  symmetrical 
path.     Figures    III.  1 1    to    III.  14    show  the  function     F(0d)     for  values  of     s     from    0.01    to     1, 
for  each  of  four  values  of     N   .     For  values  of     0d   <    10,     the  effect  of  asymmetry  is  negligible, 
but  increases  with  increasing     9d,   particularly  when     s  <    0.5. 

For  values  of     s     between    0.7    and    1,   the  function     F(9d)     may  be  computed  as 
follows: 


where 


■where 


where 


0.01  <   6d  <    10,         F(9d)  =  K(N   )  +  O.340d  +  30  log(ed)  (in.46) 


K(N  )  =  139.5  +  0.06  N     -  2.4  X  10~^  N  ^ 

s  s  s 


10  <  9d  <   70  ,         F(ed)  =  K  (N    )  +  K   (N   )ed  +  K   (N   )  log(ed)  (in.  47) 

OS  X        s  ^        s 


K„(N   )  =   140.7  +  0.06  N     -  3.2  X  10   ^  N  ^ 
0      s  s  s 

K,(N   )  =  3.925  +  3.3  X  lO"^  N     -  1.9  log  N 
Is  s  "      s 

K^(N   )  =  22.5  +  0.05  N 
2      s  s 


ed  >  70  ,         F(9d)  =  K(N    )  +  0.158  ed  +  44  log  (0d)  (in.  48) 


K(N   )  =   123.6  -    0.008    N 


Using  equations    (111.46)    to    (111.48),   the  function     F(9d)     may  be  computed  for  any 

[     N   .     Fc 
s 

following  values: 


value  of     N   .     For  the  values  of     N       shown  in  figures    IIT.  1 1    to    III.  14,     K(N    )     has  the 
s  s  ^  s 


N 
s 

250 

301 

350 

400 

0.01  <  9d  <  10 

K(N^) 

139.5 

135. 82 

131.  10 

125.  10 

10  <  ed  <  70 

V^s) 

135.70 

129.77 

122.50 

113.50 

K^(N^) 

0.  194 

0.212 

0.246 

0.  301 

K^(N^) 

35 

37.5 

40 

42.5 

9d  >  70 

K(N^) 

121.6 

121.2 

120.  8 

120.4 
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The  frequency  gain  fionction,      H   ,      for  the  special  case     h       =  h  frequently  used 

in  systems  design,    is  shown  as  a  function  of     r     on  figures    III.  15    to    III.  19    for      vj     =   1,     2, 

3,     4,      5,     and  for      s  =  0.75     to     1,     0.5,      0.25     and     0.1.     In  this  case,   no  correction 

factor     AH       is  required, 
o 

The  fixnction     H       for      -q     =  0,     shown  on  figure    9-5  corresponds  to  the  assumption 
of  a  constant  atmospheric  refractive  index.     This  function  may  be  computed  as  follows; 

2  ,       2 
2(1 -h        /h        ) 

re       te 


H  (t)     =  0)  =   10  log  ^ I  (in.  49) 

^    r2'[h(r^)-h(r2)]J 

/here     r,   =  4tt  9  h     /\    ,  r^  =  4iTeh      /\ 

1  te  '2  re       ' 

h(r^)  =  r^  f(r^)    ,         f(r  ^)  =  Ci(r  ^j  sin  r  ^  +   [  u/2  -  Si(r  ^)  ]   cos  r  ^  (IH.  50) 


and 


Mr^)  =  r^  i{T^)   .         i{T^)  =  Ci{T^)  sin  r^  +  [^/2  -  Si{T^)  ]   cos  r^ 


r  r 

C       cost     ,  „.,    ,         C       sint 


Ci(r)=     \       — -—    dt   ,         Si(r)  =     \      — -—  dt  (in.51) 

^00  ^  "'  0 

Values  of  the  sine  integral     Si{r)     and  the  cosine  integral     Ci(r)     for  arguments  from 
10    to     100    are  tabulated  in  volume  32  of  the  U.   S.   NBS  Applied  Math  Series  [1954].     See 
also  [NBS  AMS  1964].     The  function     h(r)     is  shown  graphically  in  figures    111.20    and    III.21. 

For  the  special  case  of  equal  effective  antenna  heights,      h      =  b     ,     equation    (111.49) 
is  not  applicable.     In  this  case     H  {-q     =0)     is  computed  as; 

H  An    =  0)  =  10  log    I  -J \  (in.  52) 

L  r''[h(r)-rg(r)]    J 

where 

g(r)  =  Ci(r)  cos  r  -  [it/2  -  Si(r)  ]   sin  r  (IH.  53) 

When  the  effective  height  of  one  antenna  is  very  much  greater  than  that  of  the  other, 
the  computation  may  be  simplified  as  follows: 

For     r    <  <  r  H   (t|     =  0)  =   10  log    J— I.  (in.  54a) 

^'[l-h(r^)] 

For     r2>>r^,        H^tj^  =  0)  =   10  log    i  — ^ j.  (lU.  54b) 


^1    [  1  -  h(r^)  ] 
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THE    FUNCTION     F(0d)    FOR    Ns=  250 
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III.  6     Transmission  Loss  with  Antenna  Beams  Elevated 
or  Directed  Out  of  the  Great  Circle  Plane 

The  methods  of  section  9  may  be  modified  to  calculate  a  reference  value  of  long-term 
median  transmission  loss  v/hen  antenna  beams  are  either  elevated  or  directed  away  fronn  the 
great  circle  path  between  antennas.      For  many  applications,    the  average  transmission  loss 
between  antennas  with  random  relative  orientation  is  about     10  db    more  than  the  basic  trans- 
mission loss,   which  assumes  zero    db    antenna  gains. 

Figure    111.22    shows  scattering  subvolumes  at  intersections  of  antenna  main  beams 
and  side  lobes.     A    "scatter"    theory  assumes  that  the  total  power  available  at  a  receiver  is 
the  sum  of  the  powers  available  from  many  scattering  subvolumes.      For  high  gain  antennas, 
the  intersection  of  main  beams  defines  the  only  important  scattering  volume.     In  general,    all 
power  contributions  that  are  within     10    db    of  the  largest  one  should  be  added. 

For  a  total  radiated  power     p    ; 

-0.1  L  -0.1  L. 

p    /p     =   10  ^^    ,         p    ./p     =   10  ^  (m.  55) 

at  ai     t 


where      L  is  the  transmission  loss  and     L.     is  the  loss  associated  with  the     ith     power 

sr  1 

contribution,     p    .  : 
ai 

^         -0.1  L. 

L       =     -10  log  (p    /pj  =  -10  log     >     10  ^  (nL56) 

sr  a      t  (_i 


L.  =  30  logf  -  20  log(d^/r    )  +  F(e.d)-F.+H.+A     -G.-G.+L.  (HI.  57) 

1  '=  "  o  ei  oi  oi  a  ti  ri  gi 

In    (111.57)    f,     d,     and     A       are  defined  as  in    (9.1)    and  the  other  terms  are  related  to 

similar  terms  in    (9.1).     If  the  effective  scattering  angle     9    .     for  the    ith    intersection  is 

equal  to  the  minimum  scattering  angle     0,   then     F(8    .d  ),     F   .,     H   .     are  equal  to     F(9d), 

ei  oi         oi 

F   1      and     H   ,     and     G  .  +  G   .  -  L    .  =   G    .      Note  that  a  term     20  log(r    /d)     has  been  added 
o  o  ti  ri  gi  p  °  *    o      ' 

in    (III.57)    to  provide  for  situations  where  the  straight  line  distance     r       between  antennas 
is  much  greater  than  the  sea-level  arc  distance     d.     Such  differences  occur  in  satellite  com- 
munication. 

Scattering  planes,   defined  by  the  directions  of  incident  and  scattered  energy, may  or 
may  not  coincide  with  the  plane  of  the  great  circle  path.     Each    "scattering  plane"    is 
determined  by  the  line  between  antenna  locations     and  the  axis  of  the  stronger  of  the  two 
intersecting  beams,   making  an  angle     t,     with  the  great  circle  plane. 


m-37 


The  free  space  directive  gain  patterns  of  the  antennas  are  replaced  by  equivalent  values 

for  ease  in  computation.     For  an  idealized  pencil-beam  antenna  with  a  half-power  beamwidth 

2 
25     and  a  circular  beam  cross-section,   the  directive  gain     g     is     4/5    >    assuming  that  all  of 

the  power  is  radiated  through  the  main  beam  and  between  the  half -power  points.     An  equiva- 
lent beam  pattern  with  a  square  cross -section  and  a  semi -beamwidth     5       has  a  gain  of 

Tr/5    ,     thus      5     =  6*^  tt/4,     and  the  maximum  free  space  gains  are 
o  o 

G^  =  10  log  g^     =  4.97  -  10  log  6^        5^  db  (m.  58a) 

t  to  two     tzo  ' 

G    =  10  log  g       =  4.97  -  10  log  6  5  db  {HI.  58b) 

r  ro  rwo    rzo  ' 


where  the  subscripts     w     and     z     refer   to  azimuthal  and  vertical  angles.     In  most  cases, 

6  and     5  may  be  replaced  by  their  geometric  mean,      5     =  (5       5      )      .     The  free  space 

wo  ZO  JO  Q  ^Q     ^Q 

directive  sain  of  a  main  beam  may  be  measured  or  approximated  as     g   ~  ir/CZ  5       5      ).     Gains 

o  wo    ZO 

for  side  lobes  are  determined  from     g        and  the  ratios      e,  /g    ,     g^/g    ,    .  .  .  ,     which  may  be 

°o  1      o         Z      o 

measured  or  calculated.  The  average  gain  g  for  other  directions  depends  on  the  fraction 
of  power  radiated  in  those  directions.  For  instance,  if  half  the  total  power  of  a  transmitter 
is  radiated  in  these  directions,   and  if  the  polarization  coupling  loss,      L      ,      is     3  db,     then 

G^     -  L       =     -  6    db 
bt  cp 

since  the  definition  of  the  directive  gain,    G,    ,     assumes  for  every  direction  the  receiving 

antenna  polarization  appropriate  for  maximum  power  transfer. 

Figure    III.23    shows  an  antenna  power  pattern  in  several  different  ways,    including  a 
Mercator  projection  of  the  surface  of  a  unit  sphere. 

The  plane  that  determines  the     "bottom"    of  a  beam  is  perpendicular  to  the  great 
circle  plane  and  forms  an  angle     4".     with  a  horizon  plane: 

4^^-  =  K^-  -  9  .  .       4j  .  =  9,    .  -  e  (m.  59) 

ti         bti         et  ri         bri         er 

where  9  is  the  angle  of  elevation  of  the  lower  half -power  point  of  a  beam  above  the  hori- 
zontal, and  0  is  defined  in  section  6.  If  an  antenna  beam  is  elevated  sufficiently  so  that 
ray  bending  may  be  neglected,    the  angles      a       and     p        are  denoted     a  and     p       : 

a       =  (a       +  4/  )  sec  ^    ,         p       ^  (^       +  qj    )  sec  t,  (IH.  60) 

eo  oo         t  eo        '^oo         r 

where     4     is  the  angle  away  from  the  great  circle  plane.     The  angles     a         and     p  are 

defined  as  in  section  6  using  the  actual  radius,      a     ii^    6370  km,     instead  of  an  effective  radius 
a . 


IU-3J 


When  ray  bending  must  be  considered,  the  equations  for     a       and     (3        are 

d^  ^sec  t,  ^  ,  J  c 

Lit  ^.    A  f  ^  d  sec  L, 


a     =  a       +  T  (    e^    . ,     N     )    -Tie,    5 ,    N      )  (IE.  6la) 

e         eo  \    bt  2  s/  \    bt  Z  s/ 


P.  =P..  +  nV      -^^'     ^)-^(^br'    ^^'    ^s^  (™-^l^) 


where     """(B    ,    d,    N   )     is  the  bending  of  a  radio  ray  which  takes  off  at  an  angle     6       above  the 

horizontal  and  travels     d     kilometers  through  an  atmosphere  characterized  by  a  surface 

refractivity     N   .     The  ray  bending     t     may  be  determined  using  methods  and  tables  furnished 
s 

by  Bean  and  Thayer    [1959  ]•     For  short  distances,      d,      or  large  angles,      6    ,     t     is  neglig- 
ible.    If     9       is  less  than    0.1    radians,   the  effective  earth's  radius  approximation  is  adequate 
for  determining     t  , 

'b'     ^^^'^s)=^     [l-a^/a(N^)]  (in.  62) 

o 

The  reference  value  of  long-term  median  transmission  loss     L  is  computed  using 

(111.56)    where  the  losses  associated  with  several  scattering  subvolumes  are  computed  using 

(111.57).     The  attenuation  function     F(d9    .)     is  read  from  figure    9.1    or  figures    111,11  -111.14 

as  a  fimction  of     9    .. 
ei 

The  generalized  scattering  efficiency  term     F   .     is 

F   .=   1.086  (ti      /h   )  (2  h     -  h,    -  h     -h^^-h^     )        db  (ni.  63) 

oi  'se      e  o         1         e         Lt         Lr 

where 

Q     =  a     +  p     ,         s=a/p,         h=sde/(l+s),        t|       =Ti(h,N)        (III64) 
e         6       ^e  e         e    "^e  e         e         e  e  "se        's      e        s  ^xj-j-.  ot; 


and  the  other  terms  are  defined  in  section  9.     In  computing  the  frequency  gain  function     H   ., 

if        a     >  a       use     r,  =    oo,      if     B     >  (3       use     r     =  oo;     then     H   .  =  H     +3  db.     If  both  antennas 
eo  1  eo  2  oio 

are  elevated  above  the  horizon  rays,    H   .  =  6  db.      Atmospheric  absorption    A       is  discussed 

oi  a 

in  section  3.     The  gains     G       and    G   .    are  the  free  space  directive  gains  defined  by   (111.58), 

and  the  loss  in  gain     L,    .     is  computed  as  shown  in  section  9  replacing    -q   ,     s,      5,      and     9 
gi  s 

by     TI      ,     s    ,     5    ,      and     9    . 
'       'se         e         e  e 
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In  computing  long-term  variability  of  transmission  loss  for  beams  elevated  above  the 

horizon  plane,   the  estimates  of     V     and     Y     given  in  section  10  should  be  reduced  b-jr  the 

factor     f(9    )     shown  in  figure    111.24,    with    9     =  9    • 
h  h         b" 

V   (50,  d    )  =  V(50,  d    )  f(9    )  (HL  65a) 

e  e  e  n 

Y^(P.  d^)  =  Y(p,  d^)  f(e^)  (in.65b ) 

The  angle     9       used  in    (111.65)    should  be  the  elevation  above  the  horizontal  of  the 
scattering  subvolume  corresponding  to  the  minimum  value  of     L.  . 
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SCATTERING    SUBVOLUMES    IN   A   SCATTERING    PLANE 


SIDE  LOBE 


SIDE  LOBE 


MAIN    BEAM 


MAIN   BEAM 


SIDE   LOBE 


D.  =  2B 


ILLUSTRATION    OF   A    SCATTERING    PLANE    CONTAINING 
THE    MAIN    BEAM    OF   THE    RECEIVING    ANTENNA 


FREE    SPACE    ANTENNA    PATTERN 


INTERSECTION    OF   AN    EQUIVALENT 
ANTENNA    POWER    PATTERN    WITH 
THE    SURFACE   OF  A   UNIT    SPHERE 


ANTENNA    VOLTAGE    PATTERN. 

POLAR    DIAGRAM    IN    THE 

EQUATORIAL    PLANE 


POLAR    AXIS 


POLAR    PLANE 

THROUGH 

MAIN   BEAM   AXIS 


PATTERN    PARAMETERS 

8^0=0.04 

Stwo^ooe 

Sjzi  =0.02 

Stwl  =0.02 

S|22  =  0.02 

Stw2  =  0.02 

^tzr°'2 

e,     ,  =0.2  0 

twi 

€,22  =  0.22 

Mw2=0.30 

gto  =654.5 

Gfo    =28.l6db 

g,|    =20.7 

G,|     =  I3.l6db 

g, 2  =3.28 

G,2    =    5.l6db 

g^,  =0.395 

Gbt    =  -4.04db 

CARTESIAN  DIAGRAM  OF  TRANSMITTING  ANTENNA  POWER    PATTERN  IN  THE    EQUATORIAL   PLANE 
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III.  7     Long-Term  Power  Fading 
Long-term  power  fading  is  discussed  in  section  10.      Figures   10.  5  to  10.  16  show  em- 
pirical estimates  of  all-year  variability  for  (1)   continental  temperate  (2)   maritime  temperate 
overland  and  (3)   maritime  temperate  oversea  climates.      The  curves  shown  on  these  figures 
are  based  on  a  large  amount  of  data.      Estimates  of  variability  in  other  climates  are  based  on 
what  is  known  about  meteorological  conditions  and  their  effects  on  radio  propagation,    but  have 

relatively  few  measurements  to  support  them. 

Figures     III.25    to    111.29     show  curves  of  variability  relative  to  the  long-term  median, 

prepared  by  the  CCIR  [  1963  f  ]     for  the  following  climatic  regions: 

(4)  Maritime  Subtrdpical,  Overland. 

(5)  Maritime  Subtropical,    Oversea. 

(6)  Desert, 

(7)  Equatorial. 

(8)  Continental  Subtropical. 

In  some  cases,    random  path  differences  have  ixndoubtedly  been  attributed  to  climatic  dif- 
ferences.    Available  data  were  normalized  to  a  frequency  of    1000    MHz,     and  the  curves 
correspond  to  this  frequency.     They  show  all-year  variability     Y   (p,  d    )     about  the  long-term 
median  as  a  function  of  the  effective  distance     d       defined  by    (10.3),     Variability  estimates 
for  other  frequencies  are  obtained  by  using  the  appropriate  correction  factor     g(f)     shown  in 
figure    in.  30: 

Y{p)  =  Y^(p,  d^)  g(f)  (in.  66) 

The  empirical  curves     g(f)     are  not  intended  as  an  estimate  of  the  dependence  of  long- 
term  variability  on  frequency,   but  represent  an  average  of  many  effects  that  are  frequency- 
sensitive,   as  discussed  in  section  10. 

Variability  about  the  long-term  median  transmission  loss     L(50)     is  related  to  the 
long-term  reference  median     L         by  means  of  the  function     V(50,  d    )     shown  on  figure     10,1. 
The  predicted  long-term  median  transmission  loss  is  then: 

L(50)  =  L       -  V(50,  d    )  (in.  67) 

cr  e 


and  the  predicted  value  for  any  percentage  of  time  is 

L(p)  -  L(50)  -  Y(p)  (in.  68) 

as  described  in  section  10. 
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in.  7.  1    Diurnal  and  Seasonal  Variability  in  a  Continental  Temperate  Climate 
The  curves  shown  in  figures     10.5    to     10.16    and    111.25    to    111.29    represent  variability 
about  the  long-term  median  for  all  hours  of  the  day  throughout  the  entire  year.      For  certain 
applications,   it  is  innportant  to  know  something  about  the  diurnal  and  seasonal  changes  that 
may  be  expected.     Such  changes  have  been  studied  in  the  continental  United  States,   where  a 
large  amount  of  data  is  available.     Measurement  programs  recorded  VHF  and  UHF  trans- 
mission loss  over  particular  paths  for  at  least  a  year  to  determine  seasonal  variations.     Data 
were  recorded  over  a  number  of  paths  for  longer  periods  of  time  to  study  year-to-year 
variability. 

As  a  general  rule,   transmission  loss  is  less  during  the  warm  summer  months  than  in 
winter,   and  diurnal  trends  are  usually  most  pronoionced  in  summer,   with  maximum  trans- 
mission loss  occurring  in  the  afternoon.     The  diurnal  range  in  signal  level  may  be  about 
10  db    for  paths  that  extend  just  beyond  the  radio  horizon,   but  is  much  less  for  very  short  or 
very  long  paths.     Variation  with  season  usually  shows  maximum  losses  in  mid-winter, 
especially  on  winter  afternoons,    and  high  fields  in  summer,   particularly  during  morning 
hours.     Transmission  loss  is  often  much  more  variable  over  a  particular  path  in  summer 
thanit  is  during  the  winter,    especially  when  ducts  and  elevated  layers  are  relatively  common. 
The  data  were  divided  into  eight    "time  blocks"    defined  in  table   III.  1,     The  data  were 
assumed  to  be  statistically  homogeneous  within  each  of  the  time  blocks.     With  more  and 
shorter  time  blocks,   diurnal  and  seasonal  trends  would  be  more  precisely  defined,   except 
that  no  data  would  be  available  in  some  of  the  time  blocks  over  many  propagation  paths.     Even 
with  the  division  of  the  year  into  winter  and  summer  and  the  day  into  four  periods  as  in  table 
UI.  1,   it  is  difficult  to  find  stifficient  data  to  describe  the  statistical  characteristics  expected 
of  transmission  loss  in  Time  Blocks  7  and  8. 

Table  in.  1 
Time  Blocks 

No.  Months  Hours 

1  Nov.  -  Apr.  0600  -  1300 

2  Nov.  -  Apr.  1300  -  1800 

3  Nov.  -  Apr.  1800  -  2400 

4  May  -  Oct.  0600  -  1300 

5  May  -  Oct.  1300  -  1800 

6  May  -  Oct.  1800  -  2400 

7  May  -  Oct.  0000  -  0600 

8  Nov.  -  Apr.  0000  -  0600 

In  some  applications,   it  is  convenient  to  combine  certain  time  blocks  into  groups,   for 
instance,   some  characteristics  of  long-term  variability  are  significantly  different  for  the 
winter  group    (Time  Blocks     1,    2,    3,    8)  than  for  the  summer  group    (Time  Blocks    4,    5,    6,    7). 


Ill -4  5 


In  other  climatic  regions,    if  the  annual  range  of  monthly  average  values  of     N        is  less 
than    20  N    units     (figure    111.31),     seasonal  variations  are  expected  to  be  negligible.      One 
would  also  expect  less  diurnal  change,   for  example,   in  a  maritime  temperate  climate  where 
changes  in  temperature  during  the  day  are  less  extreme.     In  climates  where     N       changes 
considerably  throughout  the  year,   the  consecutive  4-6  month  period  when     N       is  lowest  may 
be  assumed  to  correspond  to    "winter",   whatever  months  maybe  involved. 

The  parameter     V(50,    d    )    for  each  of  the  eight  time  blocks  and  for    "summer"  and 
"winter"  is  shown  in  figure    III.  32.      Curves  of  the  variability    Y    (p,  d    )    about  the  long-term 
median  for  each  of  these  times  of  day  and  seasons  are   shown  in  figures     III.  33  to    III.  42. 
These  curves  are  drawn  for  a  frequency  of  100   MHz  .      Figures    III.  33     and    III.  34     show  the 
1  to  99  percent  range  of    Y    (p,  d    )    for  the  winter  time  blocks,    1,    2,    3,    8  and  the  summer 
time  blocks  4,    5,    6,    7.     Each  group  of  data  was  analyzed  separately.     Some  of  the  differences 
shown  between  time  blocks   1,    2,    3  and  8  are  probably  not  statistically  significant.     Marked 
differences  from  one  time  block  to  another  are  observed  during  the  summer  months. 

The  curves  for  summer,    winter,    and  all  hours  shown  in  figures   10.  1  through  10.  10 

represent  a  much  larger  data  sample.     For  a  number  of  the  paths  time  block  information  was 

not  available. 

The  smooth  curves  of    V(50,    d    )    and    Y    (p,  d    )    versus    d       shown  in  figures   10.  1, 

e  o  e  e 

10.  2,     and    III.  3  2    to    III.  42    may  be  represented  by  an  analytic  function  of  the  general  form; 


V{50) 


n  n 

Yjl0)j=[c^d^     -f,(dj'l    exp(-  c,d^     )+f^{dJ  (in.69) 


o'      '    r         I    ^1     e  2^    e'  "  '        3     e     '    '     2^    e 


-Y    (90) 
o 


where 


r  ""2 

f^(d    )  =  f      1  1   -  (1   -  f     /f    )  exp  (  -  c,  d       )    I  (in.  70) 


The  terms      c    ,      c^,     c^,     n    .     n^,     n.,     f     ,      and     f        in    (111.69)    and    (III.70)    are 
12312Jm  00 

constants  for  any  given  time  block  and  value  of     p.      The  parameters     f         and     f        are  max- 
innum  and  asymiptotic  values,    respectively.  Tables     III. 2    to    111.4    list  values  of 

the  eight  parameters  required  in    (111.69)    to  obtain     V(50,  d    ),      Y   (10,  d    ),      and     -Y    (90,  d    ) 
for  the  eight  time  blocks  in  table  III.l,     and  for  summer,   winter,    and  all  hours.      The  constants 
given  in  Tables  III.  2  to  III.  4  for  summer,    winter  and  all  hours  were  determined  using  only 
radio  paths  for  which  time  block  information  is  available.      They  do  not  yield  the  curves  shown 
in  figures   10.  1  and  10.  2  of  section  10,    which  represents  a  much  larger  data  sample. 
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in.  7.  2    To  Mix  Distributions 

When  a  prediction  is  required  for  a  period  of  time  not  shown  on  the  figures  or  listed 
in  the  tables,   it  may  sometimes  be  obtained  by  mixing  the  known  distributions o     For  example, 
the  distributions  for  time  blocks  5  and  6  would  be  mixed  if  one  wished  to  predict  a  cumulative 
distribution  of  transmission  loss  for  summer  afternoon  and  evening  hours.     In  mixing  dis- 
tributions,   it  is  important  to  average  percentages  rather  than  levels  of  transmission  loss. 
Distributions  of  data  for  time  blocks  may  also  be  mixed  to  provide  distributions  for  other 
periods  of  time.     For  example,    data  distributions  for  time  blocks    1,    2,    3,     and    8    were 
mixed  to  provide  distributions  of  data  for    "winter".     When  averages  are  properly  weighted, 
such  mixed  distributions  are  practically  identical  to  direct  cumulative  distributions  of  the 
total  amount  of  data  available  for  the  longer  period. 

The  cumulative  distribution  of     N     observed  hourly  median  values  is  obtained  as 

follows:    (1)    the  values  are  arranged  in  order  from  smallest  to  largest,     L^,     L,_,     L.   , , 

L   , ,    L    ,     (2)    the  percentage     p     of  hourly  median  values  less  than     L       is  computed: 

(3)    a  plot  of     L,       versus      100  p(n)     for  values  of    p     from    1/{2N)     to     1  -  1/(2N)     is  the 
n 

observed  cumulative  distribution. 

To  mix  two  distributions,   the  following  procedure  is  used:    (1)    choose  ten  to  fifteen 
levels  of  transmission  loss     L   , ,    L,   ,     covering  the  entire  range  of     L{p)     for  both  dis- 
tributions,    (2)    at  each  of  these  levels,    read  the  value     p     for  each  distribution  and  average 
these  values,     (3)    plot  each  selected  level  of  transmission  loss  at  the  corresponding  average 
percentage  to  obtain  the    "mixed"    distribution.     In  this  way,   any  number  of  distributions 
may  be  combined,   if  each  of  them  represents  the  same  number  of  hours.     In  general,   a 
weighted  average  percentage     p     should  be  computed,   using  as  weights  the  number  of  hours 
represented  by  each  distribution. 
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Yo  (p.de)  IN  DECIBELS 
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in.  8     Examples 
Sample  calculations  are  given  here  to  show  the  application  of  the  prediction  methods 
described  in  this  report  to  specific  problems.      The  first  example  shows  predictions  for  a 
line-of-sight  path  over  irregular  terrain.      The  second  shows  predicted  transmission  loss  over 
a  rounded  isolated  obstacle.      The  predicted  cumulative  distribution  of  hourly  median  values 
is  compared  with  measurements  over  this  path.      The  third  example  computes  predicted  scatter 
and  diffraction  losses  over  a  transhorizon  path  and  compares  the  predictions  with  long-term 
recordings  of  data. 

III.  8.  1     Line-of-Sight  Predictions 

Attenuation  relative  to  free  space  is  predicted  for  a  short  line-of-sight  path  shown 
in  figure  III.  43.     Measurements  at  a  frequency  of  100  MHz  were  made  using  vertical  polari- 
zation.     The  transmitting  and  receiving  antennas  are  4  meters  and  9  meters,    respectively, 
above  ground. 

A  straight  line  is  fitted  by  least  squares  to  the  terrain  visible  from  both  antennas. 
Terrain  near  the  transmitter  is  excluded  because  it  is   shadowed  by  high  foreground  terrain. 
Twenty-one  equidistant  points    x.  =  x   ,    x   ,    .    .    .   x  are  chosen  as  shown  onfigure  III.  43a, 

and  the  corresponding  terrain  heights,      h.  ,    are  read.     From  (5.  15)   the  average  terrain  height 
h    is   1531.8  m,    the  average  distance    x    is     13.0  km,    and  the  slope    m    is     -6.0  meters  per 
kilometer.      The  equation  for  the  straight  line  is  then 

h(x)  =  1531.8  -  6(x-13)  m  =   1.  5318  -  6(x-13)  •    lO"'^  km. 

An  effective  earth's  radius,      a,    is  obtained  using  figure  4.  1  and  equations  (4.3)   and 

(4.4).     For  this  area  in  Colorado    N       is  280  and    a  =  8200  km.     From  (5.  16)   the  adjustment 

s 

to  allow  for  the  sea  level  curvature  is 

2 
y(x)  =  h(x)    -  x    /(16,  400)     km 

Figure  III.  43b  shows  the  curve    y(x)  vs  x    and  terrain  which  has  been  modified  to  allow  for 
the  sea  level  curvature. 

At  the  transmitter,      x  =  0     and    h(x  -  0)   is   1609.  5  m.     At  the  receiver,      x  =  d  =   19.  75 
km  and    h(x  -   19.75)  is   1491.4  m.     From  (5.  17)   the  antenna  heights  above  the  smooth  reflecting 
plane  are  then: 

h'    =  h       -  h(0)  =   1647.  1   -   1609.  5  =  37.  6  m  =  0.0376     km, 
1  ts 

h'    =  h        -  h(d)  =   1524.0   -   1491.4  =  32.  6  m  =  0.0326    km, 
2  rs 

where    h       =   1647.  1  m    and    h        =   1524.0  m    are  the  heights  above  sea  level  at  the  transmitter 
ts  rs  ^ 

and  receiver  respectively.     At  100  MHz  (X.  =  3  m),    the  criterion  that  both    h'     and    h'      must 
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exceed    0.  16  A.    is  met.     Neither    h'      nor    h'      exceeds  one  kilometer,    so  no  correction  factor. 
Ah,    is  required.     From  (5.  6)  and  (5.  7)   the  distances     d       and    d       from  each  antenna  to  the 
point  of  specular  reflection  are 

d     =  d(l  +  h'  /h')"^  =   10.  58  km,      d     =  d(l  +  h'  /h')"-^  =  9.  17    km. 


and  the  grazing  angle     i|/    is 


tan  qj  =  h'  /d     =:  h'  /d     =  0,003  554 


i(j  =  0.003  554  radians. 
From  (5.9)   the  path  length  difference,  Ar  ,    between  direct  and  reflected  rays  is 

Ar  =  [d^  +  (h'    +  h')^ 


1 

2 

d^    +    (h'^ 

1 

21^                                 -4 
-  h^)          =  1,2413  X  10         km 

- 

_ 

The  approximation 

Ar  o^  2h'  h'   /d  =   1.2413  x  lO"'^  km  =  0.  124  m  =  0.04k 

is  also  valid  in  this  case.     Note  that    Ar     is  less  than    0,  12k    and  optical  methods  including  a 
divergence  factor  may  underestimate  the  attenuation. 

One  should  note  that  important  reflections  might  occur  from  the  high  ground  near  the 
transmitter.     In  this  case  the  reflecting  plane  would  correspond  to  the  foreground  terrain 
giving    h'    =  4  m,    h'    =   50  m,    d     =   1 .  53  km,    d     =   18.22  km  and    Ar  =  0.02  m    which  is  much 
less  than    0.  16\  .     Optical  methods  would  not  be  applicable  here. 

The  attenuation  relative  to  free  space  may  be  estimated  using  one  of  the  methods 

described  in  subsection  5.2.     Of  these,    methods  4  and  5  would  apply  in  this  case.     Choosing 

heights    h     =  4  m,    h     =  25  m,    as  heights  above  foreground  terrain,    the  theoretical  smooth 

earth  curves  in  the  CCIR  Atlas  [  1959]    show  the  predicted  field  to  be  about  36  db  below  the 

free  space  value.      The  "standard"  propagation  curves  on  figure  1.7,    drawn  for   100  MHz  and 

h     =  h     =30  meters  show  the  median  basic  transmission  loss  to  be  about  15  db  below  the 
1  2 

free  space  loss.     Greater  attenuation  would  be  expected  with  lower  antennas  over  irregular 
terrain.     Method  5  using  the  empirical  curve  through  data  recorded  at  random  locations  in 
figure  I.  1  shows  the  attenuation  to  be  about  20  db  below  free  space.      These  data  were  recorded 
with  an  average  transmitter  height  of  about  250  m,    and  a  receiver  height  of  10  m. 

For  the  very  low  antennas  used  on  this  Colorado  path  one  would  expect  the  losses  to 
exceed  the  values  shov/n  on  figures  I,  7  and  I.  1,    and  also  to  exceed  the  theoretical  smooth 
earth  value  of    A  =;  3  6  db    obtained  from  the  CCIR  Atlas.     Spot  measurements  yield  a  value  of 
about  40  db. 
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If  a  prediction  were  desired  for  transmission  over  the  same  path  at    300  MHz,  X.  =   1  m, 
then    Ar  =  0.  124  m    is  slightly  greater  than    0 .  12  k    and  optical  methods  may  be  used.     Using 
the  value     Ar  =  0,  124  m  the  path  length  difference  in  electrical  radians  is     2iTAr/\  =  0.  780  5 
radians.     As  a  check,    this  quantity  may  be  computed  using  (5.  14a): 

2tt  Ar/X.  =  41.917  fh'  h'  /d  =  0.  7805  radians 

=  44. 7  degrees 
Equation  (5.4)   shows  the  attenuation  relative  to  free  space  assuming  a  single  ground 

reflection  from  the  smooth  curve    y(x),    figure  III.  43b.     Assuming  that  free  space  gains  are 
realized  so  that    G     =   10  log  g      g         the  equation  may  be  written 

p  Ol      02 


A  =  -10  log 


2  /2iTAr 

1  +  R      -  2R     cosf-^^— ^  -  c 
e  e  \      K 


where    R       is  the  effective  reflection  coefficient  defined  by  (5.  1): 
e 

a      g      .  -T               -0.  6  cr     sin  \\i 
R     .  DR/fn^t    -pf T ^^. 

With    f  =  300  MHz,    and  tan  i\i  -  0.003  554,    figure  III.  3   shows  the  theoretical  reflection  coef- 
ficient   R  =  0.97  and  the  phase  shift    c  =  0     for  vertical  polarization  over  average  ground. 
The  angle  between  the  direct  and  the  reflected  ray  is  small  so  the  ratio  of  gains  in  (5,  1)   may 

be  considered  to  be  unity.     The  divergence  factor    D    and  effective  reflection  coefficient    R 

e 

are 


B  =   fl  + 


'^^2  i 


(  i  +  —r- --I   '  =  0.865 

\.  ad  tan  4>y 


R    =    0.839  expi 


-0.6  0"     sinij'^ 
e       ""^\^  X^  )' 

The  terrain  roughness  factor,      cr    ,    is  the  root-mean-square  deviation  of  nnodified 

h 

terrain  elevations  relative  to  the  curve    y(x)     within  the  limits  of  the  first  Fresnel  zone  in 

the  horizontal  reflecting  plane.     The  first  Fresnel  ellipse  cuts  the  great  circle  plane  at  two 

points    X       and    x,      kilometers  from  the  transmitter.     The  distances    x       and    x,      may  be 
a  b  a  b 

computed  using  equations  (III.  18)   or  (III.  19)   to  (III.  21). 

B  =  0.135,     X     =   10.02,     X     =9.12 

X     =  x„   -  X,   =  0.90  km,     x,    =x„  +  x,   =   19.14  km 
a         0  1  b         0  1 

The  first  Fresnel  zone  cuts  the  great  circle  plane  at  points    0.9  and  19.  14  km    fronn 
the  transmitter  with  an  intervening  distance  of  18.24  km.      Equidistant  points  are  chosen  at 
X  =   1,    2,    ...    19  and  corresponding  modified  terrain  heights  and  values  of    y(x)     are  obtained. 
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With  height  differences  in  kilometers: 

19 


i'l'T 


h.)    /19,      0-     =  0.008222 


j=l 


The  effective  reflection  coefficient  is  then 


R     =  0.839  exp  -  0.01753  =  0.824 
e 


which  is  greater  than  0.  5  and  greater  than  »^sin  ij;  .     The  predicted  attenuation  relative  to  free 
space    A    is  then 


10  log 


2  2TTAr         \  i  r 

1  +  R^  -  2R^  cos        ^        "  '^  )  I        =  -  10  log  I  1.  6793  -  1.  6484  cos  0.  7805 


=^  3  db 


Due  to  diffraction  effects  over  irregular  terrain,    the  attenuation    A    is  often  ob- 
served to  be  much  greater  than  the  values  corresponding  to  the  ray  theory  calculations  il- 
lustrated in  this  example.     Ray  theory  is  most  useful  to  identify  the  location  and  depth  of 
nulls  in  an  interference  pattern  in  the  region  visible  to  two  antennas.     Figure  III.  44  shows  an 
interference  pattern  from  an  aircraft  at  10,000  ft.,    transmitting  on  3  28.2  Mc/s.  Measured 
values  compared  with  theoretical  curves  based  on  ray  theory  are  shown  on  the  figure. 
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III.  8.  2     Transmission  Loss  Prediction  for  a  Rounded  Isolated  Obstacle 
The  path  selected  to  provide  an  example  of  knife-edge  diffraction  calculations  is  lo- 
cated in  eastern  Colorado,    extending  from  a  location  near  Beulah,    southwest  of  Pueblo,    to 
Table  Mesa  north  of  Boulder.      The  common  horizon  is  formed  by  Pikes  Peak,    with  an  eleva- 
tion 4300  meters  above  mean  sea  level.     For  the  purpose  of  these  calculations  Pikes  Peak  is 
considered  to  be  a  single  rounded  knife  edge.      The  complete  path  profile  is  shown  in  figure 
III.  45.      Table  III.  5  gives  all  applicable  path  and  equipment  parameters  and  permits  a  compari- 
son of  calculated  and  actually  measured  values. 

TABLE  III.  5 
Path  and  Equipment  Parameters 
Carrier  Frequency 
Total  Great  Circle  Path  Distance 


Great  Circle  Distances  from  Pikes  Peak 
to  Transmitter  Site 
to  Receiver  Site 

Terminal  Elevations  above  Mean  Sea  Level 
Transmitter  Site 
Receiver  Site 

Elevation  of  Pikes  Peak 

Above  Mean  Sea  Level 

Above  Mean  Terminal  Elevation 

Transmitting  Antenna  Height  Above  Ground 

Transmitting  Antenna  Gain  Above  Isotropic 
(4.3  m  Dish) 

Receiving  Antenna  Height  Above  Ground 

Receiving  Antenna  Gain  Values  Above 
Isotropic  (3  m  Dish) 

Polarization 

Modulation 


751  MHz 
223.  3  km 

77.3  km 
146.0  km 

1,  905  m 

1,  666  m 

4,  300  m 

2,  50  7  m 

7.3  m 

26. 7  db 
20.0  m 

23 . 6  db 

Horizontal 

C  ontinuous 
Wave 

445  watts 


Transmitter  Power 

Calculations  are  given  for  single-ray  diffraction,    neglecting  possible  specular  reflections 

from  foreground  terrain. 

The  minimum  monthly  surface  refractivity    N       (referred  to  mean  sea  level)     from 
figure  4.  1  is    300  N-units.     From  Table  III.  5  the  terminal  elevations  are  190  5  and  1666  m, 
respectively.     Corresponding  surface  refractivity  values    N       are  245  and  251  N-units  (4.3), 


and  the  average  of  these  values  is     N 


248. 


In  this  example,      N       is  calculated  for  the 
s 


terminals,    as  the  antennas  are  more  than  1  50  m  below  their  4300  m  radio  horizon.      Using 
(4.4)   or  an  extrapolation  of  figure  4.  2,    the  effective  earth  radius     a    for    N     -   248     is  found 
to  be  7830  km. 
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The  angular  distance    9    in  radians  and  related  parameters  are  calculated  using  (6.  15) 

and  (6.  18a,  b):     9       =0.008581     9       =0.025953     a       =0.021827     (3       =0.041225.     In  this  ex- 
et  er  oo  oo 

ample     d        and    d         are  negligibly  small,    and  the  corrections     Aa       and    Ap       (6.  19a,  b)   can 
st  sr  o  o 

be  neglected.      Thus,     a       =   a   ,    p        =   P    .    and    9  =  0.0  630  52  radians, 
oo         o        oo  o 

The  free-space  loss  and  the  attenuation  relative  to  free  space  are  computed  con- 
sidering Pikes  Peak  to  be  a  single  isolated  rounded  obstacle.     From  a  study  of  large-scale 
topographic  maps,    the  distance    D       between  the  radio  horizons  at  the  top  of  the  peak  is  esti- 
mated to  be    0.040  km.     With    f  =  751  MHz,      d     =   146.0  km,      and    d     =  77.3  km,    we  deter- 


V  =  31.  73  (v  is  positive,    as  both    a       and    P       are  greater  than  zero)  (7.  la) 

o  o 

r  =  0. 6344  (7.  10) 

vp  =  0.858,  (7.9) 

and  therefore  p  =  0.0271 

The  test  described  in  section  7.3   shows  that  the  assumption  of  an  isolated  obstacle 
is  applicable.      The  components  of  basic  transmission  loss  are  then  determined  as  follows: 

Free-space  Loss     L,  ,=   137.0     db  (2.31) 

bf  ■ 

A(v,  0)   =  43.  0     db  (7.  2) 

A(0,  p)  =  6.0     db  (7.4) 

U(v,  p)  =   5.  1     db  (7.  5) 

Totals  are:    A(v,  p)  =   54.  1  db,    and  from  (7.3),      L,  ^  =  L,  ,  +  A(v,  p)   =   191.1  db. 

bd  bf 

The  average  atmospheric  absorption  term,     A    ,    from  figure  3.6  is  0.7  db.      The 

total  basic  transmission  loss  value     L        =   191.8  db,    which  is  equal  to  the  long-term  refer- 

dr 

ence  value     L       .      This  reference  value,    is  strictly  applicable  only  to  those  hours  of  the  year 
cr 

which  are  characterized  by  a  surface  refractivity  of  approximately  2  50  N-units.      The  ex- 
pected behavior  of  the  hourly  median  basic  transmission  loss  for  all  hours  of  the  year  over 
this  path  can  now  be  determined  in  the  form  of  a  cumulative  distribution  by  the  methods  given 
in  section  10 . 

As  this  is  a  knife-edge  diffraction  path,    it  will  be  necessary  to  calculate  cumulative 

distributions     Y(p,  d   )     separately  for  portions  of  the  path  on  each  side  of  Pikes  Peak  and  to 
e 

combine  the  results  as  described  in  subsection  10.6.     Effective  antenna  heights  are  computed 
as  heights  above  curves  fitted  to  terrain  on  each  side  of  Pike's  Peak  using  (5.  15)   and  (5.  16). 
The  curves  are  extrapolated  to  each  antenna  and  to  Pike's  Peak.      The  effective  heights  are 
then  the  heights  of  the  antennas  and  of  the  Peak  above  these  curves.     From  Beulah  to  Pikes 
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Peak  the  terrain  near  the  Peak  is  excluded  because  it  is  partially  shadowed.      Twenty-one 
evenly  spaced  points,      x.  ,    from    d  =  3.3  km    to    d  =   70  km    were  selected  and  the  cor- 
responding terrain  heights    x.     were  read.     From  (5.  15b)     h  =  2100  m,     x  =  36.  6  km,    and 
m  =  2  5.  5,    and  the  straight  line  fitted  to  terrain  is 

h(x)  =  2100  +  25.  5(x  -  36.  6)   meters 

At  the  Beulah  antenna,      x  =  0     and    h(x)  =   1167  meters,    at  Pikes  Peak    x  -  75.  5  kilometers 
and    h(x)  =  3095  meters.      The  effective  antenna  heights  are  then  738  and  1205  meters  (5.  17). 
Using  (10.  1)  to  (10.3)   the  distances     d       =   -187.  1  km,      d       =  33.2  km,    and  the  effective  dis- 

SO  SI 

tance    d      is  33.2  km. 
e 

Similarly  on  the  Table  Mesa  side  much  of  the  terrain  is  shadowed  by  the  small  peak  at 

about  122  km  and  by  the  elevated  area  at  about  202  km.      The  curve  fit  is  therefore  computed 

for  the  intervening  terrain  with    x     =   122.  5  km    and    x       =  200.  5  km.     Using  21  equidistant 

terrain  heights  between  these  points  (5.  15b)  gives    h  -  202  5  m,      x  =   161.  5  km,    and  m  =  -9.3. 

From  (5.  15a)     h(x  =  75.  5)  =  2827,     h(x  =  223.3)  =   1448  meters.      The  effective  antenna  heights 

are  then  1473  and  218  meters  (5.  17),      d       =   -77.  7,    d       -  33.2  km,    and  the  effective  distance 

so  Sl 

d     =  74.  3  km. 
e 

We  thus  have  two  paths  in  tandem  where  the  effective  distances  are  33.  2  and  74.3  km 
respectively.  Cumulative  distributions  are  obtained  using  figures  10.1,  10.2,  10.3,  and  eq- 
uations ( 10  .  4)   to  ( 10 .  7) .      The  frequency  factors  are    g(  10,  f)  =   1.  33,    g(90,  f)   =   1 .  29. 


Table  Mesa' 

'Side 

B- 

eulah  Side 

V(50,  d  ) 

0.2 

0 

figure  10. 

,  1 

Y    (10,  d  ) 
o              e 

4.7 

1.2 

figure  10. 

2 

Y   (90,  d  ) 
o             e 

-3.1 

-0.6 

figure  10. 

2 

Y(10) 

6.3 

1.6 

(iO.6) 

Y(90) 

-4.0 

-0.  77 

(10.6) 

Using  the  reference  value     L        =   191.8  db    and  the  ratios  given  in  (10.  7)   the  predicted 
"  cr 

cumulative  distributions  for  both  portions  and  for  the  entire  path  are  tabulated  below: 


percent  Table  Mesa  Beulah                                 Entire  Path 

0.01  170.6  186.5  170.1 

0.1  174.4  187.4  173.8 

1  179.0  188.6  178.6 

10  185.3  190.2  184.8 

50  191.6  191.8  191.3 

90  195.6  192.6  195.5 

99  198.9  193.2  199.0 

99.9  201.2  193.7  201.3 

99.99  203.2  194.0  203.2 
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The  cumulative  distribution  of  predicted  basic  transmission  loss  for  the  entire  path 
was  obtained  by  convoluting  the  distributions  for  each  part  of  the  path,    as  described  in  sub- 
section 10.  6.      This  cumulative  distribution  is  shown  graphically  in  figure  III.  46  together  with 
a  distribution  derived  from  measurements  over  this  path,    reflecting  10  56  hours  of  data  ob- 
tained in  I960  and  1962. 

The  confidence  limits  on  figure  III.  46  were  derived  assuming  that 

0-   (p)   :=   16.73  +  0.  12  Y^(p) 

2  2  2 

where  the  variance     cr   (50)  =   12.  73  db       given  in  (V.40)  has  been  increased  by    4  db       to  allow 

for  equipment  and  reading  errors. 
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III.  8.  3     Predicted  Transmission  Loss  for  a  Transhorizon  Path 
Predicted  forward  scatter  and  diffraction  losses  are  computed  for  a  283  km  path  from 
Dallas  to  Austin,    Texas.      The  prediction  is  compared  with  measurements  at  a  frequency  of 
104,  5  MHz.     Figure  III.  47  shows  the  great  circle  profile  of  this  path,    obtained  by  the  methods 
of  section  6,    plotted  on  a  linear  scale  with  allowance  for  the  sea  level  curvature.      The  verti- 
cal scale  in  the  figure  is  much  exaggerated  to  show  some  of  the  detail  of  terrain. 

The  effective  earth's  radius  is  determined  as  described  in  section  4.      From  figure 
4.  1  minimum  monthly  values  of    N       are  read  at  approximate  horizon  locations.     Appropriate 

terrain  heights  are  used  in  (4.3)  to  compute  horizon  values  of    N    ,    whose  average  is     N     = 

s  °  s 

306.     From  (4.4)  the  effective  earth's  radius  is  then  8580  km. 

A  reference  value,      L       ,    of  basic  transmission  loss,    which  corresponds  to  winter 
cr 

afternoon  conditions,    is  computed  as  shown  below.     For  each  parameter  the  appropriate 
equation  or  figure  is  referenced. 

Path:     Dallas  to  Austin,    Texas 

d  =  283.  1  km,    f  =   104.  5  MHz,    N     =  306,    a  =  8580  km 

s 

Transmitter  Receiver  Reference 

280.4  m  243.9  m  figure  III.  47 
135.0  m                                           9,8  m                          (6.  IL)  and  (6.  12) 

219.5  m  274,3  m  figure  HI,  47 
39.6  km                                         8,8  km  figure  III.  47 
-3.845  mr                                 +2,933  mr  (6.15) 
12.777  mr                                    19.296mr  (6.18) 

0,768  mr  3.961  mr  (6.  16) 

130.72  km  103.95  km  (6.20) 

0.057  mr  0.021  mr  figure  6.9 

12,834  mr  19.317  mr  (6.19) 

The  angular  distance     9  =  32.  151  mr  =  0.032151     radians  and  the  product    9d    is  9.  10. 

The  ratio  of    a   /  p     =  s  =  0.664.      Then  the  function    F(9d)     from  figure  9.  1  or  figures  III.  12 
o      o 

and  III.  13  is    F(9d)  =   167.0  db.     From  (9.3)  the  height    h     =2.18  km,    and    r|     =   1.06.     From 
(9.4)   and  figure  9.3  the  parameters     r     =   19.01,    H   (r   )  =  0.  15  db,    r     =   1.38,    and    H   (r   )  = 

X  O        i.  ^  O       ^ 

13,5db,     From  figure  9.4  and  (9.  5)     AH     =0.61  and    H     =7.40db.      The  correction  term 

o  o 

F       defined  by  (9.  7)   is  negligible.     Figure  3.  6  shows  that  in  this  case  the  allowance  for  ab- 
o 

sorption    A     ==^0.01  db. 
a 

The  reference  value     L  of  long-term  median  basic  transmission  loss  due  to  for- 

bsr 

ward  scatter  given  by  (9.  1)   is  then 


Parameter 

^s' 

h 

rs 

h     , 
te 

h 
re 

^t' 

\r 

^Lt, 

\r 

^et' 

e 

er 

oo 

Poo 

'of 

9 

or 

%f 

d 
sr 

Aa  , 
o 

^Po 

-D        O 
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30  log  f  =  60. 573 

20  log  d  =  49.038 

F(ed)  =  167.0     figure  9.  1 

F     =  0  (9.  7) 

o 

H     =7.40        (9.  5) 
o 

A     =0.01    figure  3.  6 


and    L,         =   186  db. 
bsr 


Although  this  is  almost  certainly  a  scatter  path,    the  diffraction  loss  for  a  transhorizon 


path  is  computed  as  an  example.     For  horizontal  polarization  over  average  ground  figures 
8.  1  and  8.  2  sh 
are  computed: 


1  and  8.2  show    K(a  =  8497)  =  0.001     and    b  =  93    .      The  following  additional  parameters 


Reference 

a     =   5808  km, 

a     =  39  51  km, 

a    =9179  km. 

a     = 
r 

5804 

(8,8)  and  (8.9) 

C      =   1.135 

Ol 

C       =   1.291 

02 

C  , =  0.975 
ot 

C 
or 

=   1 

.  136 

(8.13) 

K     =  0,001135 

K^  -  0.001291 

K    =  0.000975 

K     = 
r 

0, 

001136 

(8.13) 

B^  =   1.6059 

B     =   1.60  58 

B    =1.60  62 
t 

B     = 

1. 

60  59 

figure  8.  3 

C^(K^)  =  20.03 

C^(K^)  =  20.03 

C^(K^^^)  =  20. 

03 

figure  8.4 

x^  =  385.85 

x,  =   112.25 

La 

X     =  2452,37 
o 

(8.  12)  and  (8.  13) 

F(x  )  =  -11.0 

n-^^  -  -34.  5 

G(x  )  =   107,14 
o 

figure  8.  6 
(8.4) 

Then  from  (8.  1)  the  attenuation  relative  to  free  space    A  =   132.  61  +  A     =   132.  62  db. 

a 

The  free  space  loss  froin  (2.31)  is    L      =  121.87  db    and  the  long-term  median  basic  trans- 
mission loss  due  to  diffraction    L,   ,  =  254.  5  db.     As  expected  this  is  much  more  than  the 

bd 

predicted  scatter  loss  and  the  long-term  reference  value     L        =  L,         =   186  db. 
^  cr  bsr 

In  order  to  predict  long-term  variability  the  effective  distance    d       is  conlputed  using 

(10.1),    (10.2)  and  (10.3b),     d       =64.0  km,     d       =220.  5  km,    and    d     =286.4  km.     From 

SI  so  e 

figures  10.2,    10.3,    III.  32  and  from  (10.6)  the  following  parameters  are: 


'arameter 

Summer 

Winter 

All  Hours 

Reference 

V(50,d  ) 
e 

5.0 

1.0 

3.0 

figure  in.  32 

Y    (10,  d  ) 
o             e 

7.  75 

7.05 

7.  75 

figure  10.2 

Y   (90.  d  ) 
o            e 

-6.35 

-6.35 

-6.35 

figure  10.  2 

g(io,f) 

1.055 

1.055 

1.055 

figure  10.3 

g(90,f) 

1.0  55 

•     1.055 

1.055 

figure   10.  3 

Y(10) 

8.  18 

7.44 

8.  18 

(10.6) 

Y(90) 

-6.70 

-6.  70 

-6.  70 

(10.6) 
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Using  the  reference  value    L        -   186  db    and  the  ratios  given  in  (10,  7)  the  predicted  cumula- 
tive distributions  for  summer,    winter  and  all  hours  are  tabulated  below: 


p  Summer  Winter  All  Hours 

0.01  153.8  160.2  155.8 

0.1  158.7  164.7  160.7 

1  164.6  170.1  166.6 

10  172.8  177.6  174.8 

50  181.0  185.0  183,0 

90  187.7  191.7  189,7 

99  193.2  197.2  195.2 

99,9  197.1  201.1  199.1 

99,99  200,4  204.4  202.4 

These  cumulative  distributions  are  shown  graphically  on  figure  III.  48  together  with 
distributions  derived  from  measurements  over  the  path.      The  data  represent  more  than 
23,  000  hourly  medians  from  measurements  over  a  period  of  more  than  3  years.      The  ar- 
rows on  the  0.  1  percent  values  indicate  that  losses  were  less  than  the  values  plotted  while 
those  on  the  99  and  99.9  percent  values  indicate  that  the  losses  were  greater  than  the  values 
plotted. 

An  example  may  be  included  here  of  the  method  of  mixing  distributions  described  in 
subsection  III.  7.  2.      The  summer  and  winter  distributions  of  data  may  be  mixed  in  order  to 
obtain  the  distribution  of  hourly  median  values  for  all  hours.     Several  levels  of  transmission 
loss  from  165  to  195  db  are  chosen  and  the  value    p    is  read  from  figure  111.48  for  each  dis- 
tribution,    A  weighted  average  is  then  obtained  at  each  level  to  provide  the  mixed  distribu- 
tion of  data  corresponding  to  all  hours.      The  weights  are  the  number  of  hours  of  data  in  each 

case  and  the  weighted  average  is    — — j—    (p  x  12,  160  +  p  x  11,  134), 

'  winter  summer 
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Level  p    .  p  Ave  p  (all  hours) 

winter  summer 

165  <0.05  1.70  0.84 

170  2.1  10.5  6.1 

175  8.3  29.2  18.3 

180  24.2  59.0  40.8 

18  5  50.0  88.3  68.3 

190  86,0  98.7  92.1 

195  >99.75  >99.954  >99.85 

The  weighted  averages  of    p    are  plotted  on  figure  III.  48  at  each  of  these  levels  of 
transmission  loss  and  show  very  good  agreement  with  the  distribution  obtained  directly 
from  the  data. 
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A  LINE-OF-SIGHT  PATH   NEAR  BOULDER,  COLORADO 
f  =  IOO  AND  300  MHz  ,   htg  =  4m  ,  h^g  =9m 
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OBSERVED    INPUT  VOLTAGE   VARIATION    AT  GROUND  STATION   RECEIVER 
FROM   AN   AIRCRAFT  AT  10,000  FEET   TRANSMITTING  ON  328.2  MC/S 

TRANSMITTER   POWER-.   6  WATTS;   TRANSMITTING  AND  RECEIVING  ANTENNA  GAIN:  2.15  db  (RELATIVE  TO  AN  ISOTROPIC) 
GROUND  ANTENNA  HEIGHT:  75  FEET,  TRANSMISSION  OVER  WATER;  6  db  COMMUNICATIONS  SYSTEM  LOSS  ASSUMED  FOR  THEORETICAL  CURVES 
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Annex    IV 
FORWARD  SCATTER 

This  annex  discusses  some  of  the  similarities  and  differences  between  forward  scatter 
from  refractive  index  turbulence  and  forward  scatter  or  incoherent  reflections  from  tropo- 
spheric  layers. 

To  scatter  is  to  spread  at  random  over  a  surface  or  through  a  space  or  substance. 
Scattering  which  tends  to  be  coherent  is  miore  properly  called  forward  scatter,    reflection, 
refraction,    focusing,    diffraction,    or  all  of  these,    depending  on  the  circumstances.     Modes  of 
scattering  as  well  as  mechanisms  of  propagation  bear  these  names.     For  example,    we  may 
speak  of  the  reflection,    refraction,    diffraction,    focusing,    scattering,    and  absorption  of  a  radio 
wave  by  a  single  spherical  hailstone,    and  all  of  these  modes  can  be  identified  in  the  formal 
solutions  of  Maxwell's  equations  for  this  problem. 

The  large  voluine  of  beyond-the-horizon  radio  transmission  loss  data  available  in  the 
frequency  range  40  to  4000  MHz  and  corresponding  to  scattering  angles  between  one  and  three 
degrees  indicates  that  the  ratio     10  corresponding  to  the  transmission  loss.     A,    rela- 

tive to  free  space  is  approximately  proportional  to  the  wavelength,    \,    or  inversely  propor- 
tional to  the  radio  frequency,   f,      [Norton,    I960],    so  that  the  ratio     10  corresponding 

3  -3 

to  the  forward  scatter  basic  transmission  loss  is  approxinnately  proportional  to    \       or  to    f 

This  circumstance  is  more  readily  explained  in  terms  of  forward  scatter  from  layers  [  Friis, 
Crawford,    and  Hogg,    19  57]    or  in  terms  of  glancing  or  glinting  from  brilliant  points  on  ran- 
domly disposed  "feuillets",    [duCastel,    Misme,    Spizzichino,    and  Voge,    1962],    than  in  terms 
of  forward  scatter  from  the  type  of  turbulence  characterized  by  the  modern  Obukhov- 
Kolmogorov  theory  [  Obukhov,    1941,    1953;  Batchelor,    1947,     1953].     There  is  recent  evidence 
[Norton  and  Barrows,    1964]    that  the  wavenumber  spectrum  of  refractivity  turbulence  in  a 
vertical  direction  has  the  same  form  as  the  more  adequately  studied  spectrum  of  variations 
in  space  in  a  horizontal  direction.     Some  mechanism  other  than  scatter  from  refractivity 
turbulence  must  be  dominant  most  of  the  time  to  explain  the  observed  transmission  loss 
values  over  a  majority  of  the  transhorizon  tropospheric  paths  for  which  data  are  available. 
Scattering  from  refractivity  turbulence  and  scattering  from  sharp  gradients  are  mechanisms 
which  coexist  at  all  times  in  any  large  scattering  volume.     Sharp  gradients  always  exist  some- 
where,   and  the  atmosphere  between  them  is  always  somewhat  turbulent.      Power  scattered 
by  these  mechanisms  is  occasionally  supplemented  by  diffraction,    specular  reflection  from 
strong  extended  layers,    and/ or  ducting. 

A  tropospheric  duct  exists,    either  ground-based  or  elevated,    if  a  substantial  amount 
of  energy  is  focused  toward  or  defocused  away  from  a  receiver  as  super-refractive  gradients 
of    N    exceed  a  critical  value  called  a  "ducting  gradient."     This  gradient  is  about  -157  N-units 
per  kilometer  at  sea  level  for  horizontally  launched  radio  waves.      The  duct  thickness  must 
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exceed  about  5062  f    ^  meters  with    f    in    MHz  [Kerr,    1964]    for  a  duct  to  completely  trap 
such  radio  waves.     A  few  useful  references  in  this  connection  are  cited  at  the  end  of  section 
4. 

A  more  or  less  horizontally  homogeneous  "kink"  in  a  refractive  index  profile  may  in- 
dicate the  possibility  of  ducting  for  very  short  wavelengths,    the  presence  of  a  refracting 
layer  for  some  of  the  longer  waves,    and  merely  a  slight  and  random  perturbation  of  average 
atmospheric  conditions  for   other  frequencies,    antenna  locations,    or  antenna  beam  patterns 
and  elevation  angles.      The  layer  that  presents  a  sharp  discontinuity  for  radio  frequencies 
from  30  to  100  MHz  (X.  =   10  to  3  meters)   may  represent  a  relatively  gradual  change  of  re- 
fractive index  at  300  MHz  (\  =   1  meter)   and  higher  frequencies.     A  tropospheric  layer  or 
"feuillet"  requires  a  sufficiently  abrupt  change  in  refractive  index,    usually  associated  with 
fine  weather  conditions,    to  reflect  a  substantial  amount  of  radio  energy  at  the  grazing  angles 
and  frequencies  of  interest.      These  may  be  horizontal  changes,    in  thermals,    for  instance, 
as  v^'ell  as  changes  of  refractivity,      N,    with  height. 

Almost  specular  reflection  from  tropospheric  layers  is  often  observed  between  30 
MHz  and  200  MHz.     At  higher  frequencies,    where  focusing,    defocusing,    and  ducting  are 
common,    and  where  extensive  layers  are  not  sufficiently  abrupt  or  sufficiently  numerous  to 
provide  strong  reflections,    a  number  of  small  and  randomly  oriented  surfaces  come  into 
play.     A  recent  summary  of  the  role  of  the  layer  structure  of  the  troposphere  in  explaining 
tropospheric  propagation  [  Saxton,    Lane,    Meadows,    and  Matthews,    1964]    includes  an  exten- 
sive list  of  references.     Also  useful  are  general  discussions  of  tropospheric  propagation  by 
Bullington  [  1955]  ,    du  Castel  [  I960] ,    Crawford,    Hogg,    and  Kummer  [  1959]  ,    Fengler  [  1964] , 
Fengler,    Jeske,    and  Stilke  [  1964]  ,    Kirby,    Rice,    and  Maloney  [  196  1] ,    Johnson  [  19  58]  ,    Rice 
and  Herbstreit  [  1964],    Shkarofsky  [  1958]  ,    and  Vvedenskii  and  Arenberg  [  1957]  . 

There  are  at  least  three  distinguishing  features  in  most  theories  of  forward  scatter 
from  clouds,    precipitation,    refractive  index  turbulence,    layers,    or  feuillets.     A  calculation 
is  first  made  of  the  expected  or  average  forward  scattering  pattern,    reradiation  pattern,    or 
diffraction  pattern  of  a  scatterer  or  a  group  of  scatterers,    usually  located  in  free  space,    and 
usually  assuming  an  incident  plane  wave  and  a  distant  receiver.     Second,    a  decision  is  made 
that  the  relative  phases  of  waves  scattered  from  individual  raindrops  or  subvolumes  of  re- 
fractivity turbulence  or  feuillets  are  random,    so  that  we  may  simply  add  the  power  contri- 
butions from  these  elements  and  ignore  the  phases.      This  is  an  essential  feature  of  a  random 
scatter  theory.     And  third,    some  way  is  found  to  relate  the  actual  terrain,    atmosphere,    and 
antenna  parameters  to  the  theoretical  model  so  that  a  comparison  may  be  made  between  data 
and  theory. 

The  mechanisms  of  scattering  from  refractivity  turbulence,    reflection  from  elevated 
layers,    and  ducting  are  much  more  sensitive  to  vertical  refractivity  gradients  than  to  the 
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horizontal  gradients  commonly  observed.     The  forward  scatter  theory  used  to  develop  the 
prediction  methods  of  section  9  assumes  that  only  vertical  scales  of  turbulence  or  layer 
thicknesses  are  important.      The  radio  wave  scattered  forward  by  all  the  scattering  sub- 
volumes  visible  to  both  antennas  or  by  all  the  layers  of  feuillets  visible  to  both  antennas  is 
most  affected  by  a  particular  range  of  "eddy  sizes",     IL  ,    or  by  layers  of  an  average  thick- 
ness    SL I  Z  .     A  stack  of  eddies  of  size    i.     must  satisfy  the  Bragg  condition  that  reradiation  by 
adjacent  eddies  shall  add  in  phase.     Reflections  from  the  exterior  and  interior  boundaries 
of  a  layer  will  add  in  phase  if  the  ray  traversing  the  interior  of  the  layer  is  an  odd  number 
of  wavelengths  longer  than  the  ray  reflected  from  the  exterior  boundary.     Either  the  mech- 
anism of  forward  scatter  from  refractivity  turbulence  or  the  mechanism  of  reflection  from 
layers  or  feuillets  selects  a  wavenumber  direction    k     that  satisfies  the  specular  reflection 
condition  corresponding  to  Snell's  law  that  angles  of  incidence  and  reflection,    i|j  ,    are  equal. 
Mathematically,    these  conditions  are  represented  by  the  following  relations: 


2  sin (9/2)    ~    9     ' 


R 
o 

+  R 

1  ^ 
'     o 

^     1 

+  R| 

(IV.  1) 


where    R       and    R    are  unit  vectors  from  the  centers  of  radiation  of  the  transmitting  and  re- 
ceiving antenna,    respectively,    towards  an  elementary  scattering  volume,    or  towards  the 
point  of  geometrical  reflection  from  a  layer.      The  angle  between    R    and    R       is  the  scat- 
tering angle     9    illustrated  in  figure  IV- 1  and  is  thus  twice  the  grazing  angle    ^    for  reflec- 
tion from  a  layer: 


9  =  2t|j=  cos        (-R  •    R   )       radians 
o 


(IV.  2) 


The  plane  wave  Fresnel  reflection  coefficient    q       for  an  infinitely  extended  plane 
boundary  between  homogeneous  media  with  refractive  indices    n     and    n       and  for  horizon- 
tal polarization  [Wait,    1962]   is 


2  2 

sia.^  -  I   2(n^  -  n^)  +  (n,  -  n^)     +  sin    i); 


sinip  + 


2  Z      " 

2(nj  -  n^)  +  (n^  -  n^)      +sin    ^ 


(IV.  3) 


2  2, 

The  following  approximation,    valid  for     (n     -  n  )      <    sin  4J<<    1     is  also  good  for  vertical 

polarization: 


%- 


"2-^1 


2i)j 


exp 


(n^  -n^)^/(24;^) 


(IV.  4) 
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A  differential  amplitude  reflection  coefficient    dq    for  a  tropospheric  layer  is  next 

defined  as  proportional  to  the  difference  between  two  gradients  of  refractive  index,      m    and 

m    ,    where    m    is  the  average  refractive  index  gradient    dn/dz    across  the  layer,    and    m 
o  o 

is  the  average  refractive  index  gradient  for  the  region  in  which  the  layer  is  embedded.     Let 
the  layer  extend  in  depth  from    z  =  0  to  z  =  z       in  the  wavenumber  direction    k.     defined  by 
(IV.  1),    and  write  the  differential  reflection  coefficient  as 


dq  =  dz  (m  -  m  )  /  (Z^ 
o 


(IV.  5) 


Aphasor     exp  [ -i  z  (4tt'Jj  /X)]     is  associated  with    dq  ,    and  the  power  reflection  coefficient 

2 
q       for  a  tropospheric  layer  of  thickness     z       is  approximated  as 


\  dq  exp  [-i  z(4tt4'/X)  ] 


(4-n-)^  \^  4j'     M 


(IV.  6a) 
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M  =  (m-m)        [1  -  cos  (4Tr4JZ    /\)]/(4Tr) 


(IV.  6b) 


If    M    is  assumed  continuous  at    z  =  0     and    z  =  z    ,    somewhat  smaller  values  of    q       and    m 

o 

will  result  [Wait,    1962]. 

Friis,    Crawford,    and  Hogg  [  19  57]   point  out  that  the  power  received  by  reflection 
from  a  finite  layer  can  be  approximated  as  the  diffracted  power  through  an  absorbing  screen 
with  the  dimensions  of  the  layer  projection  normal  to  the  direction  of  propagation.      They 
then  consider  layers  of  large,    small,    and  medium  size  compared  to 


2w  =  2(\R    R/d)      ,       d  ^   R     +  R 
o  o 


(IV.  7) 


which  is  the  width  of  a  first  Fresnel  zone.     Let    b    represent  the  dimensions  of  a  layer  or 
feuillet  in  any  direction  perpendicular  to    k   .     Since     k     is  usually  nearly  vertical,      b    is 
usually  a  horizontal  dinaension.     Adopting  a  notation  which  conforms  to  that  used  elsewhere 

in  this  report,    the  available  power     p       at  a  receiver  at  a  distance    d    from  a  transmitting 

a 

antenna  radiating    p      watts  is 

2    2 
p     __     ^^t^t^r^   ^  [cV)+sV)]    [cV)+sV)] 

^  (4t:  d) 


(IV.  8) 


in  terms  of  Fresnel  integrals  given  by  (III.  33),    where 

u  =  b»7~2~/w  ,         v  =  b  i|ji^~2'/w 


(IV.  9) 
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and     g      and     g       are  antenna  directive  gains.     For  large     u     and     v, 

C^{u)  =  S^u)  =  C^(v)  =  S^(v)  =     1/4   , 
and  for  small     u     and     v,     C    (u)  =  u   ,      C    (v)  =  v    ,      and     S   (u)  =  S^(v)  =  0  . 


For  large  layers,   where     b  >  >  2w: 

4     -6    -2 
P^  =  P^g^g^^     4^      d       M  (IV,10a) 

For  intermediate  layers,   where     b  =  2w: 

P^  =  Ptg^gj,>^    4^"     (RR^d)"    b     M  (IV.  10b) 

For  small  layers,   where     b  <  <  2w; 

2-4  -2     4 

P^  =  P^g^gj,^    4^       ^^^o^       ^     ^  (IV.  10c) 


Forward  scatter  from  layers  depends  on  the  statistics  of  sharp  refractive  index 
gradients  in  the  directions     k     defined  by  (IV.  1).      The  determination  of  these  statistics  from 

radio  and  meteorological  measurements  is  only  gradually  becoming  practical.     A  study  of 

2  4 

likely  statistical  averages  of  the  meteorological  parameters     M,      b   M,      and    b    M    indicates 

that  these  expected  values  should  depend  only  slightly  on  the  wavelength    \    and  the  grazing 

etngle    4",    as  was  assumed  by  Friis,    et  al  [  1957]  .     The  expected  value  of 

[  1   -  cos  (4Tr  4^  z    /\)3 

can  vary  only  between    0     and    2    and  is  not  likely  to  be  either    0     or     2    for  any  reasonable 
assumptions  about  the  statistics  of    z     . 

Available  long-term  median  radio  transmission  loss  data  usually  show  the  frequency 
law  given  by  (IV.  10b)  for  medium-size  layers.      Long-term  cumulative  distributions  of  short- 
term  available  power  ratios  on  spaced  frequencies  rarely  show  a  wavelength  law  outside  the 

2  4 

range  from    X.       to    X.     [Crawford,    Hogg,    and  Kummer,    1959;  Norton  I960].     An  un- 
reported analysis  of  8978  hours  of  matched  simultaneous  recordings  at  159.5,    599,    and 
2120  MHz  over  a  3 10 -km  path  in  Japan  shows  that  this  wavelength  exponent  for  transmission 
loss     p    /p      is  within  the  range  from  2  to  4  ninety-eight  percent  of  the  time.      This  cor- 
responds to  a  wavelength  exponent  range  from    0     to    2     or  a  frequency  exponent  range  from 

2  4 

0     to  -2    for  attenuation  relative  to  free  space  values,    and  to  corresponding  ranges     \       to    X. 

-2  -4 

or    f         to    f         for  values  of  basic  transmission  loss,      L     . 

b 
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Figures    IV.  1(a)  and  IV.  1(b)   illustrate  forward  scattering  from  a  single  small  layer 
and  from  refractivity  turbulence  in  a  single  small  scattering  subvolume  of  the  volume  V  of 
space  visible  to  two  antennas.     Figures  IV.  1(c)   and  IV.  1(d)  illustrate  models  for  the  addition 
of  power  contributions  from  large  parallel  layers,    and  from  scattering  or  reflection  sub- 
volumes,    respectively.     Contributions  from  diffraction  or  ducting  are  ignored,    as  well  as 
returns  from  well-developed  layers  for  which  a  geometrical  reflection  point  is  not  visible  to 
both  antennas.     Combinations  of  these  mechanisms,    though  sonnetimes  important,    are  also 
not  considered  here. 

For  each  of  the  cases  shown  in  figure  IV.  1,    coherently  scattered  or  reflected  power 

p    .     from  the  neighborhood  of  a  point    R    .     is  conveniently  associated  with  a  scattering  sub- 
ai  _  oi 

volume     d   R      =  dv  =  v  (R    .),    so  that  the  total  available  forward  scattered  power  at  a  receiver 
o  1      oi 

is 

N  N 

V  v 


Pa 


=     >      p       =     >      V.  p    .    ^      \      d^R     p   (R    ,  R)    watts  (IV.  11) 

^     -^ai        l_i       X  "^vi  J  o  ^^v'    o 


where 


1=1  1=1  V 


p    .  =  p    ./v.  =  p   (R    ,  R) 

VI         ai      1         V      o 


(IV.  12) 


is  the  available  power  per  unit  scattering  volume  for  the    i         scattering  subvolume,    feuillet. 


Each  of  the  power  contributions     p    .     is  governed  by  the  bistatic  radar  equation.     Omit- 

ai 


or  layer,    and  it  is  assunned  that  only    N       such  contributions  to    p       are  important 

Each  of  the  power  contributions     p    .     is  gover 

ai 

ting  the  subscript    i,    this  equation  may  be  written  as 

^^^  (7^)  (t^' (^^  """■ 

4ttR  4tt  R 

o 

where     a     c       is  the  effective  scattering  cross -section  of  a  single  scatterer  or  group  of 
s      p 

scatterers,    including  the  polarization  efficiency    c       of  the  power  transfer  from  transmitter 
to  receiver.      The  first  set  of  parentheses  in  (IV.  13)   represents  the  field  strength  in  watts 
per  square  kilometer  at  the  point    R    ,    the  second  factor  enclosed  in  parentheses  shows  what 
fraction  of  this  field  strength  is  available  at  the  receiver,    and    )^   g    /(4tt)     is  the  absorbing 
area  of  the  receiving  antenna. 

The  key  to  an  understanding  of  scattering  from  spacecraft,    aircraft,    rain,    hail,    snow, 
refractivity  turbulence,    or  inhomogeneities   such  as  layers  or  feuillets  is  the  scattering  cross- 
section    a     c       defined  by  (IV.  13)   or  the  corresponding  scattering  cross -section  per  unit 
s      p 

volume    a    ,    defined  from  (IV.  12)   and  (IV.  13)   as 

V 
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a^  =  (47r)^  (R^  R)^  p^/(p^  g^  g^  \^)        per  km  {IV.  14) 

This  quantity  is  usually  estimated  by  isolating  a  small  volume  of  scatterers  in  free 
space  at  large  vector  distances    R       and    R,    respectively,    from  the  transmitter  and  receiver. 
If  both  antennas  are  at  the  same  place,    (IV.  13)   becomes  the  monostatic  radar  equation,    cor- 
responding to  backscatter  instead  of  to  forward  scatter. 

The  scattering  cross-sections  per  unit  volume  for  large,    medium,    and  small  layers, 
assuming  a  density  of    N       layers  per  unit  volume,    may  be  obtained  by  substituting  (IV.  10a)   to 
(IV.  10  c)  in  (IV,  14): 

For  large  layers,   where     b    >>  2w: 

a       =  w"^  i|j"^  M  N,     =  \^  4^"^  (R    R/d)^  M  N,  (IV.  15) 

For  intermediate  layers,   where     b=  2w: 

a       =  w^  ijj"^  b^  M  N„     =  \  4;"^(R    R/d)b^  M  N.  (IV. 16) 

V2  £  O  £ 

For  small  layers,   where      b  <  <  2  w  ; 

a       =  ijj"'^  b^  M  N„     =  \^  ii~     b^  M  N,  (IV.  17) 

V3  i  i 

The  modern  Obukhov-Kolmogor ov  theory  of  homogeneous  turbulence  in  a  horizontal 
direction,    when  extended  to  apply  to  the  wavenumber  spectrum  of  instantaneous  variations 
of  refractive  index  in  a  vertical  direction,    predicts  a     )^    ^     or    f ^    law  for  the  variation  with 

wavelength    X    or  carrier  frequency    f    of  either     a       or  attenuation  relative  to  free  space,    or 

5/3  -5/3  ^ 

a    \  or    f  law  for  variations  of  the  transmission  loss     p    /p    .      Theoretical  studies  of 

multiple  scattering  by  Beckmann  [  1961a],    Bugnolo  [  1958],    Vysokovskii  [  1957,    1958],    and  others 
suggest  that  single  scattering  adequately  explains  observed  phenomena.     Descriptions  of  at- 
mospheric turbulence  are  given  by  Batchelor  [  1947,    1953],    de  Jager  [  1952],    Heisenberg 
[1948],    Kolmogoroff  [  1941],    Mer kulov  [  1957] ,    Norton  [  I960] ,    Obukhov  [  1941,     1953],    Rice 
and  Herbstreit  [  1964],    Sutton  [  1955]  ,    Taylor  [  1922] ,    and  Wheelon  [  1957,    1959]. 

The  observed  wavelength  exponent  for  the  Japanese  transmission  loss  data  previously 
noted  was  below  5/3  less  than  two  tenths  of  one  percent  of  the  time,    and  an  examination  of 
other  data  also  leads  to  the  conclusion  that  forward  scatter  from  Obukhov-Kolmogorov  turbu- 
lence can  rarely  explain  what  is  observed  with  frequencies  from  40  to  4000  MHz  and  scat- 
tering angles  from  one  to  three  degrees. 
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Early  recognition  of  this  fact  by  Norton,   Rice,   and  Vogler  [  1955  ]   led  to  the  proposal 
of  a  mathematical  form  for  the  vertical  wavenumber  spectrum  which  would  achieve  agreement 
between  radio  data  and  the  theory  of  forward  scatter  from  refractivity  turbulence  [Norton, 
i960  ].     Radio  data  were  used  to  determine  the  following  empirical  form  for     a    ,   upon  which 
the  predictions  of  section  9  are  based: 

a     =  \  qj"     M  {IV.  18) 

M  =  3  <  (An)^  >/(32  jf^)  (IV.  19) 

where 

An  =  n  -    <An>  (IV.  20) 

is  the  deviation  of  refractive  index  from  its  expected  value     <An>,     and     i       is  a    "scale  of 
turbulence"     [  Rice  and  Herbstreit,    1964]. 

Values  of  the  variance     <(An)    >     of  refractivity  fluctuations  and  scales  of  turbulence 
i       obtained  from  meteorological  data  lead  to  good  agreement  between    (IV.  18)    and  radio 
data  when  an  exponential  dependence  of     M       on  height  is  assumed,    substituting  the  corre- 
sponding value  of     p       in    (IV.  11).     It  is  not  yet  clear  how  the  estimates  of     m,      m   ,     z    ,     b, 
and     N      required  by  the  theory  of  forward  scatter  from  layers  of  a  given  type  can  be 
obtained  from  direct  meteorological  measurements,   nor  how  these  parameters  will  vary 
throughout  the  large  voliime  of  space  visible  to  both  antennas  over  a  long  scatter  path.     It 
does  seem  clear  that  this  needs  to  be  done. 

Data  from  elevated  narrow-beam  antennas  that  avoid  some  of  the  complex  phenomena 
due  to  reflection  and  diffraction  by  terrain,    and  which  select  small  scattering  volumes,    sug- 
gest that  for  scattering  angles  exceeding  ten  degrees,    reflections  from  large  layers  can 
hardly  be  dominant  over  reflection  from  intermediate  and  small  layers  or  from  refractivity 
turbulence.      Preliminary  results  indicate  that  field  strengths  decrease  more  slowly  at  a 

fixed  distance  and  with  scattering  angles    9     increasing  up  to  fifteen  degrees  than  would  be  pos- 

-6 
sible  with  the     9         dependence  of    a       given  by  (IV.  15)   added  to  a  probable  exponential  decay 

with  height  of  the  expected  value  of  the  meteorological  parameter     MN       for  large  layers. 

The  wavelength  and  angle  dependence  of  forward  scatter  characterized  by  the  Obukhov- 

Kolmogorov  turbulence  theory  is  nearly  the  same  as  that  for  small  layers,    given  by  (IV.  17). 

For  scattering  from  refractivity  turbulence: 

a       =X-^/^    ^-''''    M  (IV.  21) 

vo  o 

r(ll/6)  <  (An)^> 
°        4(2.)^^/-'r(l/3)i^^/^ 
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Although  most  of  the  propagation  paths  which  have  been  studied  rarely  show  this  frequency 
dependence,    some  occasionally  do  agree  with    (IV. 21).     In  general,    the  radiowave  scattering 
cross -section  per  unit  vol\ime     a        is  a  weighted  average  of  scattering  from  all  kinds  of 
layers  or  feuillets  and  the  turbulence  between  them. 
Summarizing  the  argiiment: 

-5 
a=a       +a       +a       +a       ^   \>\i       M  (IV.  23) 

V  VO  VI  V2  V3 

-4  -2 

for     10        <   \  <  10        km,      0.01  <  4^  <   0.03     radians,    where      M     has  been  determined  from 

radio  data,    subject  to  the  assumption  that     M     decreases  exponentially  with  height  above  the 
earth's  surface.     Equation    (IV. 23)    is  intended  to  indicate  the  present  state  of  the  twin  arts  of 
formulating  theories  of  tropospheric  forward  scatter  and  comparing  these  theories  with  avail- 
able long-term  median  transmission  loss  data.     A  great  deal  of  available  data  is  not  forward 
scatter  data,    and  it  is  for  this  reason  that  estimates  of  long-term  variability  as  given  in 
section  10  and  annex  III  are  almost  entirely  empirical. 

Also,    for  this  reason,    estimates  of     L  as  given  in  section  9  are  restricted  to  long- 

gP 
term  median  forward  scatter  transmission  loss.     Available  measurements  of  differences  in 

path  antenna  gain  agree  within  the  limits  of  experimental  error  with  the  values  predicted  by 

the  method  of  section  9  whenever  the  dominant  propagation  mechanism  is  forward  scatter. 
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GEOMETRY  FOR  FORWARD  SCATTER 


(a) 


(b) 


(c) 


(d) 


Figure     Iff.l 
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Annex  V 
PHASE  INTERFERENCE  FADING  AND  SERVICE  PROBABILITY 

As  a  general  rule,    adequate  service  over  a  radio  path  requires  protection  against  noise 
when  propagation  conditions  are  poor,    and  requires  protection  against  interference  from 
cochannel  or  adjacent  channel  signals  when  propagation  conditions  are  good.      Optimum  use 
of  the  radio  spectrum  requires  systems  so  designed  that  the  reception  of  wanted  signals  is 
protected  to  the  greatest  degree  practicable  from  interference  by  unwanted  radio  signals  and 
by  noise . 

The  short-term  fading  of  the  instantaneous  received  power  within  periods  of  time 
ranging  from  a  few  minutes  up  to  one  hour  or  more  is  largely  associated  with  random  fluctu- 
ations in  the  relative  phasing  between  component  waves.      These  arrive  at  the  receiving  anten- 
na after  propagation  via  a  multiplicity  of  propagation  paths  having  electrical  lengths  that  vary 
from  second  to  second  and  from  minute  to  minute  over  a  range  of  a  few  wavelengths.     A 
small  part  of  this  short-term  fading  and  usually  all  of  the  long-term  variations  arise  from 
minute -to -minute  changes  in  the  root-sum-square  value  of  the  amplitudes  of  the  component 
waves,     i.  e.  ,    in  short-term  changes  in  the  mean  power  available  from  the  receiving  antenna. 
In  the  analysis  of  short-term  fading,    it  is  convenient  to  separate  the  effects  of  these  phase 
and  root-sum-square  amplitude  changes  on  the  distribution  of  the  instantaneous  received 
power  and,    for  a  fixed  transmitter  power,    on  the  distribution  of  the  instantaneous  trans- 
mission loss.      Multipath  or     "phase  interference  fading"     among  simultaneously  occurring 
modes  of  propagation  usually  determines  the  statistical  character  of  short-term  variability. 

Over  most  transhorizon  paths,   long-term  variability  is  dominated  by    "power  fading", 
due  to  slow   changes    in  average  atmospheric  refraction,    in  the  intensity  of  refractive  index 
turbulence,    or  in  the  degree  of  atmospheric  stratification.      The  distinction  between  phase 
interference  fading  and  power  fading  is  somewhat  arbitrary,   but  is  nevertheless  extremely 
useful.     Economic  considerations,    as  contrasted  to  the  requirements  for  spectrum  conserva- 
tion,   indicate  that  radio  receiving  systems  should  be  designed  so  that  the  minimum  practicable 
transmitter  power  is  required  for  satisfactory  reception  of  the  wanted  signals  in  the  presence 
of  noise.      Fading  expected  within  an  hour  or  other  convenient     "short"    periods  of  time  is 
allowed  for  by  comparing  the  median  wanted  signal  power     p         available  at  the  receiver  with 
the  median  wanted  signal  power     p  which  is  required  for  satisfactory  reception  in  the 

presence  of  noise.     This  operating  sensitivity     p  assumes  a  specified  type  of  fading 

wanted  signal  and  a  specified  type  of  noise,   but  does  not  allow  for  other  unwanted  signals. 

In  the  presence  of  a  specified  unwanted  signal,   but  in  the  absence  of  other  unwanted 
signals  or  appreciable  noise,   the  fidelity  of  information  delivered  to  a  receiver  output  will 
increase  as  the  ratio     r        of  wanted -to -unwanted  signal  power  increases.     The  degree  of 
fidelity  of  the  received  information  may  be  measured  in  various  ways.     For  example,   voice 
signals  are  often  measured  in  terms  of  their  intelligibility,   television  pictures  by  subjective 
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observation,    and  teletype  signals  by  the  percentage  of  correctly  interpreted  received  charac- 
ters.    A  specified  grade  of  service  provided  by  a  given  wanted  signal  will  guarantee  a 
corresponding  degree  of  fidelity  of  the  information  delivered  to  the  receiver  output.      For 
example,   a  Grade  A  teletype  service  could  be  defined  as  one  providing    99-99    percent 
error-free  characters,   while  a  Grade  B  service  could  be  defined  as  one  providing    99.9    per- 
cent error-free  characters. 

The  protection  ratio     r  required  for  a  given  grade  of  service  will  depend  upon  the 

ur 

nature  of  the  v/anted  and  un^wanted  signals;     i.  e.  ,    their  degree  of  modulation,    their  location 

in  the  spectrum  relative  to  the  principal  and  spurious  response  bands  of  the  receiving  system, 

and  their  phase  interference  fading  characteristics.      The  use  of  receiving  systems  having 

the  smallest  values  of     r  for  the  kinds  of  unwanted  signals  likely  to  be  encountered  will 

ur 

permit  the  same  portions  of  the  spectrum  to  be  used  simultaneously  by  the  maximum  number 

of  users.     For  instance,      FM     with  feedback  achieves  a  reduction  in     r  for  a  cochannel 

ur 

unwanted  signal  by  occupying  a  larger  portion  of  the  spectrum.      But  optimum  use  of  the  spec- 
trum requires  a  careful  balance  between  reductions  in     r  on  the  same  channel  and  on 

ur 

adjacent  channels,    taking  account  of  other  systenn  isolation  factors  such  as  separation  between 

channels,    geographical  separation,   antenna  directivity,    and  cross -polarization. 

Npte  that  the  operating  sensitivity     p  is  a  measure  of  the  required  magnitude  of 

the  median  wanted  signal  power  but     r  involves  only  the  ratio  of  wanted  to  unwanted  signal 

ur 

powers.      For  optimum  use  of  the  spectrum  by  the  maximum  number  of  simultaneous  users, 
the  transmitting  and  receiving  systems  of  the  individual  links  should  be  designed  with  the 
primary  objective  of  ensuring  that  the  various  values  of     r       exceed     r  for  a  large  per- 

centage of  the  time  during  the  intended  periods  of  operation.      Then  sufficiently  high  trans- 
mitter powers  should  be  used  so  that  the  median  wanted  signal  powers     p         exceed     p 

m  mr 

for  a  large  percentage  of  the  time  during  the  intended  period  of  reception  at  each  receiving 
location.     This  approach  to  frequency  assignment  problems  will  be  unrealistic  in  a  few  cases, 
such  as  the  cleared  channels  required  for  radio  astronomy,   but  these  rare  exceptions  merely 
serve  to   test      the  otherwise  general  rule  [Norton,    1950,    1962  and  Norton  and  Fine   1954]    that 
optimum  use  of  the  spectrum  can  be  achieved  only  when  interference  from  other  signals  rather 
than  from  noise  provides  the  ineluctable  limit  to  satisfactory  reception. 

This  annex  discusses  the  requirenaents  for  service  of  a  given  grade  g,  how  to  esti- 
mate the  expected  time  availability  p  of  acceptable  service,  and,  finally,  how  to  calculate 
the  service  probability      Q     for  a  given  time  availability. 

In  this  annex  all  quantities  expressed  in  decibels  are  denoted  by  capital  letters;  for  ex- 
ample,   X  =   10  log  X  ,    where    x    is  usually  a  power  or  a  power  ratio. 
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V.  1     The  Two  Components  of  Fading 

Both  the  wanted  signal  power     p.     and  the  unwanted  signal  power     p        available  to  a 

1  ui 

receiving  system  will  usually  vary  from  minute  to  minute  in  a  random  or  unpredictable 
fashion.     It  is  convenient  to  divide  the  "instantaneous  received  signal  power"  P.  =  10  log  p. 
into  two  or  three  additive  components  where     p.     is  defined  in  the  usual  way  as  the  average 

power  for  a  single  cycle  of  the  radio  frequency,    so  as  to  eliminate  the  variance  of  power 

2 
associated  with  the  time  factor    cos    (wt)  : 

P.  =  P       +  Y.  =  P     (50)  +  Y  +  Y.         dbw  {V.  1) 

1  m         1  m  1 

P         is  that  component  of     P.     which  is  not  affected  by  the  usually  rapid  phase  interference 

m  1 

fading.      P         is  most  often  identified  as  the  short-term  median  of  the  available  po'wer     P.      at 

the  receiving  antenna.      P     (50)     is  the  median  of  all  such  values  of     P     ,    and  is  most  often 
"  m  m 

identified  as  the  long-term  median  of     P..      In  terms  of  the  long-term  median  transmission 

loss   L    (50)    and  the  total  power     P        radiated  from  the  transmitting  antenna: 
m  r 

P     (50)  =  P     -L     (50)       dbw  (V.  2) 

m  r  m 

The  characteristics  of  long-term  fading  and  phase  interference  fading,    respectively,    are 
described  in  terms  of  the  two  fading  components      Y     and     Y.    defined  by    (V.  1): 

Y  =  P       -  P     (50)    ,         Y.  =  P.   -  P  (V.  3) 

m  m  1  1  m 

The  long  term  for  which  the  median  power,      P     (50),    is  defined  may  be  as  short  as 

rn 

one  hour  or  as  long  as  several  years  but  will,    in  general,    consist  of  the  hours  within  speci- 
fied time  blocks  such  as,    for  example,    the  hours  between  6  PM  and  midnight  in  the  summer 
months.     However,    it  will  be  convenient  in  the  case  of  most  continuously  operating  services 
to  consider  that    P     (50)     refers  to  the  median  power  over  a  long  period  of  time,    including 
all  hours  of  the  day  and  all  seasons  of  the  year  for  a  period  of  several  years. 

The  observations  of  long-term  variability  summarized  in  section  10  and  in  annex  III 
show  that     P         is  a  very  nearly  normally  distributed  random  variable  characterized  by  a 
standard  deviation    that     may  range  from  one  decibel  over  a  period  of  time  of  the  order  of 
one  hour  up  to  ten  decibels  for  periods  of  the  order  of  several  years;    these  values  of  stand- 
ard deviation  are  representative  only  of  typical  beyond-the -horizon  propagation  paths  and 
vary  widely  for  other  propagation  conditions. 

For  periods  as  short  as  an  hour,    the  variance  of     Y.  =  P.   -  P         is  generally  greater 

than  the  variance  of     P     .      The  long-term  variability  of     P         is  identified  in  section  10  with 
m  rn 

the  variability  of  hourly  medians,    expressed  in  terms  of 
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Y(p)  =  P^(p)  -  P^(50)  =  L^(50)  -  L^(p)  (V.  4) 


where  P  (p)  is  the  hourly  median  wanted  signal  power  P  exceeded  for  p  percent  of 
all  hours,  and  L  (p)  is  the  corresponding  transmission  loss  exceeded  for  (100 -p)  per- 
cent of  all  hours. 

Often,    data  for  a  path  are  available  in  terms  of  the  long-term  cumulative  distribution 

of  instantaneous  power.      P.;    that  is,    we  know     P.(p)     versus     p,     but  not     P     (p)     versus     p. 

11  m 

Let     q     represent  the  probability  that     P.     will  exceed     P.(q)     for  a  given  value  of     P     ; 
this  is  the  same  as  the  probability  that     Y.     will  exceed     Y.(q).      Then  an  approximation  to 
the  cumulative  distribution  function     L.(p)     versus      P     is  given  by 

L.(P)  =    L^(50)  ±   [Y^(p)  +  Y^^(q)  l''^  (V.  5) 

where      Y.  (q  =  p/100)     is  evaluated  for  the  particular  type  of  phase  interference  fading  charac- 
teristic of  a  wanted  signal  level  exceeded     p     percent  of  the  time.      The  plus  sign  in    {V.  5)    is 
used  for     p  >  50  percent,    and  the  minus  sign  for     p  <   50  percent. 
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V.2     The  Nakagami-Rice  Distribution 

For  studies  of  the  operating  sensitivity    P  of  a  receiver  in  the  presence  of  a 

mr 

rapidly  fading  wanted  signal,    and  for  studies  of  the  median  wanted  signal  to  unwanted  signal 
ratio    R      (g)     required  for  a  grade     g     service,    it  is  helpful  to  consider  a  particular   stat- 
istical model  which  may  be  used  to  describe  phase  interference  fading.      Minoru  Nakagami 
[  1940]    describes  a  model  which  depends  upon  the  addition  of  a  constant  signal  and  a  Rayleigh- 
distributed  random  signal  [  Rayleigh  1880;  Rice,    1945;  Norton,    Vogler,    Mansfield  and  Short, 
1955;  Beckmann,    1961a,     1964].     In  this  model,    the  root-sum-square  value  of  the  amplitudes 
of  the  Rayleigh  components  is    K    decibels  relative  to  the  amplitude  of  the  constant  component. 
K  =  -oo     corresponds  to  a  constant  received  signal.      For  a  Rayleigh  distribution,      K  =  +«> 
and  the  probability    q    that  the  instantaneous  power,      p.  ,    will  exceed    p.(q)     for  a  given  value 
of  the  short-term  median  power,,        p      ,    may  be  expressed: 

P.(q)  log    Z 


q[p.  >  p.(q)  |p^]  =  exp 


1     ,  [K=  +«] 


m  ' 


Alternatively,    the  above  may  be  expressed  in  the  following  forms: 

q[y.   >y.(q)]    =  q[Y.   >  Y.(q)]    =  exp[  -y.(q)  log^Z]     ,        [K=+oo] 


(V.  6b) 


Y.(q)  =   5.21390  +  10  log    {log  (1/q))       ,         [k=+=o] 

Figure  V.  1  and  Table  V.  1   show  the  Nakagami-Rice  phase  interference  fading  dis- 
tribution   Y.(q)     as  a  function  of    K    for  particular  values  of    q  .     It  is  evident  from  figure 
V.  1  that  the  distribution  of  phase  interference  fading  depends  only  on    K  .      The  utility  of  this 
distribution  for  describing  phase  interference  fading  in  ionospheric  propagation  is  discussed 
in  CCIR  report  [  1963k]    and  for  tropospheric  propagation  is  demonstrated  in  papers  by  Norton, 
Rice  and  Vogler  [  1955],    Janes  and  Wells  [  1955],    and  Norton,    Rice,    Janes,    and  Barsis  [  1955] 
Bremmer  [1959]  and  Beckmiann  [l961a]   discuss  a  somewhat  more  general  fading  model. 

For  within-the -horizon  tropospheric  paths,    including  either  short  point-to-point  ter- 
restrial paths  or  paths  from  an  earth  station  to  a  satellite,      K    will  tend  to  have  a  large 
negative  value  throughout  the  day  for  all  seasons  of  the  year.     As  the  length  of  the  terrestrial 
propagation  path  is  increased  or  the  elevation  angle  of  a  satellite  is  decreased  so  that  the 
path  has  less  than  first  Fresnel  zone  clearance,    the  expected  values  of    K    will  increase 
until,    for  some  hours  of  the  day,      K    -will  be  greater  than  zero  and  the  phase  interference 
fading  for   signals  propagated  over  the  path  at  these  times  will  tend  to  be  closely  represented 
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by  a  Rayleigh  distribution.     For  most  beyond-the-horizon  paths    K    will  be  greater  than  zero 

most  of  the  time.     For  knife-edge  diffraction  paths    K    is  often  much  less  than  zero.     When 

signals  arrive  at  the  receiving  antenna  via  ducts  or  elevated  layers,    the  values  of    K    may 

decrease  to  values  much  less  than  zero  even  for  transhorizon  propagation  paths.     For  a 

given  beyond-the-horizon  path,      K    will  tend  to  be  negatively  correlated  with  the  median 

power  level    P      ;i.e.,    large  values  of    K     are  expected  with  small  values  of    P      .      For 

"^  m 

some  within-the-horizon  paths,      K     and    P         tend  to  be  possitively  correlated. 

Discussions  of  fading  rate  and  the  distribution  of  the  duration  of  fading  are  given  in 
CCIR  Report  No.    242  for  tropospheric  propagation  [  1963  k]     and  in  CCIR  Report  No.    266  for 
ionospheric  propagation   [ 1963^  ]. 

It  is  assumed  in  this  report  that  a  particular  value  of     K     may  always  be  associated 
with  any  time  availability,      p.     However,   few  data  analyses  of  this  kind  are  presently  avail- 
able.    A  program  of  data  analysis  is  clearly  desirable  to  provide  empirical  estimates  of     K 
versus     p     for  particular  climates,    seasons,    times  of  day,    lengths  of  recording,    frequencies, 
and  propagation  path  characteristics.      It  should  be  noted  that     K     versus     p     expresses  an 
assumed  functional  relationship.      The  percentage     p     is  the  time  availability  associated  with 
the  cumulative  distribution     L     (p)     versus     p. 

An  ainalysis  of  data  for  a  single  day's  recording  at     1046    MHz    over  a  364-kilometer 

path  is  presented  here  to  illustrate  how  a  relationship  between     K     and     p     may  be  established. 

Figure    V.2    shows  for  a  single  day  the  observed  interdecile  range     P.  (0.1)  -  P.  (0.9)  = 

L.  (0.9)  -  L.  (0.1)     for  each  five-minute  period  plotted  against  the  median  transmission  loss 

L         for  the  five-minute  period.      Figure  V.  1    associates  a  value  of     K     with  each  value  of 
m 

L.  (0.9)  -  L.  (0.1),   and  the  cumulative  distribution  of     L         associates  a  time  availability     p 
11  m  /      JT 

with  each  value  of     L     .     In  figure    V.2,     K     appears  to  decrease  with  increasing     L         for 
m  rn 

the  hours     0000  -  1700,   although  the  usual  tendency  over  long  periods  is  for     K     to  increase 

with     L  The  straight  line  in  the  figure  is  drawn  through  medians  for  the  periods 

0000-1700    and     1700-2400    hours,    with  linear  scales  for     K     and     L     ,     as  shown  by  the 

m 

inset  for  figure  V.2.     The  corresponding  curve  of    K     versus     p     is  compared  with  the  data 
in  the  main  figure. 

Figure  V.3    shows  for     f  =  2  GHz     and  30-meter  antenna  heights  over  a  smooth  earth 
a  possible  estimate  of     K     versus  distance  and  time  availability.      These  crude  and  quite 
speculative  estimates  are  given  here  only  to  provide  an  example  of  how  such  information 
with    (V.  5)   may  be  used  to  obtain  curves  of     L.  (p)     versus     p.      The  solid  curves  in  the 
lower  part    of  figure   V.3    show  how     L.  (p)     varies  with  distance  for     p  =  0.01,      1,      50,      99, 
and    99.99  percent,   where     L.     is  the  transmission  loss  associated  with  the  instantaneous 
power,      P.. 
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V.  3    Noise -Limited  Service 

Let     P       =10  log   p  dbw    represent  the  median  of  that  component  of  available 

m  m 

wanted  signal  power  which  is  associated  with  phase  interference  fading  at  the  terminals  of  an 
equivalent  loss-free  receiving  antenna,    and  let     P        (g),     the    "operating  sensitivity"     ex- 
pressed in  dbw,    represent  the  minimum  value  of     P         which  will  provide  a  grade     g     service 
in  the  presence  of  noise  alone.     This  operating  sensitivity     P  assumes  a  specified  type  of 

fading  wanted  signal  and  a  specified  type  of  noise,   but  does  not  allow  for  other  unwanted 
signals.     A  detailed  discussion  of  the  effective  noise  bandwidth,   the  operating  noise  factor, 
and  the  operating  sensitivity  of  a  receiving  system  will  be  found  in  a  recent  report,     "Optimum 
Use  of  the  Radio  Frequency  Spectrum,"    prepared  under  Resolution  1  of  the  CCIR    [Geneva, 
1963c  ]. 

The  median  value     p  of  the  total  noise  power  in  watts  in  a  bandwidth     b     cycles 

mn 

per  second  at  the  output  load  of  the  linear  portion  of  a  receiving  system  includes  external 
noise  accepted  by  the  antenna  as  well  as  noise  generated  within  the  receiving  system,   includ- 
ing both  principal  and  spurious  responses  of  the  antenna  and  transmission  line  as  well  as  the 
receiver  itself.      This  total  noise  power  delivered  to  the  pre-detection  receiver  output  may  be 
referred  to  the  terminals  of  an  equivalent  loss-free  antenna  (as  if  there  were  only  external 
noise  sources)   by  dividing    p  hy    g    ,    the  maximum  value  of  the  operating  gain  of  the  pre- 

detection  receiving  system. 

The  operating  noise  factor  of  the  pre-detection  receiving  system,      f       ,    may  be  ex- 

op 

pressed  as  the  ratio  of  the  "equivalent  available  noise  power"     p        /g       to  the  Johnson  noise 

mn      o 

power     kT    b    that  would  be  available  in  the  band    b    from  a  resistance  at  a  reference  absolute 

-23 

temperature     T     =  288.37  degrees  Kelvin,    where     k=   1,  28054  X  10  joules  per  degree  is 

Boltzmann's  constant: 

P        /g 

mn      o 

op  k  T     b  ^         ' 

o 

or  in  decibels  : 

F       =    (P        -  G   )    -    (B  -  204)      db  (V.  8) 

op  mn         o 

The  constant    204    in    (V.8)    is      -10  1og(kT   ). 

o 

Note  that  the  available  power  from  the  antenna  as  defined  in  section  2  has  the  desired 
property  of  being  independent  of  the  receiver  input  load  impedance,   making  this  concept 
especially  useful  for  the  definition  and  measurement  of  the  operating  noise  factor     f         as 
defined  under  CCIR  Resolution  1    [Geneva,    1963c]. 

At  frequencies  above     100     MHz,   where  receiver  noise  rather  than  external  noise 

usually  limits  reception,      f         is  essentially  independent  of  external  noise.     In  general,      f 

is  proportional  to  the  total  noise    p  delivered  to  the  pre-detection  receiver  output  and  so 

mn  ^ 

measures  the  degree  to  which  the  entire  system,    including  the  antenna,    is  able  to  discriminate 
against  both  external  noise  and  receiver  noise. 
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Compared  to  the  total  range  of  their  long-term  variability,    it  is  assumed  that    P 

m 

and    P       (g)     are  hourly  median  values;  i.e.,    that  long-term  power  fading  is  negligible  over 
mr  

such  a  short  period  of  time.      Let    G  represent  the  hourly  median  operating  signal  gain  of 

a  pr e-detection  receiving  system,    expressed  in    db    so  that    P      +  G  very  closely  approxi- 

m  ms 

mates  the  hourly  median  value  of  that  component  of  available  wanted  signal  power  delivered  to 

the  pre-detection  receiver  output  and  associated  with  phase  interference  fading.      Let    P 

mn 

represent  the  hourly  median  value  of  the  total  noise  power  delivered  to  the  pre-detection  re- 
ceiver output.  The  median  wanted  signal  to  median  noise  ratio  availableat  the  pre-detection 
receiver  output  is  then 

R       =     [P       +  G         ]    -  P  db  (V.  9a) 

m  m         ras  -"  mn  ^       ^   ' 

and  the  nninimum  value  of    R         which  will  provide  a  desired  grade  of  service  in  the  presence 

m 

of  noise  alone  is 

R       (g)  =    [P       (g)  +  G        ]   -  P  db  (V.9b) 

mr  mr  ms  mn  ^  ' 
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V.4    Interference -Limited  Service 

Only  the  short-term  variations  of  the  wanted  signal  power       p.      and  the  unwanted 

signal  power       p    .      associated  with  phase  interference  fading  are  used  in  deternaining     r      (g), 

the  ratio  of  the  median  wanted  signal  power     p        to  the  median  unwanted  signal  power     p 

m  um 

required  to  provide  a  specified  grade  of  service    g  .     Separate  account  is  taken  of  the  varia- 
tions with  time  of  the  median  power  levels     P       =10  log  p        and    P         =10  log  p         .      This 

m  m  um  um 

separation  of  the  total  fading  into  a  phase  interference  component    Y.     and  the  more  slowly- 
varying  component    Y    appears  to  be  desirable  for   several  reasons:     (1)     those  variations  of 
that  component  of  the  instantaneous  received  power     Y.     which  are  associated  with  phase  in- 
terference alone  may  be  expected  to  occur  completely  independently  for  the  wanted  and  un- 
wanted signals,    and  this  facilitates  making  a  more  precise  determination  of  the  proper  value 

for    R      (g)  ,    (2)   the  random  variable     Y.     follows  the  Nakagami-Rice  distribution,    as  illustrated 
ur  1 

on  figure  V.  1,    while  variations  with  time  of  the  remaining  component    Y     are  approxinnately 

normally  distributed,    (3)   the  variations  with  time  of    P        and  of    P  tend  to  be  correlated 

m  um 

for  inost  wanted  and  unwanted  propagation  paths,    and  an  accurate  allowance  for  this  correla- 
tion is  facilitated  by  separating  the  instantaneous  fading  into  the  two  additive  components     Y 

and    Y    ,    and  (4)  most  of  the  contribution  to  the  variance  of    P        with  time  occurs  at  low  fluc- 

1  m 

tuation  frequencies  ranging  from  one  cycle  per  year  to  about  one  cycle  per  hour,    whereas 

most  of  the  contribution  to  the  variance  of    Y.     occurs  at  the  higher  fluctuation  frequencies 

greater  than  about  one  cycle  per  hour. 

Let    R         denote  the  instantaneous  ratio  between  the  instantaneous  wanted  signal  power 
ui 

P       +  Y      and  the  instantaneous  unwanted  signal  power     P         +  Y      : 
m  1  um  ui 

R      =  P       +  Y.  -  P         -  Y    .  =  R     +  Z.  (V.  10) 

ui  m         1         um         ui  u         1 


where 


Z.    =   Y.  -  Y    .,    and  R     =  P       -  P  (V.  11) 

1  1  ui  u  nn  um 


Note  that  the  cumulative  distribution  function    Y.(q,  K)     for     Y.     will  usually  be  different  from 
the  cumulative  distribution  function    Y.(q,  K   )     for     Y    .     since  the  wanted  signal  propagation 
path  will  differ  from  the  propagation  path  for  the  unwanted  signal.      Let    Z.  >  Z.(q,  K,  K   ) 
with  probability    q;  then  the  approximate  cumulative  distribution  function  of    Z.     is  given  by: 


Z.     (q,  K,  K    )  =   ±  jYf(q.  K)  +  Yf(l-q,  K    )  (V.  12) 

Xo-  U.  ^        X  X  U. 

In  the  above,    the  plus  sign  is  to  be  used  when    q  <  0.  5  and  the  minus  sign  when    q  >  0.  5  ;  note 

that    -Z.     (1-q,  K,  K   )  =  Z.     (q,  K,  K   )  .      This  method  of  approximation  is  suggested  by  two  ob- 
la  u  la  u 

servations:     (1)   service  may  be  limited  for     100  q  percent  of  a  short  period  of  time  either  by 
downfades  of  the  wanted  signal  corresponding  to  a  level  exceeded  with  a  probability    q    or  by 
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upfades  of  an  unwanted  signal  corresponding  to  a  level  exceeded  with  a  probability     1-q,    and 

(2)   the  standard  deviation  of  the  difference  of  two  uncorrelated  random  variables     Y.     and    Y    . 

2 

equals  the  square  root  of  the  sum  of  their  variances,    and  under  fairly  general  conditions    Z.  (q) 

2  2  ^ 

is  very  nearly  equal  to    Y.  (q)  +  Y    .(100-q)  ,     Equation  (V.  11)   is  based  on  the  reasonable  as- 
sumption that    Y      and    Y         are  independent  random  variables  and,    for  this  case,     (V.  11) 
i  ui 

would  be  exact  if    Y      and    Y         were  normally  distributed.      The  departure  from  normality  of 

i  ui 

the  distribution  of    Y.     is  greatest  in  the  limiting  case  of  a  Rayleigh  distribution,    and  for  this 

special  case,    the  following  exact  expression  is  available  [  Siddiqui,    1962]; 

Z.  (q,  «,  CO)  =   10  log        (~"0  (V-13) 

Table  V.2  compares  the  above  exact  expression     Z.  (q,  «,  oo)     with  the  approximate  expres- 
sion    Z.    (q,  «,  ").     Note  that  the  two  expressions  differ  by  less  than    0.2    db  for  any  value  of 

ia 
q     and,    since  this  difference  may  be  expected  to  be  even  smaller  for  finite  values  of     K,   it 

appears  that    (V.  12)    should  be  a  satisfactory  approximation  for  most  applications  and  for 
any  values  of     K     and     q. 

Table  V.  2 

The  Cumulative  Distribution  Function     Z.(q,  «,  «)    for 
the  Special  Case  of  the  Ratio  of  Two  Rayleigh-Distributed  Variables 


q 

Z.(q,  «>,  «) 

Z.  (q,  «>,  «) 
la 

Z  -  Z 
1    la 

db 

db 

db 

0. 0001 

39.999  57 

40. 0178 

-0.01823 

0.0002 

36.98883 

37. 0362 

-0.04737 

0.0005 

33.00813 

33. 0757 

-0. 06757 

0.  001 

29.99566 

30. 1099 

-0. 11424 

0.002 

26.98101 

27. 1216 

-0. 14059 

0.005 

22. 98853 

23. 1584 

-0. 16987 

0.01 

19. 95635 

20. 1420 

-0.  18565 

0.  02 

16.90196 

17.0949 

-0. 19294 

0.05 

12.78754 

12.9719 

-0.  18436 

0.  1 

9.  54243 

9.7016 

-0. 15917 

0.2 

6. 02060 

6.  1331 

-0.  11250 

0.5 

0 

0 

0 
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Let     R        (g)     denote  the  required  value  of     R        for  non-fading  wanted  and  unwanted 
uro  u 


signals  and  it  follows  from    (V.IO)    that  the  instantaneous  ratio  for  fading  signals  will  exceed 

R        (g)     with  a  probj 
uro"^ 

100  q     provided  that: 


R        (g)     with  a  probability  at  least  equal  to     q     or  for  a  percentage  of  time  at  least  equal  to 
uro 


\  >  V  (S'  'I'  ^'  ^u>  =  Vo^g)  -  ^i(^'  K'  ^u)  (^-  ^4) 

The  use  of    (V.  14)    for  determining  an  allowance  for  phase  interference  fading  will  almost 

always  provide  a  larger  allowance  than  will  actually  be  necessary  since     (V.  14)    was  derived 

on  the  assumption  that     R         (g)     is  constant.      For  most  services,     R         (g)     will  not  have  a 

uro  uro   ° 

fixed  value  for  non-fading  signals  and  will  instead  have  either  a  probability  distribution  or  a 

grade  of  service  distribution;    in  such  cases     R       (g)     should  be  determined  for  a  given     q 

by  a  convolution  of  the  distributions  of     R        (g)     and     -  Z..     In  still  other  cases  the  mean 

uro  1 

duration  of  the  fading  below  the  level     Z.(q,  K,  K    )     will  be  comparable  to  the  mean  duration 

of  the  individual  message  elements  and  a  different  allowance  should  then  be  made.     In  some 

cases  it  may  be  practical  to  determine    R  as  a  function  of    g,    q,    K,    and  K       in  the 

ur  u 

laboratory  by  generating  wanted  and  un'wanted  signals  that  vary  with  time  the  same  as     Y. 
and     Y    ..     This  latter  procedure  will  be  successfiil  only  to  the  extent  that  the  fading  signal 
generators  properly  simulate  natural  phase  interference  fading  both  as  regards  their  ampli- 
tude distributions  and  their  fading  duration  distributions.     Since  the  present  report  is  intended 
to  deal  only  with  general  definitions  and  procedures,    functions     R      (g)     applicable  to  parti- 
cular kinds  of  wanted  and  unwanted  signals  and  which  include  an  appropriate  phase  interfer- 
ence fading  allowance  should  be  developed  by  other  CCIR  working  groups  with  members 
having  a  greater  familiarity  with  the  requirements  of  the  wanted  service  and  the  nature  of 
the  phase  interference  fading  of  the  wanted  and  unwanted  signals  likely  to  be  encountered, 
and  the  nature  of  any  diversity  system  which  may  be  used. 

The  ratio    R       will  also  vary  with  time  since  it  has  been  defined  as  an  hourly  median 

equal  to  the  difference  between    P     ( 50)  +  Y    and    P        (50)  +  Y     : 

m  um  u 

R     =  P       -  P         =  P     (50)  -  P       (50)  +  Z  (V.  15) 

u  m         um         m  um  *  ' 

Z  =  Y  -  Y  (V.  16) 

u  *  ' 


The  random  variables     Y     and     Y       tend  to  be  approximately  normally  distributed  with  a 
positive  correlation  coefficient     p      which  will  vary  considerably  with  the  propagation  paths 
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and  the  particular  time  block  involved.      For  the  usual  time  block  involving  all  hours  of  the 

day  for  several  years  preliminary  analyses  of  data  indicate  that    p    will  usually  exceed    0.4 

even  for  propagation  paths  in  opposite  directions  from  the  receiving  point.      Z    will  exceed  a 

value     Z    (p)     for  a  percentage  of  time     p    where  the  approximate  cumulative  distribution  func- 
a 

tion    Z    (p)     of    Z    is  given  by 


J  Y^(p)  +  --^■ 


Z    (p)  =   ±^  Y   (p)  +  Y    (100-p)  +  2pY(p)  Y,(100-p)  (V.  17) 

d  ^  U  U 

In  the  above  the  plus  sign  is  to  be  used  for     p  <    50%     and  the  minus  sign  for     p  >  50%,   while 

Z    (50)  =  0.     It  follows  from    (V.15)    and    {V.17)    that     R       will  exceed     R      (g)     for  at  least  a 
a  u  ur 

percentage  of  time     p     provided  that 

P     (50)  -  P       (50)  >R      (g)  -  Z    (p)  (V.  18) 

m  um  ur  a 


In  some  cases  it  may  be  considered  impractical  to  determine  the  function    R      (g) 

ur 

by  adding  an  appropriate  phase  interference  fading  allowance  to     R        (g)  ;    in  such  cases  it 
may  be  useful  to  use  the  following  approximate  relation  which  will  ensure  that  the  instantan- 
eous ratio     R    .  >  R        (g)     for  at  least  a  percentage  of  time     p: 
ui  uro 


4-'-  -'■ 


P     (50)  -  P       (50)  >R        (g)    ±a/  Z    (p)  +  Z.  (0.01  p,  K,  K   )  (V.  19) 

m  um  uro  »      a  i  u 


In  the  above,    the  minus  sign  is  to  be  used  for     p  <   50%     and  the  plus  sign  for     p  >  50%. 
Although    (V.19),    or  its  equivalent,    has  often  been  used  in  past  allocation  studies,    this  usage 
is  deprecated  since  it  does  not  provide  a  solution  which  is  as  well  adapted  to  the  actual  nature 
of  the  problem  as  the  separation  of  the  fading  into  its  two  components      Y.     and     Y     and  the 
separate  use  of    (V.14)    and    (V.18).     Note  that    (V.19)    provides  a  fading  allowance  which  is 
too  small  compared  with  that  estimated  using    (V.14)    and    (V.18)    separately.      The  recom- 
mended method  makes  more  appropriate  allowance  for  the  fact  that  communications  at 
particular  times  of  the  day  or  for  particular  seasons  of  the  year  are  more  difficult  than  at 
other  times. 
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V.  5    The  Joint  Effect  of  Several  Sources  of  Interference  Present  Simultaneously 

The  effects  of  interference  from  unwanted  signals  and  from  noise  have  so  far  been 

considered  in  this  report  as  though  each  affected  the  fidelity  of  reception  of  the  wanted  signal 

independently.      Let    p        (g)     and    r      (g)p  denote  power  levels  which  the  wanted  median 

mr  ur  um 

signal  power    p        inust  exceed  in  order  to  achieve  a  specified  grade  of  service  when  each  of 
these  sources  of  interference  is  present  alone.      To  the  extent  that  the  various  sources  of 
interference  have  a  character  approximating  that  of  white  noise,    this  same  grade  of  service 
may  be  expected  from  a  v/anted  signal  with  median  level 

P       =  P        (g)  +      /     !■      (r)  P 
^m       ^mr^^'         /^      ur^^'  ^um 

when  these  sources  are  present  simultaneously. 

An  approximate  method  has  been  developed    [Norton,    Staras,    and  Blum,    1952  ]     for 
determining  for  a  broadcasting  service  the  distribution  with  time  and  receiving  location  of  the 
ratio 

Pm^Pmr^g)  +  I  "ur^^^  ^um 

Although  this  approach  to  the  problem  of  adding  the  effects  of  interference  will  probably 

always  provide  a  good  upper  bound  to  the  interference,    this  assumption  that  the  interference 

power  is  additive    is  often  not  strictly  valid.      For  example,   when  intelligible  cross -talk  from 

another  channel  is  present  in  the  receiver  output  circuit,    the  addition  of  some  white  noise 

will  actually  reduce  the  nuisance  value  of  this  cross -talk. 

Frequently,    however,   both     p        (g)     and     p  will  be  found  to  vary  more  or  less 

mr  °  um 

independently  over  wide  ranges  with  time  and  a  good  approximation  to  the  percentage  of  time 

that  objectionable  interference  is  present  at  a  particiilar  receiving  location  may  then  be 

obtained    [Barsis,    et  al,    1961  ]    by  adding  the  percentage  of  time  that     p         is  less  than 

p       (g)     to  the  percentages  of  time  that     p         is  less  than  each  of  the  values  of     r      (g)  p 

When  this  total  time  of  interference  is  small,    say  less  than     10%,   this  will  represent  a 

satisfactory  estimate  of  the  joint  influence  of  several  sources  of  interference  which  are 

present  simultaneously.      Thus,    when  the  fading  ranges  of  the  various  sources  of  interference 

are  sufficiently  large  so  that  this  latter  method  of  analysis  is  applicable,   the  various  values 

of     p       (g)     and  of     r      (g)  p  will  have  comparable  magnitudes  for  negligible  percentages 

of  the  time  so  that  one  may,    in  effect,    assume  that  the  various  sources  of  interference  occur 

essentially  independently  in  time. 

Minimum  acceptable  wanted-to-unwanted  signal  ratios      r  may  sometimes  be  a 

function  of     r      ,    the  available  wanted  sienal-to-noise  ratio.      When     r  is  within     3  db    of 

m  ur 

r        ,    an  unwanted  signal  may  be  treated  the  same  as  external  noise,    and,    in  a  similar  fashion, 

long-term  distributions  of  available  wanted-to-unwanted  signal  ratios  may  be  determined  for 

each  class  of  unwanted  signals  for  which     r  is  nearly  the  same. 

"  ur 
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V.  6     The  System  Equation  for  Noise -Limited  Service 

Essential  elements  of  a  noise-limited  communication  circuit  are  summarized  in  the 

following  system  equation.      The  transmitter  output     P    dbw    which  will  provide     P    =  P      -  L 

dbw    of  total  radiated  power  in  the  presence  of  transmission  line  and  matching  network  losses 

L      db,    and  which  will  provide  a  median  delivered  signal  at  the  pre-detection  receiver  output 

which  is    R       db   above  the  median  noise  power    P  delivered  to  the  pre-detection  receiver 

m  mn 

output  is  given  by 

(P„    -  L„   )  =  L       +R       +  (P  -G        )      dbw  (V.20) 

it        it  m  m  mn  ms 

in  the  presence  of  a  median  transmission  loss    L        and  a  median  operating  receiving  system 

m 

signal  gain  G        .     The  operating  signal  gain  is  the  ratio  of  the  power  delivered  to  the  pre- 
ms 

detection  receiver  output  to  the  power  available  at  the  terminals  of  an  equivalent  loss-free 
antenna.     G      is  its  maximunn  value  for  c.w.   frequencies  in  the  receiver  pass  band,    and   G 

O  -1  r-  j^g 

is  its  median  value  for  all  signal  frequencies  in  the  pass  band.     (P         -G       )    in  (V.20)  is  the 

mn         ms 

equivalent  median  noise  power  at  the  antenna  terniinals,    as  defined  in  section  V.3. 

It  is  appropriate  to  express  the  system  equation    {V.20)    in  terms  of  the  operating 

noise  factor     F  defined  by    (V.8),    rather  than  in  terms  of     P  or     (P        -  G       )     in 

op  mn  mn         ms 

order  to  separate  studies  of  receiving  system  characteristics  from  studies  of  propagation. 

For  this  reason  all  predicted  power  levels  are  referred  to  the  terminals  of  an  equivalent 

loss -free  antenna,    and  receiving  system  characteristics  such  as     F      ,     G   ,     G       ,     and 

op  o  ms 

B  =   10  log   b       are  separated  from  transmission  loss  and  available  power  in  the  formulas. 

Rearranging  terms  of    (V.8),    the  equivalent  median  noise  power     P        -  G  delivered 

mn         ms 

to  the  antenna  terminals  may  be  expressed  as 

P         -G         =F       +(G-G       )  +  (B  -  204)  (V.  21) 

mn         ms  op  o       ms 

where     G       and     G  are  usually  nearly  equal.     Assuming  that     Li„    ,     G   ,     G       ,     and     B 

are  constant,    it  is  convenient  to  combine  these  parameters  into  an  arbitrary  constant     K    : 

K     =L,.     +G     -G         +  B  -  204        dbw  (V.  22) 

o         it  o  ms  ■      ' 

and  then  write  the  system  equation  as  : 

P„,=  K+L      +R       +F  dbw  (V.  23) 

i  t        o  m         m  op 

In  general,    if  unwanted  signals  other  than  noise  may  be  disregarded,    service  exists 

whenever     R     (p)     exceeds     R       (g),   where     R     (p)     is  the  value  of     R         exceeded    p    percent 
m  mr  rn  m 

of  all  hours.    With     G  and     P  assumed  constant,    so  that 

ms  mn 

R     (p)  =  P     (p)  +  G         -  P  (V.  24) 

m  m  ms  mn 


y-u 


service  exists  whenever     P     (p)     exceeds     P        ,    or  whenever     L     (p)     is  less  than  the  max- 

m  mr  m^^' 

imum  allowable  transmission  loss     L       (g)  =  P     -  P       (e)  .     An  equivalent  statement  may  be 

mio  °  t  mr "  ' 

made  in  terms  of  the  system  equation.      The  transmitter  power    P     (p)    which  will  provide 

for    p    percent  of  allhours  at  least  the  grade     g    service  defined  by  the  required  signal-to-noise 

ratio     R       (g)     is 


mr 


^it(P)  =  K^  +  LJp)+F       +R^^(g)  (V.25) 


op 

where     L     (p)     is  the  hourly  median  transmission  loss  exceeded    (100 -p)    percent  of  all 
hours, or  not  exceeded    p    percent  of  all  hours. 

For  a  fixed  transmitter  power     P       dbw,    the  signal-to-noise  ratio  exceeded    p    percent 
of  all  hours  is 

R     (p)  =  P     -  K     -  F       -  L     (p)        db  (V.26) 

m  o  o  op  m 

for  an  "average"  or  "median"  propagation  path. 

The  maximum  allowable  transmission  loss 

L       (g)  =  P     -  K     -  F       -  R       (g)  (V.27) 

mo  o  o  op  mr  °  ^ 

is  set  equal  to  the  loss      L     (P)  Q)     exceeded    (100 -p)    percent  of  all  hours  with  a  probability 

Q.     This  value  is  fixed  when     P    ,     K   ,     and     R       (g)     have  been  determined,   and  for  each 

o  o  mr 

time  availability    p    there  is  a  corresponding  service  probability,      Q(p)-     Section    V.  9  will 
explain     how  to  calculate     Q(p). 

When  external  noise  is  both  variable  and  not  negligible,   the  long-term  variability  of 
F         must  be  considered,   and  the  following  relationships  may  be  used  to  satisfy  the  condition 

R     (P)>R     ^(g)  {V.28) 

m  mr 

R     (p)  =  R^(50)  +  Y     (p)  (V.  29) 

mm  m 

R     (50)^    P     -K     -F     (50)  -  L     (50)  (V.30) 

m  o  o  op  rn  *  ' 

Y    ^p)  ^    Y^p)  +  Y^(lOO-p)  -  2p       Y(p)  Y   (100-p)  (V.31) 

m  n  tn  n 

Y(p)  H    L^(50)  -  LJp)    .         Yjp)  H    F^p(p)  -  F^p(50)  (V.32) 

where      p         is  the  long-term  correlation  between     P         and     F      .      Though     p  could 

■^tn  *  m  op  "        "^tn 

theoretically  have  any  value  between    -  1    and    1,   it  is  usually  zero. 
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V.  7     The  Time  Availability  of  Interference -Limited  Service 

Let     p         denote  the  lone-term  correlation  between     P         and     P       ,    the  power  ex- 
■^tu  m  um 

pected  to  be  available  at  least    p    percent  of  all  hours  at  the  terminals  of  an  equivalent  loss- 
free  receiving 
respectively : 


free  receiving  antenna  from  wanted  and  unwanted  stations  radiating     p        and     p       watts, 


P     (p)  =  P     -L     (p)      dbw   ,         P       (p)  =  P     -  L       (p)      dbw  (V.33) 

m  o  m  um  u  um  ^  ' 

P     =   10  log  p       dbw  ,         P     =   10  log  p        dbw  (V.34) 

o  o  u  u  '  ' 

The  criterion  for  service  of  at  least  grade    g    in  the  presence  of  a  single  unwanted  signal  and 
in  the  absence  of  other  \anwanted  signals  or  appreciable  noise  is 

R^(P)  >  R^^{g,  P)  {V.3  5) 

■where 

R   (p)  =  R   (50)  +  Y    (p)  (V.36) 

U  u.  xv 

R   (50)  ^    P     (50)  -  P       (50)  (V.37) 

u  m  um 

Y„^P)^    Y^(p)  +  Y^(lOO-p)  -  2p       Y(p)Y(100-p)  (V.  38) 

x\,  u  lu.  u 

Y   (p)  =  P       (P)  -  P       (50)  =  L       (50)  -  L       (p)  (V.39) 

u  um  um  um  um 

If     P     ,     P       ,     and     F         were  exactly  normally  distributed,     (V,31)    and    (,V.38)    would  be 
m  um  op 

exact,   and  they  represent  excellent  approximations  in  practice. 


V-20 


V,8     The  Estimation  of  Prediction  Errors 

Consider  the  calculation  of  the  power     P     (p)     available  at  the  terminals  of  an  equiva- 

rn 

lent  loss-free  receiving  antenna  during    p    percent  of  all  hours.     P     (p)     refers  to  hourly  median 
values  expressed  in   dbw  .     For  a  specific  propagation  path  it  is  calculated  in  accordance  with 
the  methods  given  in  sections  2-10  using  a  given  set  of  path  parameters     (d,  f,  0,  h     ,    etc.).     De- 
note by    P        (p)     observations  made  over  a  large  number  of  randomly  different  propagation  paths, 

mo  i.        J.     o 

which, however , can  all  be  characterized  by  the  same  set  of  prediction  parameters.     Values  of 

P        (p)     will  be  very  nearly  normally  distributed  with  a  mean  (and  median)   equal  to    P     (p), 
mo  -,  m 

and  a  variance  denoted  by     a   (p).      This  path-to-path  variability  is  illustrated  in  Fig.    V.4  for 
a  hypothetical  situation  which  assumes  a  random  distribution  of  all  parameters  which  are  not 

taken  into  account  in  the  prediction  method. 

2 
The  variance     cr       of  deviations  of  observation  from  prediction  depends  on  available 

data  and  the  prediction  method  itself.  The  most  sophisticated  of  the  methods  given  in  this  re- 
port for  predicting  transmission  loss  as  a  function  of  carrier  frequency,  climate,  time  block, 
antenna  gains,  and  path  geometry  have  been  adjusted  to  show  no  bias,  on  the  average,  for  the 
data  discussed  in  subsections   10.3  and  10.4  and  in  Annex  I. 

Most  of  these  data  are  concentrated  in  the  40-1000  MHz  frequency  range,    and  were 
obtained  primarily  for  transhorizon  paths  in  climates   1,    2,    and  3.     Normally,    one  antenna  was 
on  the  order  of  10  meters  above  ground  and  the  other  one  was  higher,    near  200  meters.     Even 
the  low  receiving  antennas  were  either  on  high  ground  or  in  clear  areas  well  removed  from 
hills  and  terrain  clutter.     Few  of  the  data  were  obtained  with  narrow-beam  antennas.     An  at- 
tempt has  been  made  to  estimate  cumulative  distributions  of  hourly  transmission  loss  medians 
for  accurately  specified  time  blocks,    including  estimates  of  year-to-year  variability. 

A  prediction  for  some  situation  that  is  adequately  characterized  by  the  prediction 

parameters  chosen  here  requires  only  interpolation  between  values  of  these  parameters  for 

which  data  are  available.     In  such  a  case,      cr   (p)     represents  the  standard  error  of  prediction. 

c 

The  mean  square  error  of  prediction,    referred  to  a  situation  for  which  data  are  not  available,  is 

2 
cr   (p)     plus  the  square  of  the  bias  of  the  prediction  method  relative  to  the  new  situation. 

Based  on  an  analysis  of  presently  available  transhorizon  transmission  loss  data,    the 

2 
variance     cr   (p)     is  estimated  as 

o-^(p)     =   12.73  +  0.  12  Y^(p)      db  (V.40) 

c 

2 

where     Y(p)     is  defined  in  section  10,     Since     Y(50)  =  0,    the  variance     cr   (50)     of  the  difference 

2  "^ 

between  observed  and  predicted  long-term  medians  is     12.73  db   ,    with  a  corr ssponding  stand- 
ard deviation     cr   (50)  =  3.  57  db. 
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It  is  occasionally  very  difficult  to  estimate  the  prediction  error     cr   (p)     and  the 
service  probability    Q  .     Where  only  a  small  amount  of  data  is  available  there  is  no  adequate 
way  of  estimating  the  bias  of  a  prediction.     One  may,    however,    assign  weights  to  the  curves 
of    V(50,    d^)     in  figure  10.  1  for  climates   1-7  based  on  the  amount  of  supporting  data  avail- 
able: 


Cli 

mate  Number 

Weight 

1 

300 

2 

120 

3 

60 

4 

2 

5 

1 

6 

5 

7 

5 

As  an  example,      for    d     =  600  km,    the  average    ¥(50)     weighted  in  accordance  with  the  above 

®  2 

is     -0.4  db,    and  the  corresponding  climate-to-climate  variance  of    V(50)     is    3.7  db    . 

2 
If  the  variance  of    V(50)     is   12.7  db      within  any  given  climate,    the  standard  error 


of  prediction    L(50)  =  L       +  0 . 4  db,    given  by  ( 10.  4),    will  be         ^712.  7  +  3,7  =  4.0  db    if  a 

cr 

random  sarapling  of  these  climates  is  desired  or  applicable.     If  there  is  doubt  as  to  which  of 

two  climates    i    and    j     should  be  chosen,    the  best  prediction  of    L,  (p)     might  depend  on  the 

average  of    V.(50,    d  )     and    ¥.{50,    d  )       and  the  root-mean  square  of    Y.(p,    d  )     and    Y.(p,    d  ); 
1  e  J  e  X  e  J  e 


L(p)  =  L        -0.5 
cr 


V.(50,    d  )  +  V.(50,    d  )      -  Y(p,    d  )       db, 
1  e  J  e  e 


Y(p,d  )  = 


0.  5  Y^(p,    d  )  +  0,  5  Y^(p,    d  ) 
1  e  J  e 


db 


The  bias  of  this  prediction  may  be  as  large  as    0.5 


V.{50,    d   )   -  V.(50,    d  ) 
1  e  J  e 


(V.41) 


(V.42) 


db .      The 


root-nnean  square  prediction  error  may  therefore  be  estimated  as  the  square  root  of  the  sum 

2  2  2 

of    0-   (50)  =   12.7  db    ,    the  square  of  the  bias,    and    0.12  Y(p,    d  )     ,    or: 


{ 


12.7  +  0,25 


1 

V.(50,    d  )   -  V.(50,    d  )   I     +0.12  Y^(p,    d  )^^ 
1  e  J  e     I  e 


db. 


with    Y(p,    d  )     obtained  from  (V.  42)  above. 

According  to  figure  10.  1,     V(50,    d  )     is  expected  to  be  the  same  for  climates   1  and 

8.      This  conclusion  and  the  estimate  for     Y(p,    d  )     shown  in  figure  III.  55  for  climate  8  are 

e 

based  solely  on  meteorological  data.     In  order  to  obtain  these  estimates,    the  percentages  of 

time  for  which  surface-based  ducts  existed  in  the  two  regions  were  matched  with  the  same 

value  of    Y(p,    d  )     for  both  climates.     In  this  way,      Y„(p,    d  )     was  derived  from    Y    (p,    d  ) 
e  8  e  r^        e 

by  relating    p       to    p       for  a  given    Y    instead  of  relating     Y       to    Y       for  a  given    p. 
o  i  8  1 
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V.9     The  Calculation  of  Service  Probability    Q    for  a  Given  Time  Availability    p 
For  noise -limited  service  of  at  least  grade    g    and  a  time  availability    p,    the  ser- 
vice probability    Q    is  the  probability  that 

L       (g)    -  L     (p)   >0  (V.43) 

mo  m 

if  external  noise  is  negligible.     L       (g)     is  defined  by    (V.27)  .       The  criterion  for  service 

limited  by  external  noise,   where     F         as  given  by    (V.8)    is  variable,   is 

R     (p)  -  R       (g)  >  0    (from  equation  V.  28) 
m  mr 

For  service  limited  only  by  interference  from  a  single  unwanted  signal, 

R   (P)  -  R      (gx  P)  >  0     (from  equation  V.35) 
u  ur 

Combining    (V.22)        and    (V.27),     (V.43)    may  be  rewritten  as 

P     -  L„     -  G     +  G         -  B  +  204  -  F        -R       (g)  -  L     (p)  >  0  (V.  44) 

o  it  o  ms  op  mr^*'  m^^'  ^ 

where     (V.8),     (V.9),     or  section    V.6    may  be  referred  to  for  a  definition  of  the  terms  in 

(V.44).       Assuming  that  the  error  of  estimation  of  these  terms  from  system  to  system  is 

2 
negligible  except  for  the  path-to-path  variance     o-    (p)     given  by    (V.40)     for  the  last  term, 

it  is  convenient  to  represent  the  service  probability     Q     as  a  fiinction  of  the  standard  normal 

deviate 

z         =  j-^ (V.  45) 

mo  or  (p) 

c 

which  has  a  mean  of  zero  and  a  variance  of  unity.      L  is  identified  as  the  transmission 

mo 
loss  exceeded    (100  -p)     percent  of  the  time  with  a  probability     Q,    which  is  expressed  in 

terms  of  the  error  function     erf  x    as 

Q(z       )  =1+1  erf (z        /'TZ)  (V.  46) 

mo         '^        '^  mo 

Figure    V.5    is  a  graph  of     Q     versus      z 

mo 

For  the  method  described  here,    the  condition 

0.12  Y(p)  z        (Q)  <    -0-  (p)  (V.  47) 

mo  c 

is  sufficient  to  insure  that  the  service  probability     Q     increases  as  the  time  availability     p 
is  decreased.     A  less  restrictive  condition  is 

Y(p)[L         -L     (50)]   <   106      db  (V.  48) 

mo  m 
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An  example  is   shown  in  figure  V.6,    with    p    versus     Q     for  radiated  powers    P    =30   dbw 

o 

and    P     =  40  dbw,    and    L     (p,  Q)  =  P     +   140  db.     Here, 
o  mo 


50  +  50  erf 


L     (p)  -   140 


lO-sTT 


(V.49) 


corresponding  to  a  normal  distribution  with  a  mean    L     (50)  =   140  db    and  a  standard 

m 

deviation    Y(15.8)  =   10  db.     [Note  that       L     (p)     versus     p    as  estimated  by  the  methods  of 
section  10  is  usually  not  normally  distributed]. 

To  obtain  the  time  availability  versus  service  probability  curves  of  figure    Y.b,     L     (p) 

2  "^ 

was  obtained  from     p,     Y(p)     from    {V.4),      o"    (p)     from    (V.40),     Zrno   from    (V.45),     and     Q 

c 

from  figure    V.5.      This  same  method  of  calculation  may  be  used  when  there  are  additional 

2 
sources  of  prediction  error  by  adding  variances  to     or    (p).     Examining  possible  trade-offs 

between  time  availability  and  service  probability  shown  in  figure    V.5,    note  the  increase 

from    p  =  96.5%     to     p  =  99.3%     for     Q  =  0.95,    or  the  increase  from     Q  =  0.78    to     Q  =  0.97 

for     p  =  99%,    as  the  radiated  power  is  increased  from  one  to  ten  kilowatts. 

For  the  case  of  service  limited  by  external  noise  (V.  28)   to  (V.  30)  may  be  rewritten 


P     -  K      -  F      (50)    -  L     (50)  +  Y     (p)    -  R        (g)   >  0 
o  o  op       m  m  mr 


(V,  50) 


One  may  ignore  any  error  of  estimation  of    P    ,    K    ,    and    R        (g)     as  negligible  and  assume  no 
•'     "  '  o        o  mr  ^ 

path-to-path  correlation  between    F      (50)     and    L     (50).      The  variance     a      (p)     of    F      (50)  + 


op 


op 


op  J-ii-  ^t"  o 

L     (50)    -  Y     (p)     in  (V.48)   may  then  be  written  as  a  sum  of  component  variances     tr        and 


m 


cr^    (p)  =  0-^    +   12.  73  +  0.  12  Y^  (p)       db^ 
op  F  m 


(V.51) 


Very  little  is  known  about  values  for  the  variance     cr  „    of    F      (50),    but  it  is  probably  on  the 

'  i^  op 

2 
order  of    20  db    . 

The  corresponding  standard  normal  deviate  is  given  by 


op 


R    (P)  -  R      (g) 

m  mr 

op 


(V.52) 


and  the  service  probability    Q(p)     is  given  by  (V.  46)  with    z^^    replaced  by    z  The  restric- 

tion (V    47)   still  holds  with    z  and    cr       replaced  by    z         and    cr       .     A  less  restrictive  con- 

'  mo  c  op  op 

dition  equivalent  to  (V.48)    can  be  stated  only  if  a  specific  value  of    cr         is  assumed. 
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For  the  case  of  service  limited  only  by  interference  from  a  single  unwanted  radio 
signal  (V.  35)   to  (V.  39)   may  be  rewritten  as 

L       (50)    -  L     (50)  +  Y    (p)    -  R      (g,  p)   >  0  {V.  53) 

um  m  R  ur  ^  ' 

Let    p         denote  the  normalized  correlation  or  covariance  between  path-to-path  variations  of 

^"^  2  2 

P     (50)     and    P       (50).      Then  assuming  a  variance  of    25.  5(  1  -  n ,  )  +  0 .  12  Y    (p)   db       for  the 
m  um  ^iu  u 

first  three  terms  of  (V.  53)   and  a  variance     cr  for  the  estimate  of    R      (g,  p)  ,    the  total  var- 

,  ur  ur 

2 

iance     cr       (p)     of  any  estimate  of  the  service  criterion  given  by  (V,53)    rnay  be  written  as 
uc 

a-J(p)  =  25.5(1 -p^^)  +  0.12  Y^^p)  +  ^J  (V.  54) 

for  this  case,    corresponding  to  a  standard  normal  deviate 

R^(P)  -  R,,(P.  g) 

z        =   7-r fV.  55) 

and  a  service  probability     Q(p)     given  by    (V.46)    with     z  replaced  by     z      .      The  variance 

0"  may  range  from     10    db       to  very  much  higher  values.      The  restrictions    (V.47)     and 

ur  ° 

(V.48)    apply  with     z  and     o-       replaced  by     z  and     cr  and  with     106    db       in  (V,48) 

■^^  ^  mo  c  ^  ^        uc  uc  \  I 

replaced  by    (212  +  o-    ^)    db    . 
'  ur 
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TRANSMISSION    LOSS   L     IN  DECIBELS 
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V,  10    Optimum  Use  of  the  Radio  Frequency  Spectrum 

The  business  of  the  radio  engineer  is  to  develop  efficient  radio  systems,    and  the 
principal  tool  for  improving  efficiency  is  to  adjust  the  various   parameters  to  their  optimum 
values.      For  example,    it  is  usually  more  economical  to  use  lower  effective  radiated  powers 
from  the  transmitting  systems  by  reducing  the  operating  sensitivities  of  the  receiving 
systems.     Receiving  system  sensitivities  can  be  reduced  by    (a)  reducing  the  level  of  inter- 
nally generated  noise,     (b)  the  use  of  antenna  directivity  to  reduce  the  effects  of  external 
noise,     (c)  the  reduction  of  man-made  noise  levels  by  the  use  of  suppressors  on  noise  gene- 
rators such  as  ignition  systems,    relays,    power  transmission  systems,    etc.  ,   and    (d)  the  use 
of  space  or  time  diversity  and  coding.      The  use  of  more  spectrum  in  a  wide  band    FM    system 
or  in  a  frequency  diversity  system  can  also  reduce  the  receiving  system  operating  sensitivity 
as  well  as  reduce  the  acceptance  ratios  against  unwanted  signals  other  than  noise. 

Unfortunately,   lonlike  other  natural  resources  such  as  land,   minerals,    oil,    and 
water,  there  is  currently  no  valid  method  for  placing  a  monetary  value  on  each  hertz  of  the 
radio  spectrum.      Thus,    in  the  absence  of  a  common  lonit  of  exchange,    these  tradeoffs  are 
often  made  unrealistically  at  the  present  time.     It  is  now  generally  recognized  that  the  use  of 
large  capacity  computers  is  essential  for  optimizing  the  assignment  of  frequencies  to  various 
classes  of  service  including  the  development  of  optimum  channelization  schemes.     Typical 
inputs  to  such  computers  are: 

1.  Nominal  frequency  assignments. 

2.  Transmitting  system  locations,   including  the  antenna  heights. 

3.  Transmitting  system  signatures;    i.e.  ,   the  radiated  emission  spectrum  charac- 
teristics including  any  spurious  emission  spectrums. 

4.  Transmitting  antenna  characteristics. 

5.  Receiving  system  locations,   including  the  antenna  heights. 

6.  Spurious  emission  spectrums  of  the  receiving  systems. 

7.  Operating  sensitivities  of  the  receiving  systems  in  their  actual  environments 
which  thus  make  appropriate  allowance  for  the  effects  of  both  man-made  and 
natural  noise. 

8.  Required  values  of  wanted-to -unwanted  phase  interference  median  signal  powers 
for  all  unwanted  signals  which  could  potentially  cause  harmful  interference  to 
the  wanted  signal;    these  acceptance  ratios  include  appropriate  allowances  for 
reductions  in  the  effects  of  fading  achieved  by  the  use  of  diversity  reception  and 
coding. 
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9.  Long-term  median  reference  values  of  basic  transmission  loss  and  path  antenna 
gain  for  the  wanted  path  and  all  of  the  lonwanted  signal  propagation  paths;  these 
path  antenxia  gains  include  allowances  for  antenna  orientation,    polarization,    and 
multipath  phase  mismatch  coupling  losses. 

10.  Distributions  with  time  of  the  transmission  loss  for  the  wanted  signal  path  and 
all  of  the  unwanted  signal  paths. 

11.  Correlations  between  the  transmission  losses  on  the  wanted  and  on  each  of  the 
unwanted  propagation  paths. 

12.  Transmission  line  and  antenna  circuit  losses. 

13.  The  spurious  emission  spectrum  of  any  unwanted  signals  arising  from  unli- 
censed sources  such  as  diathermy  machines,    electronic  heaters,   welders, 
garage  door  openers,    etc. 

14.  Assigned  hours  of  operation  of  each  wanted  and  each  unwanted  emission. 

The  output  of  the  computer  indicates  simply  the  identity  and  nature  of  the  cases  of  harmful 
interference  encountered.     Harmful  interference  is  defined  as  a  failure  to  achieve  the  speci- 
fied grade  of  service  for  more  than  the  required  percentage  of  time  during  the  assigned  hours 
of  operation.     Changing  some  of  the  inputs  to  the  computer,   an  iterative  process  can  be 
defined  which  may  lead  to  an  assignment  plan  with  no  cases  of  harmful  interference. 

It  is  assumed  that  a  given  band  of  radio  frequencies  has  been  assigned  to  the  kind  of 
radio  service  \inder  consideration  and  that  the  nature  of  the  services  occupying  the  adjacent 
frequency  bands  is  also  kjiown.      Furthermore,    it  is  assumed  that  the  geographical  locations 
of  each  of  the  transmitting  and  receiving  antennas  are  specified  in  advance,   together  with  the 
relative  values  of  the  radiated  powers  from  each  transmitting  antenna  and  the  widths  and 
spacings  of  the  radio  frequency  channels.     In  the  case  of  a  broadcasting  service  the  speci- 
fication of  the  intended  receiving  locations  can  be  in  terms  of  proposed  service  areas.      With 
this  information  given,    use  may  be  made  of  the  following  procedures  in  order  to  achieve 
optimum  use  of  the  spectrum  by  this  particular  service  : 

(a)     The  system  loss  for  each  of  the  wanted  signal  propagation  paths  should  be  mini- 
mized and  for  each  of  the  unwanted  signal  propagation  paths  should  be  maxinnized; 
this  may  be  accomplished  by  maximizing  the  path  antenna  power  gains  for  each  of  the 
wanted  signal  propagation  paths,    minimizing  the  path  antenna  gains  for  each  of  the 
unwanted  propagation  paths,    and  in  exceptional  cases,   by  appropriate  antenna  siting. 
The  path  antenna  gains  for  the  imwanted  signal  propagation  paths  may  be  minimized 


V-30 


by  the  use  of  high-gain  transmitting  and  receiving  antennas  with  optimum  side  lobe 
suppression  and  front-to-back  ratios  and,    in  some  cases,   by  the  use  of  alternate 
polarizations  for  geographically  adjacent  stations  or  by  appropriate  shielding. 

(b)     The  required  protection  ratios     r      (g)     should  be  minimized  by    (1)    appropriate 
radio  system  design,     (2)    the  use  of  stable  transmitting  and  receiving  oscillators, 
(3)    the  use  of  linear  transmitting  and  receiving  equipment,     (4)    the  use  of  wanted 
and  unwanted  signal  propagation  paths  having  the  minimuna  practicable  phase  inter- 
ference fading  ranges;  from  band  6  to  band  9  {0,3  to  3000  MHz),    minimunn  phase  in- 
terference fading  may  be  achieved  by  the  use  of  the  maximum  practicable  transmitting 
and  receiving  antenna  heights,    and  (5)   the  use  of  space  diversity,    time  diversity,    and 
coding. 

(c)  Wanted  signal  propagation  paths  should  be  employed  having  the  minimunn  practi- 
cable long-term  power  fading  ranges.     In  bands  8  and  9,    minimum  fading  may  be 
achieved  by  the  use  of  the  maximum  practicable  transmitting  and  receiving  antenna 
heights. 

The  above  procedures  should  be  carried  out  with  various  choices  of  transmitting  and  receiv- 
ing locations,   relative  transmitter  powers,   and  channel  spacings  until  a  plan  is  developed 
which  provides  the  required  service  with  a  minimum  total  spectriim  usage.     After  the  un- 
wanted signal  interference  has  been  suppressed  to  the  maximum  practicable  extent  by  the 
above  methods  so  that,   at  each  receiving  location  each  of  the  values  of     r       exceeds  the 
corresponding  protection  ratio     r      (g)      for  a  sufficiently  large  percentage  of  the  time,   then 
the  following  additional  procedures  should  be  adopted  in  order  to  essentially  eliminate  inter- 
ference from  noise  : 

(d)  The  system  loss  on  each  of  the  wanted  signal  propagation  paths  should  be  mini- 
mized;   this  maybe  accomplished  by    (1)    the  use  of  the  highest  practicable  transmit- 
ting and  receiving  antenna  heights  in  bands  8  and  9,    and    (2)    maximizing  the  path 
antenna  power  gains  for  each  of  the  wanted  signal  propagation  paths.     The  path  antenna 
power  gains  of  the  wanted  signal  propagation  paths  may  be  maximized  by  using  the 
maximum  practicable  transmitting  and  receiving  antenna  gains  and  by  minimizing  the 
antenna  circuit  and  polarization  coupling  losses.     The  minimization  of  the  system 
loss  on  each  of  the  wanted  signal  propagation  paths  will  already  have  been  achieved  to 
a  large  extent  in  connection  with  procedures  (a),    (b),    and  (c)  above. 
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(e)  In  general,    receiving  systems  should  be  employed  which  have  the  lowest  practi- 
cable values  of  operating  sensitivity     p        (g)  . 

(f)  Finally,    sufficiently  high  transmitter  powers  should  be  used    [keeping  the  relative 
powers  at  the  optimum  relative  values  determined  by  procedures  (a),    (b),    and  (c)  ] 

so  that  the  wanted  signal  power     p         will  exceed  the  operating  sensitivity     p       (g) 
for  a  sufficiently  large  percentage  of  the  time  during  the  intended  period  of  operation 
at  every  receiving  location. 

Although  it  might  at  first  seem  impracticable,    serious  consideration  should  be  given  to  the 
use  of  auxiliary  channels  from  wanted  receivers  to  wanted  transmitters.      The  provision  of 
such  channels  might  well  be  feasible  in  those  cases  where  two-way  transmissions  are  in- 
volved and  might  lead  to  important  economies  in  both  power  and  spectrum  occupancy 
[Hitchcock  and  Morris,    1961  ]  . 

Ultimately,   when  optimum  use  of  the  spectrum  has  been  achieved,    it  will  not  be  possi- 
ble to  find  a  single  receiving  location  at  which  radio  noise  rather  than  either  wanted  or 
unwanted  signals  can  be  observed  for  a  large  percentage  of  the  time  throughout  the  usable 
portions  of  the  radio  spectrum  not  devoted  to  the  study  of  radio  noise  sources,    as  is  the  radio 
astronomy  service.     Although  everyone  will  agree  that  the  attainment  of  this  ideal  goal  of 
interference -free  spectrum  usage  by  the  maximum  number  of  simultaneous  users  can  be 
achieved  only  over  a  very  long  period  of  time  because  of  the  large  investments  in  radio 
systems  currently  in  operation,   nevertheless  it  seems  desirable  to  have  a  clear  statement  of 
the  procedures  which  should  be  employed  in  the  future  in  order  to  move  in  the  direction  of 
meeting  this  ultimately  desirable  goal  whenever  appropriate  opportunities  arise. 
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